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wave  equation  with  nonlinear  boundary  conditions.  In  addition  to  ordinary 
rigid-body  scattering.  Censor  predicted  nongrowing  waves  at  frequencies  equal 
to  the  sum  and  to  the  difference  of  the  frequencies  of  the  primary  waves.  The 
solution  to  the  nonlinear  wave  equation  also  yields  scattered  waves  at  the  sum 
and  difference  frequencies.  However,  the  nonlinearity  of  the  medium  causes 
these  waves  to  grow  with  increasing  distance  from  the  scatterer's  surface  and, 
after  a  very  small  distance,  dominate  those  predicted  by  Censor. -s 

■  '  y 

(j 

The  simple-source  formulation  of  the  second-order  nonlinear  wave  equation  for 
a  lossless  fluid  medium  has  been  derived  for  arbitrary  primary  wave  fields. 
(Westevelt's  original  derivation  in  J.  Acoust.  Soc.  Am.  35,  535-537,  1963)>*of 
this  particular ~  form  of  ther  HSHtiirear  "wave,  gquat ion  was  restricted  ^te^'pTane 
waves.  J4”* This  equation  was  used  to  solve  the  problem  of  nonlinear  scattering  of 
acoustic  waves  by  a  vibrating  obstacle  for  three  geometries:  1)  a  plane-wave 
scattering  by  a  vibrating  plane,  2)  cylindrical-wave  scattering  by  a  vibrat¬ 
ing  cylinder,  and  3)  plane-wave  scattering  by  a  vibrating  cylinder..  A  new 
technique  of  integration  useful  for  solving  definite  integrals  arising  from 
physical  problems  is  developed. 

/ 

^Successful  experimental  validation  of  the  theory  was  inhibited  by  previously 
unexpected  levels  of  nonlinearity  in  the  hydrophones  used.  Such  high  levels 
of  hydrophone  nonlinearity  appeared  in  hydrophones  that,  by  their  geometry  of 
construction,  were  expected  to  be  fairly  linear.  It  appears  that  this  is  a 
rather  general  problem  with  hydrophones.  A  new  technique  for  measuring  this 
hydrophone  nonlinearity  is  presented.  /  \ 
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1.  INTRODUCTION 


This  thesis  addresses  the  problem  of  the  nonlinear  scattering  of  acoustic 
waves  with  harmonic  time  dependence  by  a  vibrating  obstacle.  The  obstacle  is 
immersed  in  an  infinite  homogeneous  fluid  medium,  and  its  surface  deforms 
uniformly  with  harmonic  time  dependence.  The  case  in  which  the  surrounding 
fluid  medium  is  water  is  of  primary  interest. 

The  problem  of  the  scattering  of  a  plane  wave  incident  on  acoustically 
rigid  spherical  and  cylindrical  obstacles  surrounded  by  an  infinite  homoge¬ 
neous  fluid  medium  was  first  solved  by  Rayleigh  [1]  using  linear,  lossless 
theory.  To  obtain  numerical  values  for  the  rather  complicated  mathematical 
solutions,  he  considered  only  the  limiting  case  where  the  radius  of  the  scat- 
terer  is  much  smaller  than  a  wavelength.  These  results  were  extended  by  Morse 
[2],  who  obtained  a  solution  in  a  form  more  readily  evaluated.  He  provided 
tables  of  calculated  values,  Including  values  when  the  obstacle  is  not  small 
compared  to  a  wavelength.  Faran  [3]  obtained  a  solution  in  the  case  where  the 
scattering  obstacles  are  not  acoustically  rigid.  In  addition  to  studying  the 
problem  theoretically,  Faran  performed  an  experimental  investigation  of  the 
cylindrical  case.  Agreement  between  theory  and  experiment  was  excellent  and 
resulted  Ir.  the  establishment  of  a  proper  criterion  under  which  a  scattering 
object  could  be  considered  acoustically  rigid;  namely,  that  the  frequency  of 
the  incident  plane  wave  be  well  below  that  of  the  lowest  mechanical  resonance 
of  the  scatterer.  The  non-rigid  spherical  case  was  studied  experimentally  by 
Hampton  and  McKinney  [4]. 

All  of  the  studies  mentioned  above  assumed  that  the  linear  wave  equation 
was  sufficient  to  describe  the  situation  of  interest.  However,  the  exact 
equations  describing  acoustic  wave  propagation  in  a  fluid  (as  well  as  for 
solids  and  plasmas)  are  actually  nonlinear.  The  linear  wave  equation  is  only 
an  approximation  that  is  valid  for  small  amplitude  behavior.  Deviations  In 
behavior  from  that  predicted  by  the  linear  wave  equation  can  become  signifi¬ 
cant  when  the  Mach  number  (ratio  of  particle  velocity  to  phase  velocity  or, 
equivalently,  the  ratio  of  the  change  in  mass  density  to  the  density  of  the 
undisturbed  fluid)  is  not  much  less  than  unity  [5J.  In  this  case,  the  non¬ 
linear  wave  equation  is  required  to  accurately  represent  the  behavior.  The 
study  of  acoustical  behavior  requiring  use  of  the  nonlinear  wave  equation  is 
called  nonlinear  or  finite-amplitude  acoustics. 
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Most  solutions  of  the  nonlinear  wave  equation  for  fluids  have  been  re¬ 
stricted  o  plane  waves  In  homogeneous  Infinite  media.  In  1860,  Earnshaw  [6] 
obtained  an  implicit  solution  to  the  lossless,  one-dimensional,  nonlinear  wave 
equation  subject  to  a  boundary  condition  at  the  origin.  That  solution  is 
valid  at  propagation  distances  small  relative  to  the  plane-wave  discontinuity 
distance  (that  propagation  distance  at  which  the  solution  to  the  lossless 
nonlinear  wave  equation  for  a  wave.^stnusolda  I  .it  Its  origin  becomes  multiple¬ 
valued.  This  occurs  because  points  of ‘high  particle  velocity  in  the  wave  also 
have  higher  propagation  velocities  and  hence  tend  to  overtake  the  points  of 
low  particle  velocity,  causing  the  waveform  to  approach  a  sawtooth  shape.  It 
Is  proportional  to  the  Mach  number  and  is  a  convenient  measure  of  the  non¬ 
linearity  of  the  problem.*)  Earnshaw's  solution  illustrated  that  points  of 
high  particle  velocity/pressure  in  the  time  waveform  (i.e.,  the  variation  with 
time  of  the  particle  velocity/pressure  at  a  fixed  position)  move  more  rapidly 
than  points  of  low  particle  velocity/pressure.  This  causes  the  time  waveform 
of  a  finite  amplitude  wave  to  change  its  shape  (i.e.,  distort)  as  the  wave 
propagates.  Investigators  of  this  problem  usually  assume  an  acoustic  wave 
that  Is  sinusoidal  (i.e.,  harmonic)  at  Its  point  of  origin  and  utilize  a 
harmonic  time  analysis  of  the  waveform  to  describe  the  subsequent  propagation 
of  the  wave.  The  distortion  of  the  waveform  from  its  initial  state  manifests 
Itself  in  the  generation  of  harmonic  component  waves.  In  effect,  a  wave  of 
angular  frequency  w  generates  waves  of  angular  frequencies  2u>,  3 u,  etc.  as  it 
propagates.  These  harmonic  components  usually  gain  energy  at  the  expense  of 
the  fundamental  component  of  angular  frequency  w. 


*In  a  real  situation,  the  solution  never  actually  becomes  multiple  valued.  It 
is  prevented  from  doing  so  due  to  energy  loss  due  to  the  viscous  terms  that 
are  no  longer  negligible  when  the  discontinuity  distance  is  approached. 
(Nonetheless,  the  onset  of  shock-wave  formulation  occurs  near  the  disconti¬ 
nuity  distance.)  In  higher-dimensioned  problems  there  is  no  corresponding 
discontinuity  distance  since  geometric  spreading  is  sufficient  to  prevent  this 
catastrophic  growth  of  the  nonllnearly  generated  wave.  Nonetheless,  it  repre¬ 
sents  a  conservative  estimate  of  the  distance  to  which  the  lossless  theory  may 
be  applied.  It  is  given  by  Beyer  (Reference  5,  p.  104)  as: 

(l/l)  *  [LKB/2A)|[(um0)/c02l, 

where  A  *  P0(3P/3p)s  >p,po=  P0c02 ,  B  =»  PQ2(32P/3p2)s  ,  l  => 

discontinuity  distance.  For  water  at  2 0°C,  B/A  is  approximately  5.0. 
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la  1935,  Fublni-Ghiron  [7]  obtaiaed  an  explicit  form  for  Earnshaw's 
implicit  solution  for  the  case  of  an  initially  harmonic  wave.  He  expressed 
his  solution  as  a  Fourier  harmonic  series  with  Bessel  function  coefficients. 

In  I960,  Kech  and  Beyer  [8]  obtained  a  solution  to  the  problem  of  the 
propagation  of  an  initially  harmonic  plane  wave  including  linear  absorptive 
losses.  To  obtain  this  solution,  they  assumed  that  each  of  the  acoustic 
variables  can  be  written  in  a  perturbation  series  in  which  the  order  of 
magnitude  of  each  term  is  smaller  than  that  of  the  preceding  term  by  a  factor 
equal  to  the  Mach  number.  They  then  calculated  in  succession  the  first  six 
terms  of  the  perturbation  series  solution  for  the  particle  displacement.  This 
gives  an  approximation  to  the  exact  solution  that  is  useful  for  distances  of 
propagation  somewhat  less  than  the  plane-wave  discontinuity  distance. 

In  1958,  Hayes  [9]  succeeded  in  putting  the  equations  of  motion  (includ¬ 
ing  thermoviscous  losses)  into  the  form  of  a  Burgers’  equation  [10].  Hayes 
does  not  make  clear  in  this  derivation  to  what  order  in  Mach  number  this  equa¬ 
tion  is  valid.  In  1963,  however,  Blackstock  [11]  again  put  the  equations  of 
motion,  including  thermoviscous  losses,  into  the  form  of  a  Burgers'  equation. 
This  treatment  makes  it  quite  clear  that  the  equation  is  valid  to  second  order 
in  Mach  number.  In  Reference  11,  Blackstock  introduces  the  term  "substitution 
corollary"  for  the  standard  procedure  used  to  identify  the  ordering  of  terms. 
In  essence,  this  corollary  states  that  in  obtaining  a  second-order  approxima¬ 
tion,  the  individual  acoustic  variables  involved  in  forming  any  term  comprised 
of  a  product  of  acoustic  variables  may  freely  be  replaced  by  their  first-order 
equivalents.  A  more  precise  substitution  would  generate  terms  higher  than 
second  order.  This  formulation  has  great  utility  in  that  exact  solutions  of 
Burgers'  equation  exist.  Blackstock  [11,12]  succeeded  in  obtaining  such  a 
solution  for  an  initially  harmonic  plane  wave. 

In  1964,  Blackstock  [13]  showed  that  when  the  spatial  coordinate  is  large 
relative  to  a  wavelength,  Burgers'  equation  can  also  be  used  to  solve  the 
problem  of  propagation  of  spherical  and  cylindrical  waves  in  a  lossless  medium 
(if  the  frequency  is  not  large,  losses  may  still  be  ignored  large  propagation 
distances).  In  1981,  Trivett  and  Van  Buren  [14]  developed  a  numerical  method 
of  solving  the  Burgers'  equation  for  plane,  cylindrical,  and  spherical  waves 
including  losses. 
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Significant  nonlinear  generation  of  acoustic  waves  can  also  occur  when 
two  waves  of  different  frequency  are  present  simultaneously  In  a  fluid  me¬ 
dium.  In  1948,  Eckart  (15)  derived  a  second-order  nonlinear  wave  equation 
(re-derived  by  Westervelt  (16)  In  1957)  useful  In  obtaining  solutions  to  this 
type  of  problem.  In  1963,  Westervelt  {17}  considered  the  problem  of  two  col- 
linea,  plane  waves  (called  primaries)  with  different  Initial  harmonic  time  de¬ 
pendences  at  the  origin.  He  started  with  Llghthiil's  [18]  equations  of  motion 
and  retained  terms  up  to  the  quadratic  In  Mach  number.  In  addition  to  the 
harmonics  of  the  primaries  predicted  when  a  plane  wave  propagates  and  distorts 
In  a  fluid  medium,  waves  at  frequencies  equal  to  the  sum  and  difference  of  the 
primary-wave  frequencies  were  also  predicted.  These  sum-  and  difference- 
frequency  waves  tend  to  grow  with  increasing  distance  from  the  origin.  In 
this  paper,  Westervelt  essentially  transformed  Eckart's  second-order  nonlinear 
wave  equation  into  a  form  known  as  the  simple-source  formulation.  This 
transformation,  however,  was  restricted  to  the  case  of  plane  wave  primaries. 

Westervelt's  theory  was  confirmed  experimentally  by  Beilin  and  Beyer 
[19],  who  produced  a  1-MHz  difference-frequency  wave  by  driving  a  2 . 54-cm~d iam 
circular  piston  source  at  the  primary  frequencies  13  and  14  MHz.  Agreement 
between  theory  and  experiment  was  good. 

In  1962,  Dean  [20]  presented  a  solution  for  the  sum-frequency  wave  pro¬ 
duced  by  two  outgoing,  concentric  cylindrical  waves.  He  again  started  from 
the  basic  conservation  equations  and  derived  a  coupled  set  of  differential 
equations  (accurate  to  second  order  in  Mach  number)  in  terras  of  a  new  set  of 
variables  (defined  in  terms  of  operations  on  the  first-  and  second-order 
acoustic  quantities).  Dean  stated  that  the  solution  he  presented  to  these 
equations  for  this  particular  case  was  exact.  Lauvstad  [21]  later  stated  that 
Dean's  solution  to  tills  problem  was  Incorrect  except  in  the  farfleld.  He 
stated  that  direct  substitution  of  Dean's  solution  into  the  equations  demon¬ 
strated  their  incorrectness  except,  as  previously  mentioned,  in  the  asymptotic 
limit  approaching  the  farfleld.  This,  however,  is  not  so.  Dean  later  per¬ 
formed  measurements  that  were  in  fair  agreement  with  the  qualitative  aspects 
of  his  theory  [22].  (No  attempt  was  made  to  demonstrate  quantitative 
agreement  with  the  theory.)  These  measurements  were  performed  under  farfleld 
conditions.  Lauvstad  offered  his  own  general  expression  for  the  sum-frequency 
component;  however,  he  used  the  Green's  function  corresponding  to  the 
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unphysical  case  of  a  soft  boundary.  A  correct  expression  for  the  difference- 
frequency  pressure  for  this  case  will  be  obtained  in  this  report. 

In  1972,  Censor  [23]  solved  the  linear  wave  equation  for  the  problem  of 
the  scattering  of  an  acoustic  plane  wave  (of  angular  frequency  <o)  by  an  ob¬ 
stacle  whose  surface  deforms  harmonically  (at  angular  frequency  £1).  In  his 
calculation.  Censor  included  nonlinearities  only  in  the  boundary  conditions  he 
imposed.  He  predicted  that  in  addition  to  the  usual  rigid-body  scattered 
field,  waves  at  angular  frequency  =  w  ±  £1  would  be  created  at  the 
boundary.  These  waves  would  then  propagate  outward  from  the  boundary  with  a 
behavior  described  by  the  linear  wave  equation. 

Shortly  after  the  appearance  of  Censor's  article,  Rogers  [24]  pointed  out 
that  waves  at  angular  frequencies  arising  from  medium  nonlinearities  would 
also  be  predicted  by  the  nonlinear  wave  equation.  Rogers  also  stated  that  the 
effects  predicted  by  Censor  depend  on  the  Mach  number  in  the  same  way  as  the 
effects  arising  from  nonlinear  theory.  Hence,  the  propagation  problem  cannot 
be  linearized  in  any  physically  meaningful  way  without  also  eliminating  the 
boundary  effect  predicted  by  Censor.  In  other  words,  a  solution  to  the 
problem  of  the  generation  of  sum-  and  difference-frequency  waves  must 
necessarily  involve  solving  the  nonlinear  wave  equation. 

In  the  present  study,  the  simple-source  formulation  of  the  second-order, 
nonlinear  wave  equation  for  a  lossless  medium  is  derived  for  arbitrary  primary 
wave  fields  of  harmonic  time  dependence.  This  equation  was  previously  derived 
by  Westervelt  [25];  however,  his  treatment,  as  previously  mentioned,  was  re¬ 
stricted  to  plane-wave  primaries.  The  assumption  that  no  linear  losses,  such 
as  those  due  to  viscosity  or  heat  conduction,  exist  in  the  fluid  medium  places 
restrictions  on  the  subsequent  solution.  These  restrictions  tend  to  increase 
with  increasing  viscosity  and  frequency.  However,  the  restrictions  are  not 
expected  to  be  significant  for  the  case  of  a  water  medium  and  the  frequencies 
to  which  the  solution  will  be  normally  applied  (less  than  about  200  kHz).  In 
obtaining  this  equation,  all  terms  up  to  the  quadratic  in  Mach  number  are  re¬ 
tained  in  the  acoustic  variables.  A  perturbation  approach  is  not  used  until 
the  final  step  in  the  derivation;  hence,  the  equation  upon  which  the  simple- 
source  formulation  is  based  remains  valid  even  when  the  second-order  effects 
cause  a  significant  energy  drain  on  the  primaries  (which  invalidates  a 
perturbat ion-series  approach). 


; 


In  obtaining  solutions  to  this  equation,  a  Born-approximation  type  of 
perturbation  analysis  is  used.  Here,  first-order  acoustic  wave  fields  are 
calculated  as  solutions  to  the  linear  wave  equation  and  then  used  to  determine 
the  inhomogeneous  term  of  the  second-order  equation.  The  perturbation  analy¬ 
sis  is  in  terms  of  the  Mach  number.  (Perturbation  expansions  in  terms  of  Mach 
number  to  solve  acoustic  problems  have  been  used  extensively  before 
(15,16,26,27].) 

The  problem  of  the  generation  of  sum-  and  dif ference-f requency  waves  via 
the  nonlinear  scattering  of  acoustic  waves  by  vibrating  obstacles  is  then 
addressed  for  three  geometries: 

1.  Plane  wave  normally  incident  on  a  uniformly  vibrating  infinite  plane. 

2.  Cylindrical  wave  incident  on  an  infinitely  long  cylinder  vibrating 
uniformly  in  the  radial  direction.  (The  symmetry  axes  of  the  incident 
wave  and  the  scattering  cylinder  are  assumed  parallel  but  not 
coincident . ) 

3.  Plane  wave  normally  incident  on  an  infinitely  long  cylinder  vibrating 
uniformly  in  the  radial  direction. 

The  first  case  above  is  readily  solved  aEter  expressing  the  second-order 
nonlinear  wave  equation  in  one-dimensional  form,  due  to  the  resulting  simplic¬ 
ity  of  the  calculations.  Solution  of  the  last  two  cases  is  much  more  diffi¬ 
cult.  The  approach  taken  is  to  formulate  the  solutions  in  terms  of  a  Green's 
function.  Care  must  be  taken  to  choose  the  proper  Green's  function  for  evalu¬ 
ating  the  appropriate  Born  Integral;  l.e.,  the  one  corresponding  to  the 
boundary  surfaces  Involved  (this  requirement  has  been  discussed  previously  in 
a  paper  presented  by  the  author  [28]).  The  resulting  expressions  for  the 
acoustic  pressure  of  the  sum-  and  difference-frequency  components  involve  some 
rather  complicated  integrals.  A  new  integration  procedure  was  developed  that 
allows  the  evaluation  of  these  integrals  in  closed  form  for  the  case  of  two 
high-frequency  primaries.  This  procedure  is  described  in  the  Appendix.  Nu¬ 
merical  results  are  obtained  by  Gaussian  quadrature  integration  for  the  more 
general  case. 
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The  results  obtained  by  the  Censor-method  approach  are  also  calculated 
and  presented  for  the  two  cases  for  which  his  theory  are  applicable — those 
involving  incident  plane  waves.  It  is  demonstrated  that  Censor's  prediction 
is  of  the  same  order  as  pseudosound*. 

For  the  case  of  a  plane  wave  incident  on  a  vibrating  cylinder,  the  re¬ 
sults  are  presented  graphically  for  both  the  Censor  theory  and  the  nonlinear 
theory. 

An  experimental  investigation  was  undertaken  to  confirm  the  theoretical 
prediction  for  the  difference-frequency  pressure  produced  in  the  case  in  which 
a  plane  wave  is  incident  on  a  vibrating  cylindrical  surface.  Although  the 
experiment  was  unsuccessful  in  confirming  the  theoretical  predictions,  it  was 
nonetheless  successful  in  identifying  several  of  the  difficulties  that  arise 
In  nearfield,  nonlinear  difference-frequency  experiments  and  solving  all  but 
one  of  those  identified.  The  most  significant  of  these  difficulties  (not 
encountered  in  previously  published  work)  are: 

•Inadvertent  direct  radiation  of  the  sources  at  the  difference  fre¬ 
quency:  This  will  tend  to  be  a  greater  source  of  error  when  the  measur¬ 
ing  hydrophone  is  near  the  sources  (since  difference-frequency  pressures 
produced  by  medium  nonlinearittes  tend  to  grow  with  distance  while 
directly  radiated  different-frequency  pressure  tends  to  decay  with 
distance). 

•Electrical  filtering  problems  due  to  experimental  constraints:  The 
difference  frequency  was  only  about  one  half  the  lowest  primary 


*Pseudosound  is  an  effect  arising  from  the  uncertainty  in  the  motion  of  the 
measuring  hydrophone;  i.e.,  the  uncertainty  in  the  motion  of  the  hydrophone  is 
of  the  same  order  as  the  difference  between  Langrangian  and  Euler ian  coord¬ 
inates.  Although  pseudosound  is  of  second-order  in  Mach  number  (as  are  the 
source  terms  of  the  second-order  nonlinear  wave  equation),  it  nonetheless  re¬ 
mains  a  minor  component  of  the  second-order  field.  This  is  because  contribu¬ 
tions  to  the  pressure  predicted  by  the  acoustic  second-order  nonlinear  wave 
equation  are  cumulative  with  respect  to  propagation  distance  and  hence  tend  to 
overwhelm  pseudosound  (which  is  a  function  only  of  the  magnitude  of  the  acous¬ 
tic  variables  at  the  observation  point)  within  propagation  distances  that  are 
a  fraction  of  a  wavelength.  Example  calculations  of  pseudosound  relevant  to 
the  current  research  are  provided  in  Sections  II.  E.  3  and  IV.  F. 
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frequency.  In  addition,  the  pulse  lengths  had  to  be  less  than  about  10 
cycles  at  the  difference  frequency  to  avoid  Interfering  reflections  from 
neighboring  surfaces.  Hence  the  usual  passive  methods  employed  for 
electrical  filtering  in  previous  farfleld,  nonlinear  measurements  were 
inappropriate . 

•  Difference-frequency  voltage  generated  nonllnearly  In  the  hydrophone: 
This  effect,  due  to  nonlinear  mixing  of  the  primaries  in  the  hydrophone, 
provided  larger  difference-frequency  voltages  than  those  produced  by  the 
difference-frequency  pressure  generated  by  nonlinearities  of  the  fluid 
medium.  The  effect  was  observed  for  a  wide  range  of  available 
hydrophones. 

Solutions  found  to  the  first  two  of  the  above  difficulties  will  be 
discussed  in  the  section  concerning  experimental  results.  Although  the  third 
difficulty  has  not  been  resolved,  several  valuable  observations  were  made  and 
are  presented  in  the  section  on  the  determination  of  hydrophone  nonlinearity 
(Section  IV.  F). 

In  light  of  the  insights  gained  by  these  hydrophone  nonlinearity  measure¬ 
ments,  the  program  at  the  Underwater  Sound  Reference  Detachment  (USRD)  of  the 
Naval  Research  Laboratory  has  been  significantly  expanded.  The  technique  to 
measure  hydrophone  nonlinearity  developed  during  this  research  represents  a 
substantial  improvement  over  the  technique  originally  intended  to  be  used  by 
USRD  to  analyze  the  linearity  of  their  standard  hydrophones.  In  the  current 
development  program,  it  is  intended  to  use  the  nonlinearity  measurement  tech¬ 
nique  of  this  research  to  develop  a  hydrophone  nonlinearity  standard. 

Chapter  II  will  present  the  nonlinear  theory  as  well  as  Censor's 
theory.  The  numerical  results  of  each  theory  will  be  given  in  graphical  form 
In  Chapter  III,  along  with  a  discussion  of  the  numerical  techniques  used. 
Chapter  IV  will  discuss  the  experiment,  including  the  choice  of  experiment  and 
resolution  of  experimental  difficulties,  and  will  present  a  new  technique  for 
determining  hydrophone  nonlinearity.  Lastly,  Chapter  V  will  give  the 
conclusions.  The  new  technique  of  integration  developed  during  the  course  of 
this  research  is  presented  in  the  Appendix. 


1 1 .  THEORY 


A.  Introduction 

This  chapter  presents  some  concepts  useful  in  finite-amplitude  acoustics, 
gives  the  derivation  of  the  simple-source  formulation  of  the  second-order 
acoustic  wave  equation  in  a  manner  that  is  of  quite  general  applicability  to 
problems  involving  interacting  acoustic  pressure  fields,  and  presents  the 
theoretical  development  of  the  problem  of  the  nonlinear  scattering  of  acoustic 
waves  from  vLbratlng  obstacles  for  certain  specific  geometries* 

Section  II*  B  discusses  the  two  reference  frames  used  in  acoustics — 
namely,  the  Lagrangian  (or  material)  coordinates  and  the  Eulerian  (or  spatial) 
coordinates*  Section  II.  C  discusses  the  orders  of  acoustic  variables  and 
expressions. 

In  1963,  Westervelt  [17]  obtained  a  simple-source  formulation  of  the 
second-order  nonlinear  wave  equation.  This  formulation  is  an  analogy  to  the 
simple-source  wave  equation  of  linear  acoustics,  which  is  essentially  the 
inhomogeneous  wave  equation  for  a  volumetric  distribution  of  monopole  point 
sources  of  sound  [29]. 

The  analogy  drawn  by  Westervelt  is  that  each  elementary  volume  element  In 
a  fluid  in  which  two  waves  of  different  frequency  are  simultaneously  present 
may  be  viewed  as  an  elementary  source  of  nonllnearly  generated  waves.  The 
mathematical  form  of  Westervelt's  second-order  nonlinear  wave  equation  is 
similar  to  that  of  the  simple-source  equation  of  linear  acoustics  if  the 
proper  identification  of  variables  is  made. 

Westervelt's  derivation  of  this  equation  is  based  on  the  assumption  that 
the  interacting  waves  are  planar.  Hence,  it  is  unclear  that  his  equation  is 
applicable  to  any  other  wave  geometry. 

The  inhomogeneous  term  is  quadratic  in  nature.  Hence,  if  one  attempts  to 
obtain  the  solution  for  arbitrary  wave  fields  by  decomposing  the  individual 
waves  contributing  to  the  source  term  into  sums  of  plane  waves,  the  solution 
must  be  represented  as  a  sum  over  pairs  of  waves.  It  is  often  more  convenient 
to  use  a  closed-form  representation  of  the  source  term  (but  this  requires  a 
demonstration  that  the  equation  is  valid  for  non-planar  primaries).  In  the 
current  work,  the  simple-source  formulation  of  the  second-order  nonlinear  wave 

9 


i 


equation  is  re-derlved  for  primary  waves  of  arbitrary  geometry.  This 
derivation  Is  given  in  Section  II.  D. 

Sections  II.  E,  F,  and  G  present  the  solutions  to  the  second-order 
nonlinear  wave  equation  for  three  specific  geometries: 

•Plane-wave  scattering  from  a  vibrating  planar  surface. 

•Cylindrical-wave  scattering  from  a  vibrating  cylindrical  surface. 

•Plane-wave  scattering  from  a  vibrating  cylindrical  surface. 

In  solving  the  case  of  plane-wave  scattering  from  an  infinite  uniformly 
vibrating  planar  surface,  the  one-dimensional  form  of  the  nonlinear  wave  equa¬ 
tion  is  used.  This  equation  is  expressed  in  terms  of  Lagrangian  coordinates 
since  the  first-order  boundary  conditions  can  most  naturally  be  satisfied  in 
this  reference  frame.  The  equation  is  solved  by  a  substitution  of  perturba¬ 
tion-series  expansions  in  Mach  number  for  the  acoustic  variables  such  that  the 
solution  is  accurate  to  second  order  (in  the  sense  normally  associated  with 
perturbation  approximations).  As  noted  by  Beyer  [30],  such  an  expansion  was 
first  considered  historically  by  Airy  in  1845  in  studying  tidal  motion.  A 
careful  analysis  of  the  results  of  Censor's  theory  is  presented  for  the  planar 
case.  (Although  Censor  presented  no  solution  for  this  case,  his  method  can  be 
applied  in  a  straightforward  manner  to  obtain  one.)  It  Is  shown  that  Censor's 
theory  predicts  sum-  and  difference-frequency  pressures  that  are  of  the  same 
order  as  pseudosound. 

In  obtaining  a  solution  to  the  second-order  nonlinear  wave  equation  for 
the  two  cases  considered  involving  cylindrical  geometry,  again  the  method  of 
expanding  the  acoustic  variables  in  a  perturbation  series  in  Mach  number  Is 
used. 

The  actual  second-order  sum-  and  difference-frequency  pressures  are 
calculated  by  solving  a  related  Green's  function  equation.  It  has  been  noted 
[28]  that  care  must  be  exercised  in  the  choice  of  Green's  function.  This  is 
d Iscussed  in  detail  in  the  paragraphs  following  Eq.  (47)  of  Section  II.  D. 

Finally,  both  Censor's  theoretical  results  and  the  results  of  the  non¬ 
linear  theory  are  presented  graphically  to  facilitate  comparison  of  the  two 
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theories.  (Although  Censor  gave  analytical  expressions  for  the  scattered  sum- 
and  difference-frequency  pressures,  no  numerical  values  were  presented.) 

tn  obtaining  the  numerical  results  for  the  nonlinear  theory,  the  inte¬ 
grals  were  performed  on  a  high-speed  computer  (TI-ASC-11)  using  the  method  of 
Gauss  Quadrature  [32].  (Closed-form  solutions  to  these  Integrals  were  obtain¬ 
ed  in  the  high-frequency  limit  based  on  the  new  integration  technique  pre¬ 
sented  in  the  Appendix.)  The  sums  obtained  from  the  theory  were  also  carried 
out  on  this  computer  until  three  significant  figures  were  obtained. 

B.  Coordinate  Systems  of  Finite-Amplitude  Acoustics 

Two  different  kinds  of  coordinate  systems  are  used  to  specify  acoustic 
wave  fields  in  fluids — namely,  Lagrangian  (or  material)  coordinates  that  move 
with  the  fluid  and  Eulerian  (or  spatial)  coordinates  that  are  fixed  in 
space.  The  relationship  between  these  two  systems  is  illustrated  in  Fig.  1. 

When  displacement  is  expressed  in  terms  of  Eulerian  coordinates,  the 
displacement  is  that  of  the  fluid  element  that  happens  to  be  located  at  x  at 
the  time  of  observation.  On  the  other  hand,  when  the  displacement  is 
expressed  In  terms  of  Lagrangian  coordinates  a,  it  refers  to  a  fluid  element 
that  had  the  initial  rest  position  a. 

I  <«  _  £ 


_ X  _ 

REFERENCE 

POSITION 

x  =  a  +  £ 

a  -  LAGRANGIAN  COORDINATE 
*  REST  POSITION  OF  AN  INDIVIDUAL  PARTICLE 

x  «  EULERIAN  COORDINATE 
=  POSITION  OF  A  FIXED  POINT  IN  SPACE 

f  *  DISPLACEMENT  OF  THE  PARTICLE  FROM  ITS 
REST  POSITION 

Fig.  1  -  Geometry  of  Lagrangian  and  Eulerian  Coordinates 
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C.  The  Order  of  Acoustic  Variables  and  Expressions 


Careful  consideration  must  be  given  to  the  meaning  of  "orders"  of  acous¬ 
tic  variables  and  expressions.  in  nonlinear  acoustics,  the  appropriate  param¬ 
eter  to  determine  this  order  Is  the  Mach  number.  (This  will  be  justified  pre¬ 
sently.)  Important  dimensionless  quantities  arise  when  the  relevant  equa¬ 
tion/s  are  put  into  dimensionless  form.  The  resulting  equation  was  given  for 
the  acoustical  case  by  Blackstock  [13].  A  relationship  (analogous  to  the 
Reynolds  number)  was  obtained  by  Blackstock  in  this  report,  which  may  be  used 
to  determine  whether  the  loss  terms  are  significant  for  a  nonlinear  propaga¬ 
tion  problem  for  Initially  plane  waves.  This  criterion  involves  a  quantity  P, 
where 


r  *  BecQXc/(l/2)y[V  +  (Y-1)/Prl 

(B  Is  a  measure  of  medium  nonlinearity  and  is  approximately  3.5  for  water  at 
20°C,  V  is  the  viscosity  number,  y  is  the  ratio  of  specific  heats,  v  is  the 
kinematic  viscosity,  Pr  is  the  Prandtl  number,  c  Is  the  Mach  number,  and  xc  i? 
a  characteristic  length).  alackstock  noted  that  for  initially  sinusoidal 
waves,  xc  may  be  taken  to  be  xc  *  c0/u>  =  A,  and  showed  that  (at  low  fre¬ 
quencies)  T  =  Be/aA,  where  a  =  the  attenuation  coefficient.  Blackstock  noted 
that  the  quantity  Be/aA  has  been  referred  to  as  the  Reynolds  number  in  the 
Russian  literature.  (He  further  noted  that  this  quantity  is  certainly  not  a 
measure  of  Inertial  to  viscous  effects,  which  is  the  traditional  interpreta¬ 
tion  of  the  Reynolds  number.  Instead,  Blackstock  interprets  this  quantity  as 
a  measure  of  the  importance  of  nonlinearity  to  dissipation).  The  determina¬ 
tion  of  whether  a  given  frequency  is  "low"  can  be  ascertained  using  yet  an¬ 
other  dimensionless  quantity  known  as  the  "frequency  parameter"  [33], 

X  =  (u>nV)/poco  ,  where  V  =  viscosity  number  =  2  +  n’/n,  n  =  shear  coefficient 
of  viscosity,  n'  =  volume  coefficient  of  viscosity.  Using  a  typical  frequency 
of  Interest  of  100  kHz  along  with  the  approximate  expression  n'  a  3n  for  water 
[34],  we  get  a  value  X  =  1.4x10“**,  clearly  indicating  that  this  frequency  may 
properly  be  considered  "low"  (and  hence  Blackstock' s  low-frequency  expression 
for  T  Is  appropriate). 

An  Important  criterion  established  by  Blackstock  In  Reference  13  is  that 
the  results  of  the  propagation  of  a  plane  wave  that  would  be  obtained  via  the 


inclusion  of  loss  In  the  equations  will  closely  approach  the  results  obtained 
excluding  loss  when  T>10.  Using  typical  values  of  interest  in  the  current 
research  (  3  *  3.5,  e  »  5*10-^ ,  a  "  10  ^ra  ^ ,  and  A  "  10  ^m) ,  we  obtain 
T  »  1100.  Although  Blackstock's  expressions  are  actually  valid  only  for 
initially  plane  waves,  it  is  fairly  clear,  due  to  this  rather  substantial 
value,  that  losses  are  relatively  unimportant  here. 

The  negligibility  of  viscous  effects  for  the  frequencies  and  propagation 
distances  of  interest  in  this  report  can  also  be  demonstrated  by  computation 
of  the  Reynold's  number  in  the  two  physical  regions  of  concern:  1)  the  region 
close  to  the  scatterer's  surface  (where  a  viscous  boundary  layer  forms),  and 
2)  the  propagation  region  wherein  (at  some  point)  the  primary  fields  will 
start  to  form  shock  waves. 

Consider  first  the  required  thickness  of  a  viscous  boundary  layer.  We 

define  this  thickness  as  corresponding  to  a  Reynold's  number  of  unity.  The 
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Reynold's  number  (R)  can  be  computed  using  R  =  ULp/p  or  R  =  wL  p/u,  where 
U  =  characteristic  velocity,  L  ■  characteristic  length,  p  =  coefficient  of 
viscosity,  and  m  =  characteristic  angular  frequency.  Using  a  frequency  of 
100  kHz,  a  Reynold's  number  of  unity  corresponds  to  a  characteristic  (boundary 
layer)  thickness  of  about  10~^  cm.  This  result  makes  clear  the  fact  that  at 
positions  near  the  scatterer's  surface  viscous  effects  may  reasonably  be  ne¬ 
glected  (since  this  distance  is  «A/2n  for  the  frequencies  of  interest  here). 

At  larger  distances  from  the  scatterer's  surface,  this  analysis  breaks 
down,  and  a  more  suitable  interpretation  of  the  Reynold's  number  is  needed. 

As  an  initially  sinusoidally  shaped  acoustic  wave  propagates,  the  nonlinear 
distortion  it  suffers  causes  its  waveform  to  approach  the  shape  of  a  sawtooth 
wave.  If  viscous  effects  are  completely  neglected,  the  waveform  becomes  a 
triangular  wave  at  a  propagation  distance  equal  to  the  discontinuity  distance 
(see  the  discussion  of  this  matter  in  Chapter  I).  Using  this  effect  as  a 
guide,  it  is  clear  that  a  convenient  choice  for  the  characteristic  length  L 
used  to  compute  the  Reynold's  number  is  the  spatial  distance  separating  the 
point  of  maximum  particle  velocity  and  adjacent  point  of  zero  particle 
velocity.  This  "adjacent"  point  of  zero  particle  velocity  is  the  one  occurlng 
ahead  of  the  point  of  maximum  particle  velocity  (in  the  sense  of  the  direction 
of  propagation).  This  characteristic  length,  as  defined  here,  is  initially 
equal  to  one  quarter  of  a  primary  wavelength  and  approaches  zero  as  the 
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propagation  distance  approaches  the  discontinuity  distance*  Such  a  definition 
makes  clear  the  fact  that  viscosity  can  properly  be  neglected  in  the  propaga¬ 
tion  region  only  at  propagation  distances  small  relative  to  the  discontinuity 
distance  (since  the  Reynold's  number  will  approach  zero  at  that  point). 

Hence,  we  conclude  that  within  a  region  farther  than  a  fraction  of  a 
millimeter  from  the  scatterer's  surface  and  yet  not  closely  approaching  the 
discontinuity  distance,  viscous  effects  may  safely  be  neglected  and  the  Mach 
number  taken  as  the  appropriate  dimensionless  parameter  to  identify  the  order 
of  acoustic  variables  and  expressions. 

Due  to  the  extremely  complex  nature  of  the  equations  of  nonlinear 
acoustics,  some  type  of  approximation  method  is  usually  required  to  obtain  a 
solution.  There  are  two  primary  approximations  that  are  traditionally  used: 

•The  exact  equations  can  be  put  into  an  approximate  form  more  readily 
solvable.  It  is  the  local  Mach  number  dependence  of  the  acoustic 
variables  (discussed  above)  that  forms  the  basis  for  this  approximation. 

•Alternatively,  the  acoustic  variables  can  be  expanded  directly  in  a 
fc rmal  perturbation  series  in  Mach  number,  which  may  then  be  substituted 
Into  the  exact  equations.  In  such  a  series,  "ordering”  Is  determined  by 
where  the  boundary  conditions  are  imposed  (i.e.,  the  “Mach  number"  is  not 
the  local  Mach  number  but  rather  the  Mach  number  at  the  point  at  which 
the  boundary  conditions  are  Imposed).  This  process  results  in  an 
infinite  set  of  differential  equations  (one  equation  associated  with  each 
term,  i.e.  order  of  the  expansion).  Each  of  these  equations  can  then  be 
individually  solved,  starting  with  the  first-order  equation  and  ending  at 
whatever  order  yields  the  required  degree  of  accuracy. 

Let  us  first  consider  the  approximation  involved  in  obtaining  the  second- 
order  nonlinear  wave  equation  from  the  exact  fundamental  equations.  Although 
this  equation  is  an  approximate  one,  when  It  Is  applied  to  problems  involving 
sound  propagation  In  water,  no  measurable  difference  between  the  pressure  pre¬ 
dicted  by  it  and  that  predicted  by  the  exact  equations  results.  This  is  be¬ 
cause  in  such  problems  the  contributions  to  the  sound  field  by  terms  of  third 
and  higher  orders  are  negligible.  Justification  for  the  negligibility  of 
third-order  terms  is  given  in  Appendix  A  of  Reference  35,  which  considers  the 
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consequences  of  retaining  higher-order  terms  in  the  fundamental  equations  and 
third-order  terms  in  the  equation  of  state.  In  Reference  36,  the  equation  of 
state  was  expressed  as 


* .  r\  .  •  /'V  .  c  /'v3 

p0  + » <— >  +  j  (-5—)  +  s  <-— >  • 

000 


As  reported  in  Reference  35,  Van  Buren  wrote  a  computer  program  that 
computed  the  distortion  occurring  during  the  propagation  of  an  initially 
sinusoidal  wave  of  amplitude  0.7xl05  Pa  and  frequency  2  MHz.  (The  program 
included  the  effects  of  all  orders  in  the  fundamental  equations  and  in  the 
equation  of  state  to  third  order.  Absorption  effects  were  also  included.) 

The  wave  was  allowed  to  propagate  104  cm  (one  discontinuity  distance).  The 
results  of  this  program  were  compared  to  results  obtained  when  only  terms  up 
to  second  order  were  retained  in  both  the  fundamental  equations  and  in  the 
equation  of  state.  The  results  of  comparing  these  two  solutions  were:  With 
C/A  ■  105 ,  the  second  harmonic  deviated  by  about  -0.2%,  the  third  harmonic  by 
-1.3%,  fourth  harmonic  by  +0.06%,  etc.  In  liquids,  C/A  is  approximately  (36] 
3/2  (B/A)^,  which  gives  C/A  *  40  for  water.  It  becomes  clear  from  these 
numerical  results  that  the  effects  of  the  higher-order  terms  (at  least  in 
water)  are  completely  negligible. 

Although  additional  effects  may  arise  in  non-planar  geometries,  the 
plane-wave  should  represent  the  worst  case  (since  higher-dimensioned  geome¬ 
tries  result  in  spreading  of  the  waves  and  a  reduction  in  field  amplitude). 

One  last  argument  can  be  advanced  regarding  the  negligibility  of  third- 
order  terms:  If  one  starts  with  the  exact  wave  equation  in  Lagrangian 
coordinates  for  plane  waves  and  performs  straightforward  Taylor-Series 
expansions,  one  can  demonstrate  that  the  third-order  source  terms  are  a  Mach 
number  smaller  than  the  second-order  source  terms.  Hence,  even  if  third-order 
source  terms  result  in  cumulative  effects  (as  do  the  second-order  terms),  the 
contributions  to  the  answer  from  the  third-order  terms  will  be  a  Mach  number 
lower  than  contributions  from  second-order  terms  (the  Mach  number  in  water  is 
rarely  greater  than  10"^).  This  holds  as  long  as  dispersion  is  negligible, 
which  it  is  in  fresh  water  for  frequencies  up  to  10^  Hz.  Again,  the  plane- 
wave  case  may  be  regarded  as  a  worst-case  situation,  since  geometrical 
spreading  will  dilute  the  effects  in  higher-order  geometries. 
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We  begin  this  discussion  by  referring  to  the  exact  wave  equation  In 
Lagrangian  coordinates  {see  Reference  37,  Eq.  (35)] 


« -  c2^ar\a 


An  expansion  for  the  wavespeed  c,  accurate  to  third  order,  has  been  given 
by  Van  Buren  [38].  It  is 

2  2 
. .  .  .  B  ,  , .  u  ,  .  c  B  .  B  ,  1 .  u  , 

C  =  Co11  +  (2A  +  1}  —  +  (5a  '  7T+  2A  +  2>  ~ 21’ 

o  4A  c 

o 

In  this  equation  A,  B,  and  C  have  their  usual  meanings — namely, 
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3P 

A  =  Po(~Tp)S,p=p  "  PoCo 
o 


B  °  Po2^)S,p=p  * 
dp  o 


J  „  3,aV 

and  C  =*  p  ( — y)  • 

°  V  S’p’po 

Substituting  this  expansion  for  the  wavespeed  into  the  exact  wave  equation  in 
Lagrangian  coordinates,  expanding  the  term  (l+£a)-^>  and  retaining  only  terms 
of  th’rd  order  or  less  yields: 


f  - 


c  2S  =  c  2[2(~  +  1)  —  -  2£  ]£ 
o  aa  o  2A  c  a  aa 
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The  quantities  (u/cQ),  £a ,  and  £aa  may  all  be  regarded  as  being  of  the 
order  of  the  Mach  number.  Hence,  the  above  form  of  the  wave  equation  can  be 
Interpreted  as  follows:  The  first  set  of  square  brackets  (witli  Its  coef¬ 
ficients)  may  be  regarded  as  a  source  of  second-order  waves.  The  second  set 
of  square  brackets  may  be  regarded  as  a  third-order  source  expression.  As 
previously  mentioned,  since  the  third-order  source  terms  are  all  a  Mach  number 
smaller  than  the  second-order  source  terms,  their  contributions  will  always  be 
a  Mach  number  smaller  than  the  contributions  from  the  second-order  terms. 
Hence,  even  if  these  third-order  contributions  are  cumulative  (as  are  the 
second-order  contributions),  they  will  always  remain  small  relative  to  the 
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second-order  contributions.  Hence,  tor  all  practical  purposes,  a  solution  to 
the  second-order  nonlinear  wave  equation  may  be  regarded  as  sufficient  for 
specifying  a  solution  to  a  nonlinear  underwater  sound  propagation  problem*. 

It  is  frequently  the  case,  however,  that  exact  solutions  to  even  this  more 
simplified  equation  are  too  difficult  to  obtain.  In  such  cases,  a  solution  is 
usually  obtained  by  using  a  perturbation-series  expansion  in  Mach  number  of 
the  acoustic  pressure.  This  expansion  is  generally  substituted  directly  into 
the  second-order  nonlinear  wave  equation  (rather  than  returning  to  the  still 
more  complicated  fundamental  equations).  Such  a  substitution  can  be  used  to 
obtain  what  is  known  as  the  simple-source  formulation  of  the  second-order 
nonlinear  wave  equation. 

When  one  uses  a  perturbation  expansion,  however,  one  must  proceed  with 
extreme  caution.  Although  the  starting  equations  may  be  regarded  as  suf¬ 
ficiently  accurate  (whether  starting  from  the  fundamental  equations  or  the 
second-order  nonlinear  wave  equation),  the  solution  obtained  by  a  perturbation 
series  may  become  inaccurate  if  the  series  is  truncated  too  early.  In  fact, 
most  such  perturbation  analyses  are  carried  only  to  second  order  (since  higher 
orders  become  exceedingly  complicated).  In  this  case,  inaccuracies  arise  due 
to  the  failure  of  the  fundamental  assumption  made  in  such  treatments,  namely 
that  the  second-order  field  remains  small  relative  to  the  first-order  field. 
Since  the  Mach  number  is  rarely  greater  than  10"^  in  water,  it  might  appear 
this  assumption  would  never  become  questionable  (in  fact,  it  may  appear  pecu¬ 
liar  that  second-order  effects  ever  become  measurable).  The  reason  they  do  so 
is  that  such  effects  tend  to  act  cumulatively  with  propagation  distance. 

Hence,  the  second-order  pressure  (for  example)  at  a  given  observation  point  is 
not  simply  a  consequence  of  the  value  of  the  Mach  number  at  that  point. 

Rather,  it  is  a  consequence  of  the  entire  integrated  history  of  the  fields 
between  the  sources  and  the  observation  point. 

In  essence,  this  means  that  the  important  second-order  contributions  to 

2 

the  pressure  are  not  themselves  second  order  (e  ),  but  rather  of  the  order 
2 

of  e  times  an  enhancement  factor.  What  this  factor  is  can  be  determined  for 
plane,  cylindrical,  and  spherical  geometries  due  to  the  fortuitous 


♦Hydrophones  can  rarely  be  calibrated  more  accurately  than  to  within  l  dB  of 
relative  error.  This  corresponds  to  more  than  12%  experimental  error. 
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circumstance  that  the  Bessel  function  expansion  used  by  Fubini-Ghi ron  in  the 
planar  case  ( 7  J  can  also  be  used  (when  therinoviscous  losses  are  negligible)  to 
solve  the  higher  dimensional  cases  (see  Reference  13,  eq .  (13)].  After 
Fubinl-Ch i ron ,  we  can  use  the  first  terms  of  this  series  to  estimate  the 
growth  behavior  of  the  nonlinear  fields.  in  order  to  do  so  conveniently,  we 
introduce  notation  similar  to  Blackstock' s  in  Reference  13: 

x 

O  *  — 

X 

a  (plane) 

f  =  2/a  (/a  -  /a~  (cylindrical) 

o  o 

Oq  log  (a/o^)  (spherical) 

where  again  e  *  Mach  number,  X  =  discontinuity  distance,  and  x  =  propagation 
distance.  (This  choice  of  Blackstock's  dimensionless  quantity  a  is  made  in 
order  to  render  his  Bessel-f unction  expansion  equivalent  to  that  of  Fubini- 
Ghi  ron.)  We  can  estimate  the  appropriate  "enhancement  factor"  via  the  ratio 
of  the  second-order  to  the  first-order  contributions  to  the  field  from  this 
Bessel -funct  ion  series.  This  ratio  is  equal  to  [J2(2f)]/[2  J  j_  ( f ) ) «  which  (for 
small  arguments  of  the  Bessel  function)  is  approximately  equal  to  V2  f. 

Hence,  the  enhancement  factors  for  each  of  the  three  geometries  becomes: 


\  £  (plane) 


V?  (Vs-Vf) 


(cylindrical ) 


1  x 

j  log  (£-)  (spherical) 
o 


where  x{|  represents  the  local  ion  of  the  scatterer's  surface  in  the  nonlinear 
scattering  problem  of  interest  here. 


Confirmation  of  the  growth  of  the  second-order  fields  relative  to  the 
primary  fields  Is  provided  by  the  many  successful  measurements  of  nonlinearly 
generated  field  effects  described  In  Chapter  I.  Therefore,  it  is  seen  that  In 
an  underwater  nonlinear  propagation  problem,  the  primary  fields  tend  to  decay 
(due  to  geometrical  spreading,  energy  loss  to  secondary  field  generation,  and 
energy  loss  due  to  absorption  when  high  frequencies  and/or  large  propagation 


distances  are  involved)  while  the  secondary  fields  tend  to  grow.  Eventually, 
the  secondary  fields  can  become  comparable  to  or  even  exceed  the  primary 
fields.  When  this  occurs,  solutions  based  on  second-order  perturbation 
methods  are  no  longer  valid. 

In  the  following  section,  the  second-order  nonlinear  wave  equation  is 
derived  from  the  fundamental  equations.  The  acoustic  pressure  is  then  ex¬ 
panded  in  a  perturbation  series  to  second  order.  This  perturbation  expansion 
Is  used  to  obtain  the  simple-source  formulation  of  the  second-order  non-linear 
wave  equation.  This  equation  together  with  the  first-order  equation  of  linear 
acoustics  constitutes  essentially  a  type  of  Born-approximation.  It  is  used  to 
solve  the  problem  of  the  nonlinear  scattering  of  acoustic  waves  for  vibrating 
obstacles  for  three  different  geometries  in  Chapter  II,  Sections  E,  F,  and  G. 
Since  these  solutions  are  obtained  via  a  perturbation  method,  their  validity 
is  restricted  to  small  propagation  distances  from  the  scatterer's  surface. 

D.  Second-Order  Nonlinear  Wave  Equation 

Any  investigation  of  the  behavior  of  finite-amplitude  acoustic  waves  in  a 
fluid  begins  with  the  basic  equations  of  motion.  These  can  readily  be  derived 
by  apply li. g  mass  and  momentum  conservation  laws  to  the  fluid  under  considera- 
t  ion.  The  resulting  equations  are  expressed  below  in  Eulerlan  coordinates: 

l.  The  Equation  of  Continuity  (Mass  Conservation) 

-§£+  $  •  (p£)  =  s.  (1) 

In  this  equation,  S  is  a  mass  source  term  representing  the  rate  at 
which  mass  is  introduced  into  the  region  of  interest.  It  can  be  used  to 
represent  raonopole  sources  of  sound  as  well  [39].  In  the  current  work, 
however,  sources  of  sound  will  be  handled  Instead  via  specification  of 
appropriate  boundary  conditions  (l.e.,  by  specifying  the  normal  velocity 
of  the  surface  of  the  source),  and  the  solutions  will  be  restricted  to 
regions  outside  the  surface  of  the  sources.  Hence,  the  term  S  will  be 
taken  to  be  zero. 
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where  is  as  defined  above. 

In  the  present  work,  frequencies  on  the  order  of  100  kHz  and  propagation 
distances  on  the  order  of  100  cm  In  fresh  water  will  be  considered.  The  ef¬ 
fects  of  viscous  attenuation  (for  a  plane  wave)  may  be  summarized  by  the  equa¬ 
tion  P  =»  P0e-ax.  At  a  frequency  of  100  kHz  in  Eresh  water,  the  constant  a  has 

a  value  of  ~10  ^  m  It  is  clear,  then,  that  viscous  losses  for  the  present 

work  are  completely  negligible  and  P^j  may  be  replaced  by  PS^j.  (In  fact,  the 
equations  for  sound  propagation  in  fresh  water  need  not  Include  the  effects  of 
viscosity  until  frequencies  approaching  1  MHz  or  distances  on  the  order  of  ki¬ 
lometers  are  considered.  Hence,  the  equations  herein  derived  have  a  broad  ap¬ 
plicability;  however,  viscous  effects  will  be  important  In  the  boundary  layer. 
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Substituting  Into  Eq.  (3)  yields: 


3  (Pu.u  ) 

- +  v  P 


itj  3x>3,1j 


which  may  be  used  to  show  that* 

.2 


□  2P 


at 


(Ap  -  c  AP)  -  l 

°  t,j 


3  (Puiu^> 
3x^  3xj 


(6) 


(7) 
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The  usual  wave  equation  of  linear  acoustics,  □  P  =  0,  follows  if  we  only 
retain  terms  that  are  linear  in  the  field  variables.  On  the  other  hand,  we 
need  to  retain  terms  up  to  the  quadratic  in  the  field  variables  in  order  to 
obtain  a  nonlinear  wave  equation  that  is  accurate  to  second  order.  As  stated 
earlier,  third  and  higher  order  terms  in  the  wave  equation  do  not  measurably 
contribute  to  nonlinear  acoustic  behavior  in  liquids  (such  as  water)  where 
these  results  will  be  applied.  Thus  the  term  pu^uj  can  be  replaced  by  its 
second-order  approximation  PgU^u j ,  thereby  neglecting  the  third-order  term  (p 
-  pQ)u^Uj.  The  last  term  in  Eq.  (7)  can  now  be  written: 


l 

i » j 


3  (PUjUj) 


{(u*$)($»u)  +  ($*U)2 
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i  J 


3u  .  3u 
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Now,  it  can  easily  be  shown  that 


3U  ,  3U  ,  y  .  y  y  .  y  y 

l  =  $  •  [(u*$)u]  -  (uM,(*.5>. 

i  *  j  ^  j 


Also,  by  vector  identity: 


(u*$)u  =  i  $(u2)  -  u  x  ($xu). 


(8) 


(9) 


(10) 


*In  obtaining  Eq.  (7)  a  factor  of  p  has  been  replaced  by  Ap  and  P  by  AP.  This 
may  be  freely  done  due  to  the  presence  of  the  differential  operators  acting  on 
these  quantities.  This  freedom  will  frequently  be  used  in  several  of  the 
equations  obtained  throughout  the  remainder  of  the  current  section. 
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It  can  bo  shown  that  In  linear  acoustics  the  particle  velocity  Held  is 
ir rotat Iona  1 ;  l.e.,  =  0  [tlie  validity  of  this  assumption  is  discussed  in 

Hunt  (37)  in  reference  to  his  Eq.  (56)].  Although  the  particle  velocity  is 
irrotational  only  to  first  order,  Rlackstock  [40]  has  pointed  out  that  any 
factor  in  a  second-order  term  may  be  replaced  by  its  first-order  equivalent 
since  a  more  precise  substitution  will  result  in  terms  of  third  (or  higher) 
order.  (Blackstock  calls  this  fact  the  substitution  corrollary.)  Since  the 
last  term  in  Eq.  (10)  is  clearly  of  second  order,  we  may  freely  use  the 
l rrotat ional ity  condition  in  this  term. 

Combining  the  equation  that  results  from  Eq .  (10)  by  using  this 
substitution  with  Eqs.  (8)  and  (9),  we  obtain: 
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i,j 


3  (pu^) 

3x^3xj 


Po[U-$)($*«i)  +  ($*u)2  +  V2(|  £2)]. 
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In  the  first  term  on  the  right-hand  side  of  Eq .  (11),  we  may  freely 
replace  $*u  by  equations  accurate  to  first  order,  since  the  overall  term  will 
remain  accurate  to  the  second  order  due  to  the  presence  of  u  dotted  into  the 
remainder  of  the  term.  Of  use  here  is  the  first-order  continuity  equation  as 
obtained  from  Eq.  (1): 


1 _ 1  3p 

P  3t 


(12) 


(The  "1”  over  the  equal  sign  denotes  first-order.  This  notation  shall  be 
adopted  for  the  remainder  of  this  thesis.  Similarly,  a  ”2"  over  the  equal 
sign  will  denote  second  order.) 

Therefore  the  first  term  on  the  right-hand  side  of  Eq .  (11)  becomes 


<ti't)($*u)  i  -  ~  ($p). 


(13) 


To  complete  the  analysis  we  need  an  equation  of 
medium.  Since  any  thermodynamic  quantity  in  systems 
volume,  and  temperature  are  thermodynamic  parameters 


state  for  the  fluid 
in  which  pressure , 
can  be  represented  as  a 
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function  of  any  two  state  variables,  we  can  obtain  an  equation  of  state  by 
expanding  the  pressure  in  a  Taylor  series  in  the  state  variables  p  (density) 
and  S  (entropy).  This  yields 


P  *  Po  + 


2 

(P-PQ)  +  \  (P  -P0)‘ 

S,P-pQ  3pZ  S,P=Pq 


+  •••  +  (-gf)  (S-S  )  +  ...  .  (14) 

s,p-Po 

It  is  usual  to  simplify  Eq.  (14)  under  the  assumption  of  adiabatic 
compressibility.  According  to  Morse  &  Ingard  [41],  adiabatic  compressibility 
is  achieved  under  the  condition  that  the  highest  frequency  component  in  the 
acoustic  field  is  significantly  less  than  (cQ^P0cp)/k,  where  k  is  the  thermal 
conductivity  and  cp  is  the  specific  heat  at  constant  pressure.  For  water, 
this  expression  gives  a  frequency  of  about  10^  Hz,  a  value  well  above 
anything  of  interest  in  the  current  work.  We  will,  therefore,  neglect 
contributions  to  the  pressure  in  Eq.  (14)  due  to  changes  in  entropy,  giving 
(to  second  order) 


P  2  Po  +  <£> 


S  ,p»p_ 


<«.>  ♦  i  <£> 


(P-P  > 


3p  S  ,p=p 


If  we  solve  Eq.  (15)  for  P  -  PQ  to  first  order,  we  obtain 


P  -  Po  -  AP  1  <£>  (p-po>  1  Co2*P 

S>P-P0 


or,  in  terms  of  the  "del”  operator. 


$P  -  c  2$p, 


where  Eq.  (4)  has  been  used  for  cQ. 


The  equation  of  momentum  conservation  to  first  order  is 
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(18) 


~^p  “  po  It  • 


Combining  Eqs.  (13),  (17),  and  (18)  yields  the  following  relation, 
accurate  to  second  order. 


<;•«)(«> 2-  -4  •  -i-  • 

c  at 
o 


(19) 


Combining  the  elementary  relation 


1  a2  o  ^  a2+  a*  2 

13.2  ♦  3  u  , 

5  - <T  (u  )  -  U  • - a-  +  (tt) 

2  3tz  at2  3t 


with  Eq .  (19)  gives 


,»  *.,*  *,  2  i  s2(u)2  i  ,a!.2 

(a.})(5.u)  .  - -  (sT)  . 

2c  <>t  C 


(20) 


(21) 


Also,  by  definition  of  the  O' Alembert lan  operator, 


2*2 

l  3  ii  v2*2  _2 »2 

~i  --r  ■ 7  u  -  °  u  • 

c 

o 


(22) 


Combining  Eqs.  (21)  and  (22)  gives 


2  1  t72*2  1  --2*2  1 

(u*v)(V*u)  =  j  7  u  -  jf  □  u - 2  ’ 


(23) 


We  now  consider  the  ($*u)2  terra  in  Eq .  (11).  By  using  first-order 
approximation  for  the  equation  of  continuity  [Eq.  (12) J,  we  obtain 

,X  *2  2  1  3Ap. 2 

(^*u) - J  ("St)  * 

po 


(24) 


Combining  Eq .  (24)  with  the  elementary  relationship 
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(25) 


L  a2(AP)2  _  aq  32aP  .  f3Ap  2 

J  IT2- "  -.7*  ‘  St> 

2 

and  with  the  D'Alembertian  operator  acting  on  (&p)  , 


Q  W  -  V2(A»)2  -  1  >2'-,0)2 


2  a,-2 

c  3 1 

o 


(26) 


gives 


<«)2  2- 1%.  tm2.  -  -v  o2(ap)2  -  -a?-  % 

p  2p  c  3t 

o  o  o 


(27) 


Substituting  Eqs.  (23)  and  (27)  into  Eq.  (11)  yields 


l  *  P^2«+2-  D2  £+ 


l.j  3X13XJ 


2P, 


1  .  3u. 2  Co  „2..  .2  A p  32Ap, 

2  <3F)  +  ~  2  7  (Ap)  "  2  “2}< 

Co  2Po  Po  3t 


(28) 


Using  the  elementary  fact  that 


j  V2(Ap)2  =  ApV2Ap  +  ($Ap)2, 


(29) 


we  obtain  for  the  last  two  terras  in  Eq.  (28) 


Co  1,2,.-  x2  Vp  32Vp 


2p 


~2  V*(Vp )*  - 


p  -  3t 

O 


Ap  V2Ap  +  -2—  ($Ap)2  -  ~~ 


Po  3t2 


(W  +  '>  ,,2ap  -  ‘  ^  , 

p  c  3t 

o  o 


(30) 
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substitution  corollary  allows  simpl If Ication  of  Eq.  (30)  to  give,  accurate  to 
second  order. 


V2(Ap)2 


Ap  a2Ap  2 


(31) 


Now,  Eqs.  (17)  and  (18)  may  be  combined  to  give  the  first-order 
approximation 


$Ap  » 


P  ,-*• 
o  3u 

2  3t  * 


c 

o 


(32) 


Finally,  Eqs.  (31)  and  (32)  may  be  combined  to  show  that  the  last  three  terms 
in  Eq.  (28)  vanish.  Hence,  Eq.  (28)  becomes 


i 

i.j 


3  <APUi Ujl  2  v2+2 

— 5 — V — =  P  V  u 

3x^3xj  o 


_2rl  +2 ,  l  2..  .2,  . 

O  ljPQo  +TCo  (Ap)  /Pq). 


(33) 


From  Eq.  (16)  we  have  the  first-order  approximation 

AP  i  cQ2 (Ap ) .  (34) 

Using  this  in  Eq,  (15)  yields,  accurate  to  second  order, 

2 

AP  -  co  Ap  +  i  co  (— j)  (AP)\  (33) 

3p  S,p=p 

O 

or 

Ap  -  co"2AP  l  -  i  cQ~6  (~y)  (AP)2.  (36) 

3p  S,p=pQ 

If  Eqs.  (33)  and  (36)  are  substituted  into  Eq.  (7),  the  following 
equation  is  obtained: 
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ahpi  3A«)!  + 

3p  S,P=P0  3t 


_2,  L  +2,1  2..  s2.  ,  „2  2 

□  I  2  PoU  +  2  Co  <AP)  /po!  "  poV  U  * 


An  equation  equivalent  to  Eq.  (37)  was  first  derived  by  Eckart  [15]  in 
1948,  and  later  by  Westervelt  [16]  in  1957.  In  1963,  Westervelt  [17]  trans¬ 
formed  this  equation  into  what  is  commonly  referred  to  as  the  simple-source 
formulation.  However,  in  so  doing,  he  used  an  expression  valid  only  for  plane 
waves.  In  what  follows,  Eq.  (37)  will  be  transformed  into  the  simple-source 
formulation  without  recourse  to  plane  wave  properties. 

In  order  to  carry  Eq.  (37)  further,  we  require  some  additional  first- 
order  relationships.  Since  it  is  irrotatlonal  to  first  order,  a  scalar 
function  <}»  known  as  the  velocity  potential  can  be  assumed  such  that  ,  accurate 
to  first  order, 


5  i  -$♦. 


The  scalar  potential  is  a  solution  to  the  linear  homogeneous  wave  equation 

0^<p  =  0. 

Now,  by  rearranging  the  definition  of  the  D’Alembertian  operator  acting 
on  u^  and  by  use  of  Eq.  (38)  we  obtain,  accurate  to  second  order: 


„2  2  2  2  2 

V  u  *  u  +  c 


N2- 


The  last  terra  can  be  rewritten  using  the  identity 


1 2  1  „2 .2  ,a2. 

I  •  j  v  ^  -$v$. 


From  the  definition  of  the  D’Alembertian  we  have 
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(41) 


02 .2  n2  2  .  -2 
7  (j>  *  □  <$>  +  cq  — ^ 

3 1 


a2  .2 


Also , 


V24> 


q2.  +  c  -2  3?i  _  c  -2  a2* 
°  at2  °  at2 


(42) 


since  q2(J>  =  0. 

Substituting  Eqs.  (40,  41,  &  42)  into  Eq.  (39),  the  following  second- 
order  expression  for  72u2  is  obtained: 


<,2  2  2  n2  2  -2  3 

V  u  =  □  u  +  c  — =r 

0  at2 


r1  „2,2  J.  -2  ,3*x2, 
Il°  *  +  co  (a7>  1 


(43) 


If  Eq.  (43)  is  substituted  into  Eq .  (37),  the  following  equation  results: 


□  2P  l  -  L  c  -ylt)  _  p  c  -2  [c  -2(31)2, 

2  °  ap2  s ,p=po  at2  ”  °  at2'  °  <3t>  1 


.  r-,2  .1  -l  2,.  .2  1  2  1  -2D4  ,.2., 

+  □  lj(>o  co(Ap)--pou  — 2  (♦  >1- 

O  t 


(44) 


We  next  wish  to  re-express  the  second  term  of  Eq .  (44)  in  terms  of  the 
acoustic  pressure.  We  begin  to  do  so  by  combining  Eqs.  (18)  and  (38)  to 
obtain  the  first-order  approximation 

$P  '  $  (PQ  |f>-  (43) 

Equation  (38),  which  defines  the  velocity  potential,  allows  a  certain 
freedom  in  the  choice  of  the  function  4> ,  since  it  is  only  this  function's 
gradient  that  is  therein  defined.  (This  is  analogous  to  the  freedom  of  choice 
of  gauge  in  electrodynamics.)  We  choose  to  impose  this  additional  freedom  in 
such  a  way  as  to  allow  Eq.  (45)  to  possess  the  "solution",  accurate  to  first 
order 
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(46) 


-  P  = 
o 


AP 


IF 


It  will  be  noted  that  Eq.  (46)  in  no  way  contradicts  Eq .  (45)  or 
F.q .  (38);  so  choosing  (J>  such  that  the  usual  additive  constant  obtained  in 
solving  an  equation  such  as  Eq.  (45)  to  be  equal  to  PQ  is  a  consistent,  and 
hence  permissible,  choice.  A  second  way  of  viewing  this  situation  is  to 
consider  the  combination  of  Eqs.  (38)  and  (46)  to  constitute  a  (consistent) 
definition  of  <£.  [Although  Eq.  (46)  is  only  one  of  an  infinity  of  possible 
choices  of  "gauge".]  Using  Eq.  (46)  in  Eq.  (44)  gives 


3p  S,P=PQ  3t 


,  _2,1  -1  2..  .2121  -2  3^  .  .  2. , 

+  D  l2  Po  Co  (Ap)  ”  2  pou  ~  1  poCo  7T  (*  )] 


o  2  o  o 


3t 


(47) 


Now  Introduce  the  perturbation  expansion*  P  -  Pq  =  eP^  +  e^P2  where 

2 

□  =  0  +  obtain  secondary  waves  as  solution  to  Eq.  (47). 

Equation  (47)  may  be  simplified  by  moving  the  terms  under  the 
O'Alembertian  operator  on  the  right-hand  side  of  this  equation  to  the  left- 
hand  side.  On  the  left-hand  side  of  this  equation  we  then  have  P2  +  addi¬ 
tional  terms  under  the  D' Alembertian  operator.  This  new  equation  can  now  be 
solved  for  this  new  combination,  subtracting  the  additional  terms  from  the 
solution  to  obtain  P2 •  In  practice,  the  terras  under  the  D'Alembertian 
operator  on  the  right-hand  side  of  this  equation  are  very  small  and  can 
actually  be  neglected.  In  any  case,  these  terms  will  clearly  be  non-growing 
contributions  to  the  solution  and  will  quickly  be  overwhelmed  by  the  growing 
contrlbut ions. 


*It  is  important  to  note  that  up  to  this  point  in  the  derivation  no  pertur¬ 
bation  analysis  has  been  used.  Hence,  Eq.  (47)  remains  valid  even  under 
conditions  that  invalidate  perturbation  analyses. 
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One ''further  remark  Is  worthwhile  in  discussing  the  D'Alembert  lan  terms  of 
Eq.  (  47).  If,  in  fact,  these  terras  are  not  negligible  in  comparison  to  the 
predicted  value  for  obtained  in  solving  these  equations,  the  result  thus 
predicted  will  most  likely  be  In  error.  This  is  due  to  the  fact  that  when  the 
D'Alembert lan  terras  are  lumped  onto  the  left-hand  side,  an  appropriate 
integral  term  must  be  included  (28]  to  reflect  the  fact  that  they  satisfy 
different  boundary  conditions  than  P2.  Hence,  in  solving  Eq.  (47)  via  this 
"lumping"  technique,  the  values  given  by  the  D' Alembert ian  terras  on  the  right- 
hand  sides  must  always  be  compared  with  the  predicted  value  for  P2  in  order  to 
Insure  consistency  of  l he  solution*.  This  fact  has  not  been  previously 
mentioned  in  the  literature. 


At  this  point,  we  assume  the  D'Alembertlan  terms  on  the  right-hand  side 
of  Eq.  (47)  are  negligible. 

If  we  make  the  definitions: 


B/A  =  ~  (^) 

c  3p  S,P=P 
o  o 


r  =  nonlinearity  parameter 


l  + 


_B 

2A 


q  =  simple  source  strength 


-in-4  i>7<pi 

P0  Co 


we  can  cast  Eq.  (47)  into  the  form 


CTP 


2P  '  2  _c  .3c 


-po 


(48) 


*The  calculations  relevant  to  the  problems  considered  in  this  thesis  are 
performed  In  Chapter  IV.  B,  where  estimates  of  the  errors  introduced  In  the 
relevant  surface  and  volume  integrals  are  discussed. 


The  prime  has  been  added  to  the  symbol  for  the  second-order  pressure  to  denote 
that  certain  second-order  quantities  have  actually  been  neglected*. 

E.  Plane-Wave  Scattering  from  a  Vibrating  Planar  Surface 

In  this  section,  the  problem  of  the  generation  of  sum-  and  difference- 
frequency  waves  will  be  addressed  for  the  case  of  a  plane  wave  normally 
incident  on  a  surface  deforming  uniformly  and  harmonically. 

L .  Censor-Method  Solution 

Although  he  treated  several  different  geometries,  Censor  never  considered 
the  simplest  possible  case.  This  case  is  a  plane  wave  of  angular  frequency  ui 
normally  incident  on  an  infinite  plane  vibrating  uniformly  with  angular 
frequency  SI  (see  Fig.  2).  We  obtain  the  solution  to  this  problem  following 
the  procedure  presented  by  Censor  [23]. 

We  represent  the  incident  plane-wave  acoustic  pressure  as: 

Ut  =  Pte~la,((x/co)+tl  ,  (49) 

where  x  is  the  Eulerian  position  coordinate**.  The  scattered  wave  is  assumed 
to  be  of  the  form: 


U 


s 


=  /+<”dva(v)e1Vl(x/co)"tl. 


(50) 


U|+Ug  is  the  total  wave  field.  Using  the  method  of 
the  vanishing  of  the  normal  particle  displacement  at  the 


Censor,  we  require 
planar  surface  x  =  0 


*In  solving  Eq.  (48)  in  this  report,  we  will  not,  of  course,  obtain  a  solution 
accurate  to  second  order.  Pn'  simply  reflects  corrections  to  the  primary 
field  resulting  from  the  "enhancement"  factors  discussed  in  Chapter  II, 

Section  C,  arising  from  second-order  quantities. 


**Censor  does  not  specify  that  x  is  the  Eulerian  position  coordinate,  but  this 
choice  appears  consistent  with  the  way  he  treats  this  quantity  in  his  paper. 
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FREQUENCY  Q 


— ►  RAOIATED  PLANE  WAVE  OF 
ANGULAR  FREQUENCY  fl 

— ►  SCATTERED  WAVE  OF  ANGULAR 
FREQUENCY  u 

—  INCIDENT  PLANE  WAVE  OF 
ANGULAR  FREQUENCY  ut 


Fig.  2  -  Geometry  of  plane-wave  scattering  from  a  vibrating  plane 


for  the  rigid-body  scattering  problem.  Since  we  consider  only  linear  wave 
fields  in  the  Censor  method,  the  acoustic  pressure  satisfies  the  equation 


3P  3^x 

'  "  P°  at2  * 


For  assumed  plane  waves,  we  obtain 


(51) 


x  -  PQ(v-2  P. 


(52) 


Hence,  the  operator  0  used  by  Censor  becomes 


2 


and  Censor’s  Eq.  (4)  for  the  boundary  condition  becomes 


(53) 


J 
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{_  p  e-iwl(x/co)+tl 

Co 


+ao 

+  /  dva( v) 

— oo 


_w_1eiu[(x/co)-t ] j 
c 

o 


0. 


(54) 


We  now  let  x  =  tsinflt  on  the  surface  and  perform  the  following  expansions 

±(lve/co)sinflt  B  L  ±  t  *£  sinSit  +  ...  (55) 

c 


a(v)  =  A(v)  +  e3(v)  +  •  ••  •  (56) 

We  substitute  Eqs.  (55)  and  (56)  Into  Eq .  (54),  retaining  terms  to  zeroth 
order  In  e: 

_  "lp<e"iU>t  +  f  dvA(v)  — "le"lVt  =  o.  (57) 

Cl  ^ 

o  -°°  O 

Multiplying  Eq.  (57)  by  eiV  £,  integrating  over  t,  and  letting  v'*v 
yields 


A(v)  =  pj6(v-w).  (58) 

Taking  next  the  terras  of  first  order  In  e  when  Eq.  (55)  is  substituted 
into  Eq.  (54),  we  obtain 

2  ip  +00  —1 

- j  sinftt  +  /  dvB(v)  ~  e  iVt  -  0.  (59) 

c  -•  o 

o 

Solving  Eq.  (59)  for  B(v)  gives 

B(v)  ■  — -  (d(v-B-oj)  -  6(v-Hl-w)J.  (60) 

o 
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(65) 


K  -  c 


2 

*  3x2 


,32C, 


-cQ  <1+r><;£><7i) 

3x 


[In  this  section  x  refers  to  a  Lagrangian  coordinate-] 

2 

We  now  represent  f,  in  a  perturbation  series  F,  *  cf,^  +  e  +  •••>  whore 
c  is  of  the  order  of  the  Mach  number. 


We  then  substitute  into  Eq .  (65)  and  equate  equal  orders  of  e. 
provides  the  first-order  equation 


3x 


2 


=  0, 


This 


(66) 


(the  equation  of  linear  acoustics),  and  the  second-order  equation 


_a 

*>  3x  1  J  ' 


3x 


(67) 


To  solve  the  problem  of  nonlinear  scattering  of  a  plane  wave  by  a  vibrat¬ 
ing  plane  using  the  perturbation  approach  outlined  in  Section  II.  C,  we  must 
first  solve  the  first-order  (linear)  Eq.  (66).  The  physical  boundary  condi¬ 
tion  to  be  met  is  that  there  is  no  relative  displacement  between  the  planar 
surface  and  the  fluid  particles  in  contact  with  the  surface*.  This  condition 
can  be  naturally  met  in  Lagrangian  coordinates  by  equating  the  displacement  of 
a  particle  at  the  surface  to  the  displacement  of  the  surface  itself. 

The  first-order  solution  is  clearly  the  sum  of  the  incident,  scattered, 
and  radiated  waves: 


C(l)  -  CjtsinmCt  +  ^)  -  sin<v(t  -  ^-)]  +  f^sin^t  -  ^) . 


(68) 


*Since  we  are  solving  second-order  equations  here,  there  should  actually  be  a 
second  boundary  condition.  In  this  case  this  condition  takes  the  form  of 
requiring  that,  apart  from  the  quantity  C^sinutlt  +  (x/cQ)],  there  are  no 
Incoming  waves. 
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The  first  two  terms  In  this  expression  represent  the  particle 
displacement  for  the  rigid-body  scattering  solution,  in  which  the  reflected 
amplitude  is  the  negative  of  the  incident  amplitude  to  insure  the  vanishing  of 
the  displacement  at  the  rigid-body  surface.  The  third  term  represents  the 
radiated  wave. 

The  first-order  solution  is  now  substituted  into  the  right-hand  side 
of  the  second-order  Eq.  (67)  and  the  resulting  linear  inhomogeneous  equation 
solved  for  f^.  If  we  retain  only  those  terms  that  contribute  to  to^ ,  we 
obta in: 

Aim  H+rHA" 

4(2)±  "  Co  2  2c 

dX  O 

to_x 

+  X 

f-u)_sin(oj,  t  + - )  +  o>  sinu>  (t  -  — )] 

o  o 


where  again  <u  =  u>±fl  (and  w>tt  is  assumed). 

X 

Equation  (69)  may  be  solved  by  the  usual  methods  for  ordinary 
inhomogeneous  equations  with  constant  coefficients.  The  result  is: 


-(i+r)f,,f.  wil 

.  irrx  ,  x. 

f’(2)±  “ - 2c - {2“  cos(V  ~  "t  c~> 

o  o  o 


+  — 2 - 2  {-sin(w±t-hd_  ~)  +  sinCw^t-w^  ~) ) } 


u>  -to_ 
±  + 


where  a  homogeneous  solution  has  been  added  that  causes  =  0  at  x  =  0. 

We  obtain  the  corresponding  acoustic  pressure  by  performing  the 
integrat ion 

P  »  -  [  3xo  _ LI)* 

P(2)±  >  dxpo  3t2  ’ 


since  this  relationship  is  exact  to  all  orders  in  Lagrangian  coordinates. 
Also,  using  the  fact  that  the  acoustic  pressure  and  displacement  amplitudes 
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(for  plane  waves)  are  related  to  first  order  by  5  *  P/(pQmc0),  we  obtain  the 
following  expression  for  the  sum-  and  difference-frequency  pressure  waves 


P(2)± 


(i+r)p  p  u>  x 

- sin(u)  t-u>  ~)  -  cos(u>  t-0)  -~) 

,  ^  2  c  ±±c  ±tc 

4p  c  o  o  o 

o  o 


u>_ 

[cos(aJ,  t+m_  — )  +  —  cos(w.t-u),  — )]}. 
1  ±  +  c  w  ±  ±  c  ' 1 

o  ±  o 


(71) 


It  should  be  noted  that  Eq .  (71)  could  have  been  stated  in  dimensionless 

2 

form.  Specifically,  the  ratio  P~ , / { f ( l+r)P . P  ]/(p  c  )}  depends  only  on 

zx  I  r  o  o 

x  *  w  t,  C,  *  (a  x/c  ,  the  frequency  ratio,  and  the  boundary  condition. 

±  x  o 

Comparing  Eq.  (71)  to  Censor's  solution  [Eq.  (64)],  the  most  striking 
distinction  between  them  is  the  presence  of  the  "x"  coefficient  in  the  first 
term  inside  the  bracket  of  Eq.  (71).  Censor's  result,  being  a  boundary-effect 
solution,  does  not  grow  with  distance  from  the  scattering  surface.  Contribu¬ 
tions  fror’  medium  nonlinearities,  being  a  cumulative  volume  effect,  do  grow 
with  distance  from  the  scattering  surface.  Hence,  the  nonlinear  effect  pre¬ 
dicted  by  Eq.  (71)  will  overwhelm  that  predicted  by  Censor  within  a  small 
distance  from  the  scatterer. 

At  this  point  we  put  these  remarks  on  a  more  quantitative  basis  as  well 

as  calculate  a  region  of  validity  for  Eq .  (71).  First  we  consider  relative 

contribution  of  the  terra  arising  from  satisfying  the  boundary  condition  to  the 

term  representing  the  growing  contribution  from  the  virtual  volume  sources 

(i.e.,  the  term  with  the  "x”  coefficient).  The  term  that  arises  from 

satisfying  the  boundary  condition  is  the  last  term  in  Eq.  (71).  Hence,  a 

quantitative  estimate  of  the  relative  contributions  can  be  obtained  via  the 

2  2 

ratio  of  the  coefficients:  ( 2u)tcq)/ [  (id^  -  w_  )x].  Using  primary  frequencies 

of  160  and  100  kHz  (which  are  of  experimental  interest  later  in  this  thesis), 

.  *> 

this  ratio  gives  approximately  (1.94x10  m)/x  in  the  difference-frequency 
case.  Hence,  the  boundary  effect  becomes  less  than  8%  of  the  volume  effect  at 
just  one  difference-frequency  wavelength  distance  from  the  scattering  surface 
(this  wavelength  is  approximately  2.5  cm). 


Also  of  Interest  In  this  problem  Is  the  distance  to  which  the  solution 
represented  by  Eq.  (71)  remains  valid.  This  can  be  estimated  by  comparing  the 
energy  density  of  the  secondary  waves  to  the  energy  density  of  the 
primaries.  The  "secondary  waves"  include  not  only  the  sum-  and  difference- 
frequency  waves  but  the  second  harmonic  waves  as  well. 

We  can  estimate  the  second-harmonic  pressure  using  the  formula  obtained 

by  Fubini-Ghi  ron  [7]:  =  [  ( Po^u>)/(4p  )  1  [  2+(B/A)x  ]  ,  where  the  parameter 

Y  of  Fubinl-Ghiron's  original  expression  has  been  replaced  by  1  +  (B/A)  (see 

Reference  5,  p.  99).  We  can  estimate  the  sum-  and  difference-frequency 

pressures  using  the  coefficient  of  "x”  in  F.q .  (71).  Using  100-kHz  primary  of 

lO-’-Pa  amplitude,  Fubini-Ghlron's  formula  reduces  to  p  =  (744  Pa/m)x. 

Similarly,  if  the  second  primary  is  taken  to  be  of  160-kHz  frequency  and  also 

10^-Pa  amplitude,  the  sum-  and  dif ference-f requency  pressures  are 
3  3 

P+  =  (5.46x10  Pa/m)x  and  P_  =  (1.26x10  Pa/m)x.  We  can  estimate  the  energy 

densities  of  each  of  the  relevant  waves  using  the  elementary  plane-wave  energy 
2  2 

density  formula  P  /(2p  c  )•  We  estimate  the  energy  density  of  the  primaries 

by  inserting  10^  Pa  for  VQ  and  multiplying  by  3  (to  account  for  the  incident, 

reflected,  and  radiated  waves),  giving  approximately  6.7  J/m  for  the 

primaries.  We  estimate  the  energy  density  of  the  secondary  waves  by  applying 

this  equation  separately  to  each  of  the  four  secondary  waves  In  turn  and 

-3  5  2 

adding.  This  gives  approximately  ]7xU)  (J/m  )jx  .  Hence,  the  energy  density 
of  the  secondary  waves  becomes  1%  of  the  energy  density  of  the  primaries  at  an 
approximate  distance  of  x  =  3.1  m.  Therefore,  it  is  reasonable  to  expect  the 

solution  to  be  reliable  out  to  a  distance  of  3  m. 

Finally,  to  establish  the  fact  that  viscous  terms  do  not  become  important 
prior  to  this  distance,  we  calculate  the  discontinuity  distance  (see  Reference 
5,  p.  104]:  l/i  =  [  (  1+(B/2A)  ]  [  (u>uo/co^)  ] .  For  a  frequency  of  160  kHz  and 
amplitude  of  10^  Pa,  this  formula  gives  approximately  9.6  m.  Hence,  if  the 
solution  Is  restricted  to  distances  less  than  3  m,  viscous  terms  can 
reasonably  be  expected  to  play  a  minor  role. 

One  ftnal  comparison  of  Interest  is  the  volume  effect  term  of  F,q .  (71) 
relative  to  the  boundary-effect  predicted  by  Censor  as  represented  by 
Kq .  (64).  We  represent  this  as 
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mOim 


^  ^i^Censor 


[P±] 


nonlinear  volume 


4c 

o 

8(1 l+f)x* 


For  a  typical  planar  surface  frequency  of  100  kHz,  the  two  effects  become 
equal  at  a  distance  x  on  the  order  of  0.1  cm.  Due  to  the  presence  of  the  "x" 
factor  in  the  nonlinear  volume  term,  these  volume  effects  continue  to  grow 
from  this  point,  while  Censor's  surface  effect  remains  constant. 

One  should  not  be  disturbed  by  the  presence  of  the  11  term  In  the  denom¬ 
inator  of  this  ratio.  In  the  case  of  low  SI,  the  relevant  factor  to  scale  the 
distance  Is  the  wavelength  associated  with  Si,  which  is  cQ/Sl.  Hence,  if  we  let 
x  =  fcQ/0,  we  can  determine  the  fraction  (f)  of  a  wavelength  at  which  the  non¬ 
linear  volume  effect  overtakes  Censor's  surface  effect.  This  occurs  for  x  = 
0.89  (c0/fl).  Therefore,  even  in  the  limiting  case  in  which  the  frequency  of 
vibration  of  the  planar  surface  approaches  zero  (maximizing  Censor's  effect 
relative  to  the  nonlinear  volume  effect),  the  difference-frequency  pressure 
generated  by  the  fluid  medium  exceeds  that  produced  by  Censor's  surface  effect 
within  a  propagation  distance  less  than  the  longest  wavelength  involved  in  the 
problem. 


3.  Some  Comments  Regarding  the  Censor  Approach  to  the  Problem  of  the  Plane 

Censor  states  In  Section  (4)  of  his  paper  that  the  fundamental  boundary 
condition  for  the  problem  is  the  vanishing  of  the  normal  displacement.  It  is 
interesting  to  note  that  such  a  treatment  is  equivalent  to  simply  recasting 
the  incident  wave  into  Lagrangian  coordinates  (a  system  that  follows  fluid 
motion),  treating  the  problem  as  a  simple  rigid-body  scattering,  and  then 
transforming  the  result  back  into  Eulerian  or  fixed  coordinates.  We  proceed 
to  demonstrate  the  validity  of  this  interpretation  of  Censor's  approach. 

We  begin  by  transforming  Censor's  incident  plane  wave  (apparently  written 
in  Eulerian  form)  into  Lagrangian  coordinates.  An  arbitrary  Eulerian 
function,  fE(x,t),  may  be  transformed  into  its  associated  Lagrangian  function, 
fL(a,t),  by  an  expansion  of  the  form 
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f  (a  ,t)  =  f  *(x,t ) 


x=a+f. 


<.,o  ' 


x=a 


+  C(a,t)  +  ...  , 

3x  x=a 


where  a  =  Lagranglan  position  coordinate 
x  =  Eulerian  position  coordinate 
C  =  displacement  (common  to  both  systems). 


(72) 


Letting  fE(x,t)  *  PE(x,t)  =  PQe  iwl (x/c)+t ] ^  which  is 
Censor's  incident  plane  wave,  and  performing  the  expansion 
boundary  surface  of  the  plane),  where  £($,t)  *  esinflt,  one 
terras  of  order  e2  or  greater) 


the  form  of 
around  a  =  0  (the 
obtains  (neglecting 


PL(o,t) 


P  e 
o 


-iwt 


2c  Po[e  "e  1* 

o 


(73) 


To  compare  Eq .  (73)  to  Censor's  result  (Eqs.  63),  we  must  re-express 
Eq .  (73)  in  Eulerian  coordinates.  This  may  be  done  by  constructing  the 
function  PE(a,t),  where  the  Lagrangian  coordinate  "a"  is  inserted  into  the 
right-hand  side  of  Eq.  (73)  in  the  appropriate  places  to  form  outgoing  plane 
waves.  One  can  then  expand  the  resulting  function  PL(a,t)  in  a  series  of  the 
form 


fE(x,t) 


fL(a,t) 


a»x-£(x  ,t) 


fL(a,t) 


a=x 


a=x 


S(x,t). 


(74) 


We  construct  PE(a,t)  from  the  form  of  Eq .  (73)  as 
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PL(a,t) 


P  ei«[(a/co)-t] 


“ePo  ,  i(u)-ft)  [  (a/c  )-t  ] 

~2c~  {e  0 

o 

i(u-Kl)I(a/c  )-t  J, 

C  o  /  I 


(75) 


and  now  expand  PL(a,t)  in  a  series  of  the  form  of  Eq.  (74)  (neglecting  terms 
of  order  e2  or  greater)  obtaining: 


E.  .  io)f(x/c  )-tl 

P  (x,t)  -  PQe  1  o'  1 


a>ePo  {ei(u)-$))[(x/co)-t]  _  ei(artO)[(*/co)-t]} 

2c 

o 

fiaiPo  i»((x/c  )-t]  l(u)-fl)  i(a>-0)[(x/c  )-t). 

e  °  2c  e  °  J 

o  o  o 


Hiit«).ei(0,+n)[(x/co)'tl>C(x>t). 


(76) 


We  now  let  £  »  esinflt  and  again  neglect  terms  of  order  e2  or  higher.  The 
result  is 

PE(x,t)  =•  P 

o 


_  Tl°  {ei(«o-fl)I(x/c)-t)  _el(oi+fl)[(x/co)+t]} 

2c 

o 

_  CU>Po  {ei[(«x/c0)-(«-a)tJ  _  ei{(wx/co)-(u)+ft)t]  J  ^  (7 

2c 

o 

We  now  recognize  that  the  above  expansion  is  valid  only  at  x  =  0  (since 
this  is  the  only  place  where  ?  ■  esinflt).  Evaluating  Eq.  (77)  for  x  **  0  and 
using  the  fact  that  e  ■  Pr/Poco0,  we  obtain 
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(78) 


E.  _  -itot  PoPr“  ,  -i(oi-n)t  -i(to+Q)t, 

P  (o,t)  =  P  e - j-  [e  v  "  e  v  ] 

P  c  « 
o  o 

Comparing  the  sum-  and  difference-frequency  components  of  Eq .  (78)  with 
Censor’s  result  (Eq.  (63)],  we  note  that  for  il<< a),  they  are  identical  [P  of 
Eq.  (78)  is  equivalent  to  of  Eq.  (63)]. 

Hence,  it  is  seen  that  Censor's  result  is  clearly  of  the  same  order  as 
the  difference  between  Eulerian  and  Lagrangian  coordinates.  Such  effects 
normally  are  not  even  experimentally  measurable,  since  no  presently  available 
measurement  hydrophone  is  either  completely  rigid  (and  hence  measures  in 
Eulerian  coordinates)  or  moves  completely  freely  with  the  fluid  (and  hence 
measures  in  Lagrangian  coordinates).  Any  presently  available  hydrophone  will 
have  an  uncertainty  in  its  measuring  capability,  in  an  experiment  design-ed  to 
measure  difference-frequency  waves,  of  the  order  of  the  difference  between  the 
pressure  predicted  in  a  Lagrangian  frame  and  the  pressure  predic-ted  in  an 
Eulerian  frame.  This  difference  is  known  as  “ pseudosound''  and  is  treated  more 
fully  in  the  section  describing  the  experimental  results.  Hence,  the  effect 
predicted  Ly  Censor  cannot  be  measured  with  present-day  technology. 

F.  Cylindrical-Wave  Scattering  from  a  Vibrating  Cylindrical  Surface 

The  present  section  considers  the  problem  of  the  generation  of  sum-  and 
difference-frequency  waves  when  a  cylindrical  wave  (of  angular  frequency  w") 
is  normally  incident  on  a  cylinder  whose  surface  deforms  radially  in  a  uniform 
and  harmonic  fashion.  It  is  assumed  that  the  waves  at  frequency  u>  =  u>"  ±  u)' 
are  outwardly  propagating  waves  in  the  limit  r+®. 

I .  Solution  Using  Second-Order  Nonlinear  Wave  Equation. 

We  now  consider  the  problem  of  a  vibrating  cylindrical  surface.  The 
geometry  of  the  problem  Is  indicated  in  Fig.  3. 
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In  solving  the  second-order,  nonlinear  wave  Eq.  (48),  we  assume  the 

scattered  pressure  P„  and  the  simple-source  term  p  (3q/3t)  may  be  represented 

s  o  ' 

as*: 


p(2)  “  Re  I  Pn(r,9)e  iuJnC 


(79) 


and 

PQ  *  “Re  l  Bn(r,0)e  ^n11 ,  (80) 


VIBRATING  RADIALLY 
AT  ANGULAR  FREQUENCY  cj7 


Fig.  3  -  Geometry  of  cylindrical  wave  scattering 
from  a  vibrating  cylinder 


♦Specific  estimates  of  errors  introduced  into  the  solution  of  Eq.  (48)  by  the 
neglect  of  the  D'Alembertian  terms  of  Eq.  (47)  are  provided  in  Chapter  IV.  B. 


where  the  subscript  n  is  used  to  distinguish  between  sum-  and  difference- 
frequency  components:  n  ■  1  refers  to  the  sum  frequency  and  n  =  2  refers  to 
the  difference  frequency.  Also  present  in  this  expansion  are  terms  corres¬ 
ponding  to  the  harmonics  of  the  primaries;  but  since  only  the  sum  and 
difference  frequencies  are  of  interest  here,  no  special  notation  will  be 
provided  for  these  terms.  At  the  conclusion  of  this  analysis,  the  numerical 
subscripts  will  be  replaced  by  the  more  descriptive  “+"  and  ”- "notation* 

Hence,  for  example,  P^(r  ,9)+P+(r)  and  P2(r ,9)+P_(r).  These  replacements  will 
also  help  avoid  confusion  between  P2(r,0)  (the  difference-frequency  pressure) 
and  P2(r)  (the  perturbation  solution  accurate  to  second  order). 

By  substituting  Eqs.  (79)  and  (80)  into  Eq.  (48),  we  obtain  the  following 
equation  for  the  time-independent  amplitudes 


V2?  + 


+  k  P  =  B 
a  n  a  n 


with  kn  =  mn/cQ. 

We  define  an  associated  Green's  function  gn(r,r')  such  that 

^2gn(  t  ,r '  )  +  kn2gn(r,r')  -  -<S(r-r')  (82) 

A  -►  -► 

subject  to  the  boundary  condition  r  vg  (r,r')  =  0  on  the  cylindrical  surface, 

a  n 

where  r  is  a  unit  normal  vector,  directed  outward  from  that  surface.  Here 
<5(r-r')  is  the  Dirac  delta  function.  We  also  have  the  condition  that  as  r 


gn  ~ 


lknr/rV2 


)  x  function  of  (r't  0). 


A  representation  of  Pn  may  be  obtained  by  multiplying  Eq.  (82)  by  Pn  and 
Eq.  (81)  by  gn(r,r'),  subtracting  the  resulting  equations,  and  Integrating 
over  primed  variables.  The  result  is: 


Pn(r*) 


-  / dt'B  <r')g  fr.r')  +  //  dS'  $'P  (r ' ,9 ' )g  (r  ’ ,r) 


In  which  the  vanishing  of  the  normal  gradient  of  the  Green's  function  has  been 
Imposed.  Formally,  the  volume  Integral  of  Eq.  (83)  is  taken  over  all  space, 
excluding  the  volume  interior  to  the  cylinder  vibrating  at  angular  frequency 
a)'.  The  surface  integral  in  this  equation  represents  integrals  over  the 
surfaces  of  both  cylinders  as  well  as  over  the  surface  at  infinity.  We  will 
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consider  the  surface  integral  over  the  <a"  cylinder  later. 

In  carrying  Eq.  (83)  to  a  final  solution,  the  surface  integral  term  will 
actually  be  dropped,  and  the  volume  integral  will  be  analyzed  only  between  r' 

=  a  and  r'  =  r  (justification  for  this  will  be  presented  shortly  as  well  as  in 
Chapter  IV).  In  summary,  what  this  means  is  that:  1)  the  surface  integral  at 
r '  =  a  is  neglected,  2)  the  surface  integral  over  the  u>"  cylinder  is 
neglected,  and  3)  the  volume  integral  from  r'  =  r  to  r*  =  00  is  neglected. 
Estimates  of  the  errors  arising  from  some  of  the  neglected  terms  will  be  made 
presently  (the  rest  being  postponed  until  Chapter  IV,  Section  B) . 

We  now  consider  the  surface  integral  in  Eq .  (83)  in  somewhat  greater 
detail.  In  the  current  problem  we  are  considering  rigid-body  scattering, 
although  the  surface  is  permitted  to  deform  harmonically  at  frequency  u>’.  The 
requirement  of  rigid-body  scattering  manifests  itself  in  the  handling  of  this 
surface  integral.  From  the  first-order  equation  of  momentum  conservation  (Eq. 
(18)],  it  is  clear  that  for  harmonic  time  dependence,  the  gradient  of  the 
pressure  field  is  proportional  to  the  velocity  of  the  fluid  (and  hence  the 
velocity  of  the  surface).  Therefore,  in  the  surface  integral  of  Eq .  (83)  the 
"V'P  "  term  may  be  viewed  as  the  component  of  the  surface  velocity  at  the  sum 
and  difference  frequencies*.  We  interpret  the  "rigid-body  oscillation"  of  the 
surface  as  constraining  the  surface  to  vibrate  only  at  the  frequency  at  which 
it  is  being  driven  (i.e.,  m').  Hence  this  term  and  the  surface  integral  of 
Eq.  (83)  vanish. 

We  next  consider  what  influences  viscosity  might  have  on  this  surface 
integral  (and  hence  on  the  sum-  and  difference-frequency  pressures).  We  have 
already  demonstrated  in  Section  II.  C  that  effects  of  viscosity  on  the  propa¬ 
gation  of  acoustic  waves  are  insignificant  in  this  type  of  problem.  Hence,  it 
is  reasonable  to  separately  consider  a  "boundary  region"  and  a  "propagation 


*This  statement  is  based  on  the  fact  that  the  equation  VP  =  -p0(3u/3t)  is 
correct  to  all  orders  (in  the  lossless  case)  in  Lagrangian  coordinates.  In 
Eq.  (83),  of  course,  the  expression  V'Pn  is  evaluated  in  Eulerian 
coordinates.  However,  the  difference  between  evaluating  a  function  in  these 
two  reference  frames  is  of  the  order  of  pseudosound.  The  validity  of  this 
statement,  as  well  as  an  estimate  of  the  effect  of  an  error  of  the  order  of 
pseudosound  on  the  surface  Integral  of  Eq.  (83)  is  provided  in  Chapter  IV.  B. 


region"  (see,  for  example,  the  discussion  on  pages  281-286  of  Reference  29) 

For  a  frequency  of  100  kHz  in  water,  this  boundary  layer  is  of  approximately 
5. 6-microns  thickness  (using  equation  6.4.31  of  Reference  29).  We  assume  that 
outside  this  boundary  layer,  the  lossless  equations  apply  and,  hence,  the 
Green's  function  solution  represented  in  Eq.  (83)  is  appropriate.  We  note 
that  only  the  normal  velocity  component  is  present  In  the  surface  integral . 
Hence,  any  tangential  velocity  component  arising  from  viscous  boundary-layer 
effects  will  not  influence  the  radiation  field  in  a  significant  way. 


It  Is  assumed,  then,  that  only  the  volume  integral  from  r'  =  a  to  r*  •=  r 
contributes  significantly  to  the  solution.  At  this  point,  we  neglect  the 
surface  integral  of  Eq.  (83)  and  simply  represent  Pn  as 


P>V  -  I  dT'VJ')8n<E';'> 


(84) 


(volume  integration  only  between  r*  *  a  and  r’  =  r) 


The  rigid-body  Green's  function  appropriate  to  cylindrical  geometry  is 
required  in  Eq .  (84).  It  is  well  known  (42jXaTsd  is  given  by 


8n(r,?')  l  (2-6mo)cosm(*-0')  [Hm(  1  >(knr ) Jm(knr '  )\ 
m=o 


ra  n  m  n 


~77T”~ —  H  r)H  (1)(k  r’>l 

H  (l)'(k  a)  m  n  m 
m  n 


r>r  ’ 


(83) 


where  &mQ  is  the  Kronecker  delta.  Also,  gn(r,r')  =■  gn(r’,r). 

The  Green's  function  gn(r,r')  is  not  needed  for  the  region  r<r'  since 
contributions  from  this  region  (which  is  further  from  the  source  than  the 
point  of  interest)  tend  to  phase  cancel  against  one  another  and,  hence, 
contribute  very  little  to  the  overall  pressure  at  the  observation  point  r. 
Contributions  from  the  region  r>r',  however,  tend  to  add  constructively  and 
give  by  far  the  majority  of  the  pressure  at  the  observation  point.  Hence,  in 
both  the  current  geometry  and  the  one  In  the  next  section  (Involving  a  plane 
wave  incident  on  a  vibrating  cylinder),  the  contributions  to  the  Green's 
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function  Integral  between  the  observation  point  r  and  infinity  will  be  assumed 
negligible . 

To  complete  the  solution  to  the  current  problem,  we  must  now  calculate 
the  functions  Bn(r').  In  order  to  do  so,  we  refer  once  again  to  Fig.  2.  From 
elementary  trigonometry  we  have  r*  =*  (r2  +  b2  -  2rbcos0)  ^2  .  We  now  use  the 
summation  theorem  (see,  for  example,  Gradshteyn  and  Ryzhik  [A3])  to  obtain 


H  <1>(“^1) 
o  c 

o 


H  <1>(‘^)j  (‘HJL) 
o  c  ov  c 
o  o 


+  2  l  H  ^ 1  J  (  ~r~)cosme 

,  m  c  me 
m=l  o  o 


(86) 


for  r<b,  and 


H  ^(^  =  J(^)H o(1)(^) 
o  o  o 


+  2  l  Jm(-^)Hm(l)(^)cosm0 

o 


(87) 


,  m  e  ra 
m=l  o 


for  r>b. 

The  incident  wave  is  assumed  to  be  of  the  form  of  a  uniformly  diverging 
cylindrical  wave  originating  at  the  surface  of  the  cylinder  vibrating  at  an 
angular  frequency  u>"  (the  interiors  of  the  cylinders  are,  of  course,  excluded 
from  the  region  of  interest).  Hence,  allowing  A'  to  be  a  pressure  amplitude, 
the  incident  cylindrical  wave  may  be  written  in  the  form  A' 1  )(u>"r  ’ /cQ) . 
Equation  (86)  is  used  to  re-express  Also  note  that  Pinc  is  singular  at 

Vhtvorlgin  of  the  u>"  cylinder. 

The^e<re8pondlng  particle  velocity  may  be  obtained  by  use  of  the  first- 
order  equation  df^momentum  conservation  [Eq.  (18)],  expressed  in  the  form 
V  »  [-i/p0«] [ 3P/3r ]  (for  the  radial  component).  We  thus  have  for  the  incident 
particle  velocity  wave, 
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iA’  H  (U(w^b)J  (W_£) 


inc  pc  o  c  1  c 
o  o  o  o 


+  JA L  l  Coa(mO)lJ  (— ) 

n  ~  L  m+I  c 


P  c  , 
o  o  m»l 


o  o 


(88) 


V  »  al,  (DX'r,  “iuTt 
)  A  cosra0H  ( - )e 

m=0 


We  assume  the  first-order  scattered  acoustic  pressure  wave  to  be  of  the  form: 

(89) 

Hence,  the  scattered  particle  velocity  wave  is 

0  -  (  As.  „  <»,!£.,  +  ~ —  l  A  cosm9 

scatt  pc  o  c  2pc  +  m 
o  o  o  o  o  m*l 


[H  4.,(1)(— )  -  H  ^'r-L)]}e 
1  ro+1  c  m-1  c 


(l>,<*>"rxl,  -iw"t 


(90) 


We  assume  the  scattering  cylinder  to  be  rigid  in  solving  the  first-order 
problem;  hence,  the  boundary  condition  becomes  Ujnc  *  ~^scatt  at  r  =  a'  This 
condition  yields  the  scattering  amplitudes 


A 
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H  <l>(“-^) 

_A*  I _ Io_  j  (WJa 

(l)^  V  cq> 

I  C 
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-2A'H  <l>(^) 


'  Vl^>J 
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IH 
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^  -  v,'(i,e» 
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m  >  1 . 


(91) 
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The  first-order  scattered  acoustic  pressure  for  r<b  Is  thus  given  by  Eq . 
(89)  with  Eq.  (91)  substituted  for  A^- 

The  Incident  acoustic  pressure  for  r>b  is  given  by: 


P.  “  (A'J  (— )H  (1>(~) 
Inc  1  o  c  o  c 
o  o 


+  2  £  J  (— ^)H  ^  \ - )cosm0Jxe 

L.  m c  '  m  c  * 

m=l  o  o 


(92) 


Hence,  the  general  solution  to  the  rigid-body  scattering  problem  to  first 
order  is 


P  .  (A' H  <1>(^)J  (2^) 
tot  o  c  ov  c 

o  o 


+  A  h  (1)(^)  +  l  [A  h  ^\2JL) 
oo  c  '  Ly  1  m  m  c  ' 
o  m=l  o 


(i)/“”r> 


and , 


+2A’H  (1)(— )J  (— )]cosm0}e'it‘’  C 
m  c  m  c  ‘ 
o  o 

r<b 


P  =  {A-j  (^)h  (1)(^) 

Lot  ov  c  7  o  v  c  7 

o  o 


+  B  H  (1)(— )  +  l  IB  H  (1)(— ) 
oo  c  1  m  ra  c 

o  ra«l  o 


(93) 


+2A’J  (— )H  (1)(—  )]cosra0}e“iUI  * 
mem  c  1 
o  o 

r>b. 


(94) 


This  can  be  rewritten  as 
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/]%  w  r>  m  r< 

P_  „  =*  {A’H  A  '( - -)J  ( - -) 

tot  o  c  o  c 

o  o 


+  A  „  <l>(i£r  +  $  [A  H  (1)(^) 

oo  c  lmm  c 

o  m=*l  o 


a)  rv  ui  r 


+2A'H  (1)( - -)J  ( - — )  ]cosmO}e-iU>  t, 


c  m'  c 
o  o 


(95) 


where  is  the  greater  (lesser)  of  r  and  b  and  Am  are  given  by  Eqs.  (91). 

The  primary  field  to  be  used  in  the  simple-source  term  of  the  second- 
order  nonlinear  differential  Eq.  (48)  is  obtained  by  adding  the  cylindrical ly 
radiated  field  at  frequency  w'  to  Eq.  (95),  obtaining: 


p  *  ah  <l>(iL!)e-l**,t  +  l  [A  H  (1)(^) 

(l)  °  V  «io  nn 


.. .  u>  r.  air.  ... 

+  (2-6  )A'H  'A\— — )J  (— — •)  JcosmGe  1U  \ 
mo  m  c  me 
o  o 


(96) 


where  A  =  pressure  amplitude  of  the  cylinder  oscillating  at  angular  frequency 
a)'  and  6mo  is  the  Kronecker  delta. 

We  obtain  the  real  part  of  P^  using  the  fact  that 
Re ( p ( 1 ) I  3  lP(i)+P(i)*]/2-  Hence, 


Re(P 
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If  we  define  as  the  portions  of  that  contribute  to  the 

sum-  and  difference-frequency  pressures,  we  obtain 


ip  2i  =  V  *  h  (l>f— MH  (i>(— U 
lP(l)  4  4  2  Ho  (  c  )An 
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and 


(1) 


m=0 


2  o 


c  r  ‘  m 
o 


c  ra 
o 


(2)  W  r>  “  r< 

+  (2-6  )A’H  '*'(  )J  (  -  -- -)] 
mo  m  c  m  c  ‘ 
o  o 


-i  (a)’ -w"  )t  , 
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The  simple  source  term  of  Eq.  (48)  may  now  be  obtained  using 


<Po  If’t 
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- 4  ~2  |P(I)  4- 
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This  results  in: 
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(ioi) 


By  Eq.  (80) , 


<(,o  lf>* 


-1 

2 


-iM  t 

[B±e  +  c.c] 


It  should  be  noted  that  m_  *  ju'  -  w" |  and,  therefore,  B_  are  minus  the 
coefficient  of  e  ,  which  is  explicitly  shown  in  Eq.  (101) 

when  «’>«•'.  It  is  the  complex  conjugate  of  this  when  u>”>m' .  For  m’>m", 
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where  k_  -  ((*)'-u»”)/c0- 

The  angular  Integrals  can  be  performed  using  the  fact  that 
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This  yields  (for  r<b): 


(102) 
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(103) 


where  r<b  and  io^w"  . 

For  r>b,  the  integrals  in  Eq.  (103)  run  from  a  to  b,  plus  an  additional 
set  of  integrals  to  those  in  Eq.  (103)  is  required  in  which  the  roles  of  r' 
and  b  are  exchanged  for  the  underlined  terms  and  the  limits  run  from  b  to  r. 
However,  one  additional  difficulty  arises  when  r>b.  This  is  the  fact  that  the 
source  Hankel  function  has  a  singularity  within  the  volume  of  integration. 

This  difficulty  can  be  circumvented  by  excluding  this  singular  region  via  the 
mathematical  artifice  of  enclosing  the  u>"  cylinder  with  a  surface  of  radius  ■; 
and  analyzing  the  contribution  to  the  solution  from  the  associated  surface 
Integral.  It  is  relatively  straightforward  to  demonstrate  that  this  surface 
contribution  is  proportional  to  e£ne  and,  hence,  gives  a  vanishing 
contribution  in  the  limit  as  t  ♦  0. 


The  sum-  and  difference-frequency  pressures  can  once  again  be  specified 

in  dimensionless  form  in  terms  of  the  parameters  T  ■  u^t ,  C  *  u>±r/c q , 

0,  ut.a/c  ,  u>  b/c  ,  as  well  as  the  frequency  ratio, 
x  o  x  o 


At  this  point  we  note  that  estimating  a  region  of  validity  for  Eq.  (103) 
is  similar  to  the  case  of  plane-wave  scattering  by  a  vibrating  cylinder. 
Hence,  we  postpone  this  calculation  until  the  solution  to  this  latter  problem 
is  obtained  [see  discussion  following  Eq.  (115)]. 


2 .  Connection  With  Previous  Research 

In  1966,  Lauvstad  [21]  solved  the  problem  of  two  eccentric  cylindrical 
waves  simultaneously  present  in  a  fluid  medium.  The  present  work  differs  from 
Lauvstad's  in  two  important  ways.  First,  Lauvstad  solved  only  the  radiation 
problem;  i.e.,  no  scattering  of  the  primaries  from  the  cylindrical  surfaces 
was  considered.  Secondly,  Lauvstad  used  the  Green's  function  that  vanishes  at 
the  cylindrical  surfaces  for  the  second-order  solution.  This  corresponds  to 
the  rather  physically  unrealizable  situation  in  which  the  radiating  cylindri¬ 
cal  surface  is  acoustically  soft.  Since  the  current  work  used  the  Green's 
function  whose  normal  derivative  vanishes  at  the  cylindrical  surface  (cor¬ 
responding  to  the  more  realistic  rigid-body  case),  no  limits  can  be  taken  to 
establish  correspondence  between  the  results  of  Lauvstad  and  the  current  work. 

It  is  nonetheless  of  interest  to  obtain  an  expression  for  the  difference- 
frequency  pressure  for  the  radiation  problem  considered  by  Lauvstad  on  the 
basis  of  the  current  theory  with  rigid  boundaries  replacing  Lauvstad's  soft 
boundaries.  All  that  Is  required  is  to  drop  the  terras  in  Eq .  (103)  containing 
the  scattering  coefficients  A^*  (this  is  equivalent  to  inhibiting  the 
scattering  process  in  the  first-order  fields). 

This  results  in: 
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Although  no  general  correspondence  can  be  made  between  the  results 
presented  above  and  those  of  Lauvstad  (due  to  the  different  choices  of  Green’s 
functions),  it  is  possible  to  establish  a  connection  in  the  asymptotic 
limit  r-»®.  This  is  due  to  the  physically  reasonable  result  that  the  effects 
of  the  surface  contributions  in  this  problem  have  become  negligible  at 
distances  far  from  the  surfaces.  It  is  relatively  easy  to  demonstrate  in  this 
limit  that  Lauvstad's  Eq.  (29)  and  Eq .  (103)  above  reduce  to  the  same 
expression  in  this  limit.  It  is  necessary  again  to  discard  from  Eq.  (103)  the 
terms  corresponding  to  surface  scattering.  It  is  also  necessary  to  discard 
the  terms  which  reflect  the  different  boundary  conditions  satisfied  by  the 
Green's  functions  (namely,  the  terms  with  coefficients  involving 
derivatives).  It  is  also  necessary  in  Lauvstad's  expression  to  discard  the 
integrals  with  infinite  limits,  which  has  already  been  done  in  producing 
Eq.  (103).  Finally,  the  following  correspondence  between  constants  in  the  two 
treatments  must  be  made: 

Lauvstad's  Notation  Notation  Used  Here 


A  +  2 


2T 


A!  A'/(P0u>") 

A2  A/(pou)') 


When  these  relations  and  simplifications  are  used,  Eq.  (103)  will  reduce 
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for  t+~  and  r>b. 
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Lauvstad's  Eq.  (29)  reduces  to  V2  times  this  result.  However,  an 
apparent  algebraic  slip  occurred  when  Lauvstad  obtained  his  Eq.  (29)  from  his 
Eq.s  (26,  27,  and  28).  [His  Eq .  (29)  should  have  a  divisor  of  8,  not  16  as 
listed  in  Lauvstad’s  article.] 


G.  Plane-Wave  Scattering  from  a  Vibrating  Cylindrical  Surface 

In  the  present  section,  the  problem  of  the  generation  of  sum-  and 
di f ference-f requency  waves  arising  from  the  scattering  of  a  plane  wave  (of 
angular  frequency  up)  normally  incident  on  a  cylinder  that  deforms  radially 
and  uniformly  (at  angular  frequency  u>c )  will  be  considered. 


1 .  Censor-Method  Solution  -  „ 

The  problem  of  the  scattering  of  a  plane  wave  of  angular  frequency  Up 
normally  incident  on  a  cylinder  vibrating  radially  with  angular  frequency  wc 
(see  Fig.  4)  was  a  problem  considered  by  Censor.  Substitution  of  his  Eqs. 
(24)  Into  his  Eqs.  (10)  and  (6)  results  in  his  solution  for  the  sum-  and 
difference-pressure  waves 
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(104) 


where  q  2  pa/c0,  Pp  Is  the  plane-wave  pressure  amplitude  and  where  Pc  is  the 
cylindrical-wave  pressure  amplitude  defined  by  P  *  PcH0(kt,r),  p  =  kpa  ,  k(,  = 
wave  number  associated  with  the  cylindrical  wave,  a  =  radius  of  cylinder,  and 


“±  3  w 


In  this  expression,  Censor's  small  parameter  e  has  been  replaced 
by[Pc|H1(1)(kca)| I/(Ponaco),  obtained  by  requiring  the  pressure  and 
displacement  at  the  surface  of  the  cylinder  to  be  consistent  with  Eq.  (51). 
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INCIDENT  PLANE  WAVE  OF 
ANGULAR  FREQUENCY  wp 


INFINITE  CYLINDER 
VIBRATING  RADIALLY 
AT  ANGULAR  FREQUENCY  ojc 


Fig.  4  -  Geometry  of  plane  wave  scattering  from 
a  vibrating  cylinder 


2 .  Solution  Using  Second-Order  Nonlinear  Wave  Equation 

We  must  now  solve  the  second-order  nonlinear  wave  Eq.  (48)  using  as  P^ 
the  sum  of:  l)  solution  of  the  linear  wave  equation 


(105) 


for  the  problem  of  linear  rigid-body  scattering  from  a  cylinder  plus  2)  the 
linear  solution  for  radiation  from  a  cylindrical  source.  Thus  ®ay  be 

represented  as 


P,, .  -  P,  +  P  +  P 

(l)  inc  scatt  rad 

with  obvious  meaning  for  the  subscript  notation.  The  expressions  for  the 
well-known  linear-wave  equation  solutions  may  be  found  in  any  standard 


(106) 
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in  nm. 


acoustic  text,  such  as  Morse  and  Ingard  [44].  They  are 
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where  PQ  =  pressure  amplitude  of  incident  plane  wave 


A  =  pressure  amplitude  of  radiated  cylindrical  wave 


u)p  =  angular  frequency  of  Incident  plane  wave 


w  =  angular  frequency  of  cylindrical  radiated  wave. 


=  mt^-order  Bessel  function  of  the  first  kind 
in 


=  m^-order  Hankel  function  of  the  first  kind 
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Nm  *  mt,1-order  Neumann  function 
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a  *  cylinder  radius 


cQ  *  linear  sound  speed. 

The  boundary  condition  appropriate  for  obtaining  the  linear  rigid-body 
scattered  solution  is  that  the  fluid  particle  velocity  vanish  at  the  surface 
of  the  cylinder. 

Once  again,  in  order  to  solve  Eq .  (48),  we  resort  to  the  representations 
of  P(2)  anc*  P0<*q/»0  provided  by  Eqs.  (79)  and  (80).  Also,  as  before,  the 
surface  Integral  of  Eq.  (83)  does  not  contribute  to  the  solution  for  the  same 
reason  given  following  that  equation.  Hence,  the  solution  may  still  be  rep¬ 
resented  by  Eq.  (84).  Since  the  geometry  of  the  current  problem  is  cylindri¬ 
cal,  Eq.  (85)  still  provides  the  appropriate  Green's  function  [contributions 
to  the  sum-  and  difference-frequency  acoustic  pressures  from  regions  beyond 
the  point  of  interest  are  again  neglected,  for  the  same  reasons  given  follow¬ 
ing  Eq.  (85)]. 

As  before,  in  order  to  have  a  representation  of  the  solution,  a  represen¬ 
tation  of  Bn(r')  must  now  be  obtained.  This  is  done  by  substituting  the 
assumed  forms  [Eqs.  (79)  and  (80)]  into  the  second-order  wave  Eq.  (48).  Once 
again,  we  must  carefully  handle  the  complex  quantities  involved.  We  again  use 
the  theorem  that  Re(z)  =  (z+z*)/2  to  help  obtain  a  representation  for  t lie 
first-order  solution  P^^(r).  This  provides 
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We  note  that  P^j(r)  is  of  the  form 
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By  the  same  argument  as  given  for  the  case  of  cylindrical  waves  incident 
on  a  vibrating  cylinder,  the  Bn(r'  )  are  simply  minus  the  coefficient  of  the 
negative  time  exponential  in  Eqs.  (112)  and  (113).  Two  different  solutions 
are  obtained  for  the  difference-frequency  depending  on  whether  (Dp>u>c  or  “c>u>p> 
since  each  of  these  cases  will  give  a  different  coefficient  from  the  second  of 
Eqs.  (112)  and  (113).  Since  the  procedures  are  similar  in  each  of  the  three 
possible  cases,  we  choose  a  particular  one  to  represent  the  solution,  namely 
the  difference-frequency  pressure  when  u>  >0^  (this  is  chosen  since  it 
represents  the  case  for  which  experiments  were  performed). 


We  may  proceed  to  obtain  the  difference-frequency  pressure  P_(r)  by 
substituting  the  expression  obtained  for  Bn(r')  by  the  above-described  method 
and  the  Green's  function  gn(r,r')  of  Eq.  (85)  into  Eq.  (84).  The  integrals  on 
0'  may  again  be  performed  using  the  theorem  preceding  Eq.  (103). 

This  results  in  the  following  expression  for  the  difference-frequency 
pressure: 
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Once'&gain,  the  sum-  and  difference-frequency  pressures  can  be  stated  in 
nondlmensional  te<ms .  In  this  case  the  relevant  parameters  are 
t  =  u>^_t ,  C  “  and  t*1e  fre1uency  ratio. 

At  this  point  we  undertake^a  calculation  of  the  region  of  validity  of 

Eq.  (115).  First  we  consider  the  siftajlest  value  of  r  for  which  Eq.  (115)  may 

\ 

reasonably  be  expected  to  be  valid.  The^fev^re  several  sources  of  error  at 
radii  close  to  the  cylindrical  surface.  These^'sjre:  1)  contributions  due  to 
the  neglected  D'Alembert ian  terms  of  Eq.  (47),  2)  contributions  from  the 
neglected  Green's  function  integral  between  r  and  •»,  and  ^^  contributions  due 
to  the  neglected  Green's  function  surface  Integral.  All  of  these  are 
estimated  in  Section  IV.  B.  From  the  computations  in  that  section,  it  is 
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clear  that  an  upper  bound  on  this  error  is  102  at  r  ■  3  cm  (these  errors 
decrease  approximately  as  1/r  at  distances  beyond  this  radius).  Hence  r  « 

3  cm  can  be  taken  as  a  reasonable  lower  limit  of  validity. 

Next  we  consider  the  largest  value  of  r  for  which  Eq.  (113)  is  valid.  As 
was  done  In  the  case  of  plane-wave  scattering  from  a  vibrating  plane,  we 
effect  this  estimate  by  comparing  the  energy  densities  of  the  secondary  waves 
to  that  of  the  primary  waves.  A  conservative  estimate  can  be  made  by 
including  only  the  cylindrically  radiated  wave  and  the  rigid-body  scattered 
waves  in  the  estimate  of  the  primary  field's  energy  density.  We  choose  to 
represent  each  of  these  waves  respectively  by  the  simple  formulas  Pf  =  Ar/^r 
and  Pg  »  Ag//r.  For  the  experiment  described  in  Chapter  III,  the  empirical 
coefficients  have  the  approximate  values  Ag  »  5x10^  Pa  and  Af  *  1.36*10^  Pa* 
Again,  a  conservative  estimate  of  the  primary  energy  density  can  be  computed 
by  using  the  plane-wave  formula  and  adding  the  results.  This  results  in  an 
approximate  primary  energy  density  of  (4.66x10  m)/r. 

There  are  two  distinct  angular  regions  for  the  secondary  waves:  angles 
near  0°  and  angles  far  from  0°.  At  angles  near  0°,  the  difference-frequency 
pressure  grows  approximately  linearly.  From  the  results  of  Chapter  III,  the 
formula  P_  =  300  Pa/m  r  is  seen  to  be  approximately  followed.  Scaling  for  the 
sum-frequency  case  gives:  P+  =  1300  Pa/m  r.  Once  again,  we  use  the  very 
conservative  estimate  that  the  pressures  associated  with  the  second  harmonics 
follow  the  plane-wave  Fubini-Ghiron  formula  (of  course,  the  waves  will  actual¬ 
ly  grow  much  more  slowly  in  this  case).  Using  as  a  typical  value  the  second- 
harmonic  pressure  formula  obtained  in  the  section  on  plane  wave  and  scattering 
from  a  vibrating  plane,  we  have  P2w  -  372  Pa  r  (we  will  use  this  for  the  har¬ 
monics  of  each  of  the  primaries).  At  a  distance  of  1  m  from  the  cylindrical 
surface,  the  energy  densities  of  the  secondaries  are  less  than  IX  of  the 
primaries.  Hence,  a  reasonable  region  of  validity  of  Eq.  (115)  may  be  taken 
to  be  r  ■  15  cm  to  r  -  100  cm.  (The  discontinuity  distance  calculated  in  the 
section  on  plane-wave  scattering  from  a  vibrating  plane  can  still  be  taken  as 
a  reasonable  estimate.  Since  it  was  approximately  10  m,  viscous  terms  may 
reasonably  be  said  to  play  no  significant  role  in  the  above  estimate.) 

In  Eq.  (115)  the  linear  rigid-body  scattering  coefficients  appear 
explicitly.  In  the  Censor  solution,  given  by  Eq.  (104),  these  do  not  appear 
because  Censor  chose  to  Incorporate  the  expressions  for  the  coefficients  into 
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his  solution.  It  Is  desirable  to  re-express  his  solution  in  a  form  in  which 
these  coefficients  appear  explicitly  as  they  do  in  Eq.  (115).  This 
facilitates  the  calculation  of  the  difference-frequency  pressure  in  the  case 
in  which  the  surface  is  not  rigid  and  the  scattering  coefficients  Ag,  must  be 
empirically  determined.  In  terms  of  A^,  Censor's  solution  becomes: 
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Three  methods  were  developed  to  analyze  the  integrals  in  Eq.  (115):  1) 

Numerical  integration  by  the  use  of  Gauss  quadrature.  2)  A  new  integration 
technique  that  enables  the  calculation  of  the  integrals  in  terms  of  sums. 

3)  Expression  of  the  integrals  in  closed  form  in  terms  of  known,  although 
rarely  encountered  associated  Bessel  functions  for  the  case  kca»l. 

Method  1  proved  to  be  the  most  direct  and  efficient.  Numerical  results 
will  be  given  in  Section  III.  Methods  2  and  3  are  discussed  in  the  Appendix. 

3 .  Connection  with  Previous  Research 

In  1962,  Dean  [20]  solved  the  problem  of  two  concentric  cylindrical  waves 
interacting  nonlinearly.  He  obtained  the  following  solution  for  the  sum- 
frequency  pressure  in  the  farfleld  (for  the  case  where  a+0): 

p.  -  »PP.(2pc2r1(r)(k  k.)/2k  a2e1V,  (117) 

+  aboo  ab  + 

where  ka  .kjj  *  wavenumbers  of  primaries,  k+  ■  sum-frequency  wavenumber,  Pa,Pb  - 
constants  that  measure  the  acoustic  pressure  amplitudes  of  the  cylindrical 
primaries,  and  Dean's  1-T  has  been  replaced  by  T  to  be  consistent  with  the 
notation  used  in  this  work.  The  constants  Pa  and  P^  are  related  to  the 
primary  pressures  through  the  relations 
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p.<'>  ■  f.Ih0<1,(V>/hi<1><V>1 
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The  solution  to  the  above  problem  can  also  be  obtained  using  Eq.  (115). 

In  order  to  do  so,  the  terms  corresponding  to  the  Incident  plane  wave  must  be 
suppressed  [these  are  all  terms  In  Eq.  (115)  that  are  not  multiplied  by  the 
scattering  coefficients  A^].  Furthermore,  the  scattering  coefficients  A^  must 
be  replaced  with  the  Kronecker  delta.  Lastly,  the  following 

identif ications  must  be  made  between  constants  used  in  Dean’s  work  and 
constants  used  in  the  present  work: 

po  ■  p./ai<l)<V> 

A  -  Fb/H1(1)(kba). 

If  the  Hankel  functions  in  Eq.  (115)  are  all  replaced  by  the  first  term 
of  their  asymptotic  expansions,  an  elementary  integral  is  obtained  and  Eq. 
(117)  follows  apart  from  an  unimportant  phase  factor  of  eilT^. 

III.  NUMERICAL  RESULTS 

This  chapter  presents  in  graphical  form  the  results  of  numerical 
calculations  based  on  the  analytical  solutions  of  the  nonlinear  wave  equation 
obtained  in  Chapter  II.  Censor  gave  only  analytical  expressions  for  his 
theory.  In  order  to  compare  Censor’s  results  with  those  of  the  nonlinear 
theory,  his  analytical  expressions  were  evaluated  numerically.  The  results  of 
these  evaluations  are  also  presented  graphically  in  this  chapter. 

The  configuration  selected  for  experimental  investigation  was  the  one  in 
which  a  plane  wave  is  normally  incident  on  a  vibrating  cylindrical  surface. 
Censor's  solution  [23]  to  this  problem  is  given  by  Eqs.  (5),  (10),  and  (27)  in 
his  paper  and  by  Eq.  (104)  of  this  report. 
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In  order  Co  illustrate  Censor's  solution  (as  well  as  the  solution  to  the 
nonlinear  theory),  numerical  values  that  were  experimentally  realizable  were 
used  to  make  example  calculations.  In  these  calculations,  the  plane-wave 
frequency  was  chosen  to  be  162  kHz,  and  the  cylindrical-wave  frequency  was 
chosen  to  be  102  kHz  (giving  a  difference-frequency  of  60  kHz).  The  amplitude 
of  the  incident  pressure  wave  was  selected  to  be  1.0x10^  Pa,  and  the 
cylindrical-wave  amplitude  coefficient  (A)  was  selected  to  be  3.5*10^  Pa. 

The  angular  distribution  of  difference-frequency  pressure  at  15  cm  from 
the  symmetry  axis  of  the  cylinder  obtained  by  numerically  analyzing  Censor's 
expressions  is  shown  in  Fig.  5.  (The  reasons  for  studying  the  difference- 
frequency  case,  as  well  as  the  reasons  for  selecting  the  particular  experi¬ 
mental  parameters  indicated  in  Fig.  5  will  be  discussed  in  Chapter  IV.)  The 
maximum  pressure  at  this  radius  occurs  at  0°  and  is  0.9  Pa  for  the  parameters 
given.  (In  Fig.  5,  as  well  as  all  other  polar  plots,  the  dB  scale  is  measured 
relative  to  the  maximum  pressure  level  at  the  radius  of  interest.  The  maximum 
pressure  represented  in  a  particular  polar  diagram  is  given  in  the  information 
box  associated  with  it  and  is  referred  to  as  "Pmax"*)  As  discussed  in  Chapter 
II,  this  pressure  value  is  of  the  same  order  of  magnitude  as  pseudosound. 

In  the  figure  captions  for  the  difference-frequency  pressure,  the  value 

2 

of  the  quantity  A^  =  (P^P^O/fp^^  )  Is  listed,  since  this  factor  may  be  used 
to  obtain  a  nondlmenslonal  pressure.  Here,  Pc  (the  actual  maximum  cylindrical 
pressure  amplitude)  is  used  instead  of  the  quantity  A  of  Eq.  (115),  since  the 
pressure  represented  by  A  is  present  at  no  point  in  the  fluid.  Similarly,  the 
quantity  ■  <»>_a/co  Is  also  listed  in  the  caption. 

Also  of  interest  in  this  problem  is  the  variation  of  difference-frequency 
pressure  with  respect  to  distance  from  the  cylinder  symmetry  axis  at  fixed 
angles.  Figures  6,  7,  and  8  present  the  results  of  Censor's  theory  at  0,  90, 
and  180°,  respectively.  These  graphs  can  be  interpreted  in  the  following  way: 
Censor's  theory  predicts  the  generation  of  difference-frequency  waves  (as  well 
as  sum-frequency  waves)  due  to  the  presence  of  boundary  conditions  associated 
with  the  time-varying  nature  of  the  cylindrical  surface.  Hence,  both  the  sum- 
and  difference-frequency  waves  predicted  by  his  theory  are  created  solely  at 
the  surface  of  the  scatterer.  As  the  observation  point  is  moved  to  increas¬ 
ingly  greater  distances  from  the  boundary,  these  sum-  and  difference-frequency 
pressure  waves  must  spread  cyl indrically  (in  a  manner  similar  to  the 
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spreading  of  the  first-order  cylindrical  field).  Hence,  it  is  expected  that 
pressures  associated  with  Censor's  theory  will  (in  the  asymptotic  limit) 
decrease  in-versely  as  the  square  root  of  the  radial  distance  from  the 
symmetry  axis.  This  does  indeed  prove  to  be  the  case  for  the  pressures 
represented  in  Figs.  6  through  8. 

The  solution  to  the  nonlinear  wave  equation  for  this  problem  is  given  by 
Eq.  (115).  Unlike  Censor's  theory,  this  equation  involves  complicated  inte¬ 
grals  over  triple  products  of  Bessel  functions.  A  new  technique  of  integra¬ 
tion  is  presented  in  the  Appendix  for  treating  these  integrals  for  the  case  in 

1 

which  the  arguments  of  the  Bessel  functions  corresponding  to  the  radiated  cyl¬ 
indrical  wave  as  well  as  the  arguments  of  the  Bessel  functions  corresponding 
to  either  the  incident  plane  wave  the  difference-frequency  wave  are  suf¬ 
ficiently  large  to  be  replaced  by  their  asymptotic  forms.  Unfortunately, 
these  conditions  were  not  met  for  the  case  that  was  modeled  experimentally; 
hence,  these  Integrals  had  to  be  evaluated  numerically.  The  numerical  proce¬ 
dure  chosen  to  analyze  the  integrals  was  the  method  of  Gaussian  Quadrature 
[32].  A  32-point  quadrature  was  used.  (Suitable  abscissas  and  weighting 
factors  are  given  in  Ref.  45).  To  obtain  good  accuracy  using  the  32-potnt 
Gaussian  Quadrature,  the  radial  interval  to  be  integrated  (1  to  46  cm)  had  to 
be  subdivided  into  ten  equal  sub-intervals.  The  full  32-point  Gaussian  Quad¬ 
rature  sum  was  used  to  obtain  the  integrals  over  each  of  these  partial  inter¬ 
vals.  The  integrals  up  to  the  final  observation  point  were  obtained  by  adding 
together  the  integrals  over  all  subintervals  below  the  observation  point. 

It  is  necessary  at  this  point  to  justify  that  the  subdivision  scheme 
described  above  will  indeed  suffice  to  calculate  the  Integrals  of  interest. 
First,  it  is  essential  to  state  the  precision  of  the  Gaussian  Quadrature 
procedure.  If  m  quadrature  points  are  used,  the  integral  of  a  polynomial  of 
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degree  2m  -  1  is  represented  exactly  by  this  method  [46].  Hence,  it  is 


reasonable  to  expect  that  if  the  integrand  of  interest  contains  not  more  than 
2m-l  zeroes  over  the  interval  of  integration,  Gaussian  Quadrature  will  provide 
a  reliable  numerical  result. 

In  order  to  obtain  a  reasonable  estimate  of  the  number  of  zeroes 
occurring  in  the  integrands  of  interest  (over  the  partial  intervals  described 
above),  it  should  be  noted  that  the  general  behavior  of  the  zeroes  of  the 
Bessel  functions  of  the  first  two  kinds  can  be  surmised  by  a  careful 
examination  of  tables  listing  their  values  [47].  It  is  clear  from  these 
tables  that  the  spacing  between  these  zeroes  decreases  for  increasing 
arguments.  Hence,  since  the  greatest  density  of  zeroes  occurs  for  the 
greatest  arguments,  use  of  the  asymptotic  expansions  of  the  functions  will 
result  in  an  upper  bound  on  the  number  of  zeroes  that  occur  in  any  of  the 
subintervals  of  interest.  In  the  case  of  J^(kr),  the  location  of  the  zeroes 
may  be  approximated  by  calculating  the  zeroes  of  the  cosine  term  that  occurs 
in  the  lowest  order  of  this  asymptotic  expansion.  These  zeroes  will  occur  at 
values  of  r  that  satisfy  the  relationship 

kr  ~  i  *  ■  i  ~  (2n+1>  i> 


where  n  ,1  =  integers. 

Conversely,  the  above  relationship  may  be  used  to  obtain  an  upper  bound 
on  the  number  of  zeroes  occurring  In  a  given  r  Interval  (for  a  given  wave- 
number  k)  by  determining  the  greatest  integer  a  that  satisfies  this  relation¬ 
ship.  The  number  of  zeroes  is  then  approximately  equal  to  n+1  (since  the 
first  zero  occurs  at  n  <*  0).  It  can  also  easily  be  seen  from  the  expression 
above  that  the  greatest  number  of  zeroes  occurs  for  i  *>  0. 


Since  the  range  of  integration  is  taken  from  r  =  I  cm  to  r  =  46  cm,  the 
total  range  of  integration  is  45  cm  in  length.  Since  the  total  interval  is 


subdivided  into  10  subintervals,  the  limits  of  integration  over  which  the 
greatest  values  of  the  arguments  of  JQ(kr)  occur  (and  hence  Interval  over 
which  the  greatest  density  of  zeroes  occurs)  is  from  r  =  41.5  cm  to  r  = 

46  cm.  The  above  approximate  expression  may  now  be  used  to  determine  the 
number  of  zeroes  of  JD(kr)  corresponding  to  each  of  the  wavenumbers  of 
interest  (associated  with  the  frequencies  162,  102,  and  60  kHz)  over  this 
interval.  Rounding  all  fractional  values  obtained  in  this  way  to  the  next 
greatest  integer  (to  consider  the  worst  case),  the  results  of  Table  I  are 
obtained . 


Table  l.  Maximum  Number  of  Zeroes  of  JQ(kr)  on  the 
Interval  r  =  41.5  cm  to  r  =  46  cm 


Frequency  Maximum  No.  of  Zeroes 

162  kHz  11 

102  kHz  8 


60  kHz 
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The  maximum  number  of  zeroes  of  the  integrand  involving  the  product  of 
the  three  JQ(kr)  functions  herein  considered  is  given  by  the  sum  of  the  number 
of  zeroes  for  each  of  the  individual  functions.  This  gives  26  zeroes  in  the 
current  example. 

Since  ra  =*  32  Gaussian  Quadrature  points  are  used,  2m-l  =  63  zeroes  would 
still  result  in  an  accurate  value  for  this  Integral  via  this  numerical 
method.  Since  the  current  example  represents  the  worst  case  (in  the  sense 
that  no  other  integrand  of  interest  will  have  more  than  26  zeroes  over  any  of 


t tie  subintervals  being  considered),  It  may  reasonably  be  expected  that  the 
Gaussian  Quadrate  numerical  integration  scheme  chosen  is  adequate  to  perform 
all  the  required  integrals. 

The  results  of  this  numerical  computation  were  checked  in  the  high- 
frequency  limit  where  the  new  technique  of  integration  was  appropriate. 
Excellent  agreement  was  obtained  between  the  Gaussian  quadrature  results  and 
corresponding  results  using  the  series  given  for  these  integrals  by  Eq.  (A21>. 

In  addition  to  verifying  the  results  of  the  numerical  Integration  by 
comparison  with  numerical  results  obtained  using  Eq.  (A21),  another  verifica¬ 
tion  method  was  also  available.  A  few  of  the  ten  subintervals  were  further 
subdivided  and  then  evaluated  using  the  32-point  quadrature  over  each  of  the 
smaller  subintervals.  The  numerical  results  obtained  by  this  further  sub¬ 
division  were  always  in  good  agreement  with  the  results  obtained  with  the 
original  subdivision  scheme,  thus  showing  that  the  original  subdivision  was 
sufficient  for  evaluating  the  integrals.  (It  should  be  noted  that  this  last 
method  is  applicable  even  for  frequencies  that  are  not  large  enough  to  allow 
the  application  of  Eq.  (A21)  to  the  integrals  of  interest.  This  at  least 
provides  a  check  of  consistency.] 

The  angular  distribution  of  difference-frequency  pressure  obtained  by 
analyzing  Eq.  (115)  numerically  (at  5,  10,  and  15  cm  from  the  scatterer's 
center)  is  presented  in  Figs.  9  through  11,  respectively.  (Figures  12  through 
14  present  radial  plots  of  difference-frequency  pressure  at  0,  90,  and  180°. 
These  will  be  discussed  shortly.)  These  graphs  may  be  interpreted  qualita¬ 
tively  in  the  following  way:  For  angles  near  0°,  strong  contr ibut ions  to  the 
difference-frequency  pressure  are  obtained  both  from  the  "mixing’  of  the 
Incident  plane  wave  with  the  cyl indrical ly  radiated  wave  and  from  the  "mixing" 
of  the  rigid-body  scattered  wave  with  the  cyl Indr  leal ly  radiated  wave.  Since 
the  cyl Indrlcally  radiated  wave  does  not  vary  with  angle,  the  variation  of 
difference-frequency  pressure  with  angle  should  be  related  primarily  to  the 
angular  dependence  of  the  sura  of  the  Incident  plane-wave  pressure  and  the 
rigid-body  scattered  pressure.  Of  course,  since  the  nonllnearly  generated 
field  is  expressed  in  terms  of  an  integral  from  the  cylindrical  surface  to  the 
observation  point,  this  field  is  not  understandable  simply  in  terras  of  the 
primary  fields  that  happen  to  be  located  at  the  observation  point.  More 
illuminating  in  this  respect  is  the  evolution  of  the  primary  field  as  one 
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moves  from  the  cylinder  surface  to  the  observation  point.  Figures  15  through 
21  present  the  total  pressure  field  of  the  linear  rigid-body  scattering 
problem  (incident  plane  wave  plus  rigid-body  scattered  field)  at  2-cm 
intervals.  It  Is  interesting  to  note  that  peaks  evolve  in  this  field  at  0  and 
15°  that  correspond  closely  to  peaks  in  the  nonlinear  scattered  field.  At 
angles  not  equal  (or  close)  to  0°,  the  incident  plane-wave  no  longer 
contributes  strongly  to  the  nonlinear  field  (due  to  its  unfavorable  geometric 
relationship  to  the  cyl indrically  radiated  field). 

Hence,  for  these  angles,  the  angular  dependence  of  the  rigid-body  scat¬ 
tered  field  alone  is  more  appropriate  in  interpreting  the  angular  dependence 
of  the  nonlinear  field.  The  evolution  of  the  rigid-body  linearly  scattered 
field  at  2-cm  increments  is  presented  in  Figs.  22  through  28.  The  most  evi¬ 
dent  aspect  of  the  nonlinear  scattered  field  (between  approximately  60  and 
300°)  is  a  pressure  level  nearly  constant  with  angle.  This  corresponds  well 
with  the  linear  rigid-body  scattering  patterns  with  the  exception  of  minima 
located  at  approximately  75,  135,  225,  and  285°  in  these  patterns.  That  these 
minima  in  the  linear  rigid-body  scattered  field  do  not  contribute  signifi¬ 
cantly  to  the  nonlinear  field  may  be  qualitatively  understood  from  the  fact 
that  these  minima  are  not  very  wide  in  terras  of  angle.  (The  minima  at  75  and 
285°  start  with  ~20  degrees  of  width  and  decrease  to  ~10  degrees  of  width. 

The  minima  at  135  and  225°  remain  at  ~5  degrees  of  width  up  to  the  last  radius 
of  Interest.)  Also,  in  the  case  of  the  minima  at  75  and  285°,  the  Incident 
plane  wave  will  still  contribute  something  to  the  nonlinear  field. 

Even  though  correspondence  between  the  primary  fields  and  the  nonlinear 
field  is  not  exact,  this  Is  not  a  serious  matter  since  the  primary  fields  can 
be  used  only  as  a  very  rough  guide  to  the  behavior  of  the  nonlinear  fields. 

It  must  be  remembered  that  the  nonlinear  field  is  calculated  via  volume 
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integrals  over  the  primary  fields  (not  just  integrals  along  the  radius); 
hence,  detailed  variations  in  angle  in  the  primary  fields  can  easily  disappear 
In  the  nonlinear  field. 

The  angles  between  15  and  60°  and  between  300  and  345°  do  not  correspond 
closely  to  either  of  the  two  sets  of  patterns  given  in  Figs.  15  through  21  or 
Figs.  22  through  28  (although  some  correspondence  can  still  be  seen  between 
certain  features).  These  angles  may  be  regarded  as  a  "transition  region” 
through  which  the  effects  associated  with  mixing  between  the  incident  plane 
wave  and  cylindrlcally  radiated  wave  diminish. 

Behavior  of  the  difference-frequency  pressure  at  fixed  angles  and  varying 
distance  is  illustrated  in  Figs.  12  through  14.  Figure  12  gives  the  differ¬ 
ence-frequency  pressure  at  0°;  Figs.  13  and  14  give  the  difference-frequency 
pressures  at  90  and  180°,  respectively.  It  will  be  noted  in  Fig.  12  that  the 
difference-frequency  pressure  in  this  direction  increases  approximately  lin¬ 
early  with  distance.  This  result  is  similar  to  the  behavior  of  the  parametric 
array  [14]  (in  whLch  a  single  piston  source  Is  driven  at  two  different  primary 
frequencies).  It  Is  understandable  that  the  difference-frequency  pressure  at 
0“  In  the  current  problem  should  behave  approximately  as  a  parametric  array 
since  the  Incident  plane  wave  acts  exactly  as  one  of  the  primary  waves  does  in 
the  parametric  array,  and  the  radiated  cylindrical  wave  approximates  the 
behavior  of  the  second  primary  (although  it  diminishes  in  amplitude). 

At  angles  other  than  0°,  however,  the  geometrical  relationship  between 
the  incident  plane  wave  and  cylindrlcally  radiated  wave  is  no  longer  favorable 
for  the  nonlinear  generation  of  acoustic  waves.  Hence,  at  these  angles,  it  is 
the  mixing  of  the  rigid-body  scattered  wave  with  the  cylindrlcally  radiated 
wave  that  is  responsible  for  the  production  of  the  majority  of  the  nonlinear 
field.  Since  the  rigid-body  scattered  wave  also  spreads  cylindrlcally,  the 
interaction  between  these  two  waves  is  similar  to  the  interaction  of  two 
concentric  cylindrlcally  radiated  waves.  This  last  problem  was  considered  by 
Dean  [20],  who  showed  that  in  this  case,  unlike  the  parametric  array,  the 
nonllnearly  generated  waves  approach  a  constant  value  as  the  observation  point 
approaches  the  farfleld.  A  similar  behavior  is  apparent  in  the  present 
problem  from  Figs.  13  and  14. 

It  is  interesting  at  this  point  to  attempt  to  establish  an  actual  numer¬ 
ical  connection  between  the  asymptotic  expression  given  for  this  case  by  Dean 


(reproduced  in  Section  III.  G.  3  of  this  report)  and  the  experimental  para¬ 
meters  resulting  in  Figs.  13  and  14.  We  must,  of  course,  rewrite  the  expres¬ 
sion  so  it  represents  the  difference-frequency  case,  and  substitute  the  nota¬ 
tion  used  for  the  pressure  amplitudes  used  in  this  thesis.  This  results  in 
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We  are  now  faced  with  the  problem  of  choosing  appropriate  values  for  the 
pressure  amplitudes  Pq  and  A,  which  are  consistent  with  the  present 
experimental  parameters  and  which  realistically  represent  the  theoretical 
situation  of  concentric  cylindrical  waves.  Clearly,  one  of  these  ought  to  be 
chosen  as  3.5*10^  Pa,  the  actual  pressure-amplitude  coefficient  of  the 
cylindrical  source  used  in  the  experiment.  The  choice  of  the  other  pressure 
amplitude,  however,  is  more  subtle. 

The  scattered  pressure  plotted  in  Figs.  13  and  14  has  been  computed  using 
Eq.  (108).  It  is  this  equation  that  Is  used  as  a  guide  in  selecting  the 
second  required  pressure  amplitude  In  the  asymptotic  calculation.  We  are 
representing  our  cylindrically  radiated  wave  by  an  expression  of  the  form 
AHo^^(kr).  In  analogy  with  this  expression,  we  select  the  zeroth-order 
scattering  coefficient  as  the  required  second  pressure  amplitude.  The 
scattered  pressure  plotted  in  Figs.  13  and  14  actually  results  from  a  large 
number  of  terms  in  the  series  represented  by  Eq.  (108).  However,  the  geomet¬ 
rical  collimation  of  the  cylindrically  radiated  wave  is  clearly  strongest  with 
respect  to  the  zeroth-order  scattering  term.  Hence,  we  expect  the  majority  of 
the  difference-frequency  pressure  will  be  generated  via  the  "mixing"  of  the 
cylindrically  radiated  wave  with  the  zeroth-order  partial  wave  of  the  scatter¬ 
ed  pressure  field.  We  compute  this  coefficient  using  the  expressions 
following  Eq.  (109)  with  1.0x10^  Pa  as  the  incident  plane-wave  amplitude 
(corresponding  to  the  actual  experimental  parameter).  This  results  in  a  value 
of  approximately  3*10^  Pa  for  the  zeroth-order  coefficient.  Using  this  as  one 
of  the  required  pressure  amplitudes  (with  3.5x10^  Pa  for  the  other  pressure 
amplitude)  in  the  asymptotic  expression  for  the  difference-frequency  pressure, 
yields  a  value  of  approximately  41  Pa.  In  Figs.  13  and  14  the  difference- 
frequency  pressure  has  obtained  a  value  of  approximately  30  Pa  at  r  *  45  cm 
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and  Is  continuing  to  gradually  Increase.  This  Is  certainly  a  reasonable 
agreement  and,  hence,  gives  further  verification  that  the  computer  routines 
are  providing  accurate  calculations  of  the  difference-frequency  pressures. 

IV.  EXPERIMENT 

A.  Introduction 

An  experimental  Investigation  was  undertaken  to  confirm  the  theoretical 
predictions  for  the  case  in  which  a  plane  wave  is  normally  incident  on  a 
cylindrical  surface  that  deforms  harmonically  and  uniformly  in  the  radial 
direction.  A  discussion  of  the  choice  for  the  experiment  including  the 
principles  upon  which  the  experimental  parameters  were  selected  Is  given  below 
In  Section  IV.  B.  Although  the  investigation  was  unsuccessful  in  confirming 
the  theoretical  predictions,  it  was  nonetheless  successful  in  identifying  the 
several  difficulties  that  arise  in  nearfield  nonlinear  experiments  and 
resolving  all  but  one  of  those  identified.  A  discussion  of  these  difficulties 
and  the  solutions  that  were  achievable  are  given  in  Section  IV.  C. 

Tn  addition  to  these  positive  aspects,  the  experimental  investigation 
also  produced  significant  results  in  several  other  areas,  primarily  with 
regard  to  the  selection  and  calibration  of  the  sound  sources  and  receivers 
used  in  the  experiment.  A  description  of  these  results  is  given  in  Sections 
IV.  D,  E,  and  F. 

First,  the  selection  and  design  of  the  sound  sources  is  discussed  in 
Section  IV.  D.  Second,  the  selection  and  the  first-order  (linear)  calibration 
of  the  sound  sensors  (hydrophones)  is  presented  in  Section  IV.  E.  Lastly, 
Section  IV.  F  describes  the  nonlinear  calibration  of  these  same  (and  a  few 
additional)  hydrophones  both  by  a  previously  developed  technique  [48,49]  and 
by  a  new  method  developed  in  this  work. 

B.  Choice  of  the  Experiment 

Prior  to  describing  the  actual  experiments  performed,  it  is  worth-while 
describing  the  principles  upon  which  the  experimental  parameters  were 
selected.  Initial  decisions  were  required  as  to  which  experimental  geometry 
would  be  addressed,  whether  the  sum-frequency  component  or  the  difference- 


92 


frequency  component  would  be  Investigated  and  what  particular  primary 
frequencies  would  be  suitable  for  measurement. 

In  this  report,  three  geometries  are  considered  theoretically.  These 

are: 


1.  Plane  wave  normally  incident  on  a  uniformly  vibrating  infinite  plane. 

2.  Plane  wave  normally  incident  on  an  infinitely  long  cylinder  vibrating 

uniformly  in  the  radial  direction. 

3.  Cylindrical  wave  normally  incident  on  an  infinitely  long  cylinder 

vibrating  uniformly  in  the  radial  direction. 

Although  Case  1  would  be  the  most  straightforward  to  implement  experi¬ 
mentally,  it  is  very  similar  to  the  standard  case  of  two  infinite  plane  waves 
propagating  together  in  a  fluid  medium,  which  has  been  extensively  studied 
previously  [27,31,50-53]  and,  hence,  is  not  of  the  greatest  interest.  In 
choosing  between  the  final  two  cases.  Case  2  appears  to  be  the  better  one 
based  on  calculations  that  show  a  stgnif icantly  greater  amplitude  of  the  dif¬ 
ference-frequency  component  being  generated  than  in  Case  3.  It  will  be  appre¬ 
ciated  that  similar  cases  involving  spherical  geometry  (which  were  not  treated 
theoretically  here)  would  give  an  even  lower  amplitude  difference-frequency 
pressure  since  in  these  cases  the  energy  is  spreading  into  three  dimensions 
whereas  in  the  cylindrical  case  it  spreads  only  into  two  dimensions. 

In  choosing  between  sum-  and  difference-frequency  components,  the 
difference-frequency  component  was  selected  for  experimental  measurement.  The 
primary  reason  for  this  was  that  it  avoids  the  difficulty  of  separating  the 
sum-frequency  component  from  harmonics  of  the  primaries  in  the  hydrophone 
received  signal. 

Having  chosen  the  geometry  of  Case  2,  appropriate  frequencies  must  be 
selected.  The  cylinder  Is  the  source  that  provides  the  experimental  limits  in 
this  regard.  Since  the  cylinder  must  act  as  a  high-amplitude  sound  source 
with  vibration  in  the  radial  direction  that  is  uniform  along  its  axis,  it 
should  be  operated  near  its  lowest  radial  mode  resonance,  l.e.,  breathing 
mode.  In  addition  to  being  nonuniform,  higher  modes  of  vibration 


significantly  limit  this  amplitude  [54].  An  additional  requirement  to  obtain 
uniformity  of  vibration  along  the  cylinder's  axis  is  that  a  segmented,  rather 
than  a  single  piece,  cylinder  be  used  to  avoid  excitation  of  longitudinal 
(length)  modes  of  vibration.  Fabrication  of  such  a  cylinder  is  quite 
complicated.  Fortunately,  one  was  already  available  that  had  been  used  as  a 
Navy  standard.  It  has  a  breathing  mode  resonance  frequency  of  102  kHz.  A 
cylinder  with  a  resonance  frequency  much  greater  than  this  would  be  difficult 
to  make  and  use  since  it  would  have  an  unacceptably  small  radius. 

Assuming  the  cylinder  to  be  operated  at  about  100  kHz,  a  suitable  plane- 
wave  frequency  (and  hence  difference  frequency)  must  be  selected.  There  are 
two  possibilities: 

1.  The  plane-wave  frequency  is  less  than  the  cylindrical-wave  frequency. 

2.  The  plane-wave  frequency  is  greater  than  the  cyl indrical-wave 

f  requency . 

Case  l  is  preferable  since  the  cylinder  is  more  likely  to  behave  as  a 
rigid  body  scatterer  of  the  plane  wave  in  this  case  [3].  Unfortunately, 
certain  practical  considerations  eliminate  this  case  as  a  possibility.  First, 
in  order  to  maximize  the  ability  to  discriminate  between  the  primary  waves  and 
the  difference-frequency  waves  in  the  hydrophone  received  signal,  the 
difference  frequency  must  be  less  than  either  of  the  primary  frequencies. 

(This  allows  all  filters  that  the  hydrophone's  electrical  output  passes 
through  to  be  operated  in  a  low-pass  mode  thereby  providing  maximum 
discrimination.)  Hence,  the  lowest  possible  plane-wave  frequency  is  50  kHz. 
This  frequency  is  unacceptable,  however,  since  It  results  in  a  difference  1 re¬ 
quency  that  is  also  50  kHz.  To  avoid  this,  the  plane-wave  frequency  must  be 
chosen  closer  to  the  cylinder  frequency.  Unfortunately,  the  difference-fre¬ 
quency  pressure  generated  nonllnearly  in  the  water  varies  directly  with  the 
difference  frequency.  Hence,  the  closer  the  two  frequencies  become,  the 
smaller  the  difference-frequency  pressure  becomes.  The  optimum  frequency, 
which  minimizes  the  "closeness"  of  the  difference  frequency  to  either  primary 
frequency  as  well  as  giving  a  substantial  value  for  the  difference  frequency 
is  approximately  75  kHz.  Unfortunately,  the  25-kHz  difference  frequency 
associated  with  100-  and  75-kHz  primaries  is  still  quite  small.  Hence,  a 


plane-wave  frequency  less  than  the  cylinder  frequency  appears  to  be  a  poor 
choice. 


Selecting  a  plane-wave  frequency  greater  than  the  cylinder  frequency 
precludes  the  cylinder  from  behaving  as  a  rigid  body  In  scattering  the  plane 
wave  [3].  However,  it  Is  possible  to  account  for  this  by  modeling  the 
scattered  field  empirically.  This  is  done  by  measuring  the  scattered  pressure 
(both  amplitude  and  phase)  rather  than  assuming  rigid-body  scattering,  and 
then  calculating  how  the  coefficients  in  the  Hankel-function  expansion  of  the 
rigid-body  field  must  be  modified  in  order  to  obtain  the  measured  pressures. 
(No  modification  of  the  Green's  function  used  in  calculating  the  difference- 
frequency  component  is  required  since  the  cylinder  can  still  properly  be 
assumed  rigid  if  the  difference  frequency  is  chosen  to  be  below  the  cylinder 
resonance.)  It  remains  now  to  select  an  appropriate  plane-wave  frequency 
above  the  cylinder  frequency.  Using  the  same  reasoning  as  above,  this  results 
in  a  plane-wave  frequency  of  about  150  kHz  (in  the  experiment  actually 
performed,  the  plane  frequency  used  was  162  kHz). 

At  this  point  it  should  be  recalled  that  in  the  discussion  in  Chapter  II 
pertaining  to  Eq.  (47),  it  was  remarked  that  if  the  D'Alembertian  terms  on  the 
right-hand  side  of  this  equation  become  significant  relative  to  P2 ,  that  the 
solution  obtained  by  their  neglect  is  questionable.  It  is  essential, 
therefore,  to  estimate  the  value  of  these  terms  in  the  current  measurement. 

In  order  to  facilitate  this  estimate,  the  actual  fields  in  the  current 
problem  will  be  replaced  by  planar  fields  of  the  same  frequencies  as  those  of 
interest  (this  then  calculates  the  worst  case  possible).  Hence,  the  first- 
order  fields  will  be  represented  as: 


P0e,(y-V>  +  P,.1(V"V)  +  A.1<kcx-V>. 


In  this  equation  Pq,  A,  kc ,  wc ,  kp,  u>p  have  the  same  meaning  as  in  Chapter 
II*  The  symbol  Pg  stands  for  the  maximum  amplitude  of  the  scattered  pressure 
field.  Equation  (46)  can  be  used  to  relate  the  velocity  potential  to  the 
pressure.  This  relationship  may  be  taken  to  be: 


P  -  -iu>po<)> 
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where  sinusoidal  time  dependence  has  been  assumed.  A  complex  quantity  Z^, 
which  is  related  to  the  first-order  velocity  potential  therefore,  may  be 
defined  as: 
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The  quantity  4>^  is  then  the  real  part  of  Z^ ,  or 

^  =  (Zx+Z t*)/2. 

Also  of  interest  in  analyzing  the  D’Alembertian  terms  of  Eq.  (47)  are  the 
gradient  and  the  time  derivative  of  Complex  quantities  Z^  and  Z3  may  be 

defined,  then,  as 
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{Note  that  in  the  definition  of  Z^  the  vector  nature  of  the  gradient  has  been 
suppressed  since  it  will  have  no  effect.  This  is  due  to  the  fact  that  plane 
waves  are  being  considered,  and  It  is  actually  Z^  that  is  of  interest  in 
analyzing  the  D’Alembertian  terms  of  Eq.  (47).]  The  quantities  an<* 

3$(l)/3t  may  now  be  calculated  from  and  Z3  as 


96 


l^a)]2  -  (ReZ2)2  -  (Z2+Z2*)2/4 


3$  2 

<-g^)  -  (ReZ3)2  -  (Z3+Z3*)2/4. 

The  quantities  of  Interest  In  calculating  the  D* Alembertian  terms  need 
only  be  calculated  insofar  as  the  difference  frequency  Is  affected.  These 
contributions  to  the  difference  frequency  can  be  denoted  by  the  subscript  [-]  , 
and  can  be  related  to  Z^  ,  Z^,  and  Z^  In  the  following  way: 

a2*(n2  1  a2 

[ — i  ±  (Z  z  *)_ 

at  1  ar  11 

ll?*(l)|2|-  +T<Z2Z2*)- 

^^(l)  2  1 

*<-TT>  I-  “  I  (Z3Z3*>-‘ 

Using  the  given  equations  for  Z^,  Z^  ,  and  Zj  these  expressions  may  be  analyzed 
retaining  only  contributions  to  the  difference  frequency.  Furthermore,  since 
calculation  of  the  worst  case  Is  of  interest,  the  exponentials  will  be  dropped 
and  only  the  amplitudes  will  be  retained.  This  procedure  gives: 
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Now,  these  expressions  may  be  combined  with  the  fact  that 
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to  analyze  the  D' Alembertian  terms  of  Eq.  (47).  This  results  in 


3ij>  2  2 

r  1  -2,  *(1)  *  1  i*.  1 2  1  -2  3  .  2 . 

f  2  PoCo  (  3t  J  "2  pol^(l)l  "  2  poCo  3t2  ^(l)  J 


=  -A(PnP  )(u>  ~m  )2/( 2p  c  2a)  u)  ), 
0  s  c  p  00  c  p 


In  order  to  use  this  expression  to  estimate  the  D' Alembertian  terras,  the 
experimental  values  of  and  o>c  may  be  used  directly.  However,  it  would  not 
be  reasonable  to  directly  substitute  the  value  for  A,  since  this  pressure 
occurs  nowhere  in  the  fluid  surrounding  the  scatterer.  It  is  more  reasonable 
to  use  the  value  of  the  amplitude  of  the  pressure  of  the  cylindrically 
radiated  wave  analyzed  at  the  cylinder's  surface.  This  is  approximately 
l. 31x10^  Pa.  For  the  expression  (Pq+Ps)  it  Is  reasonable  to  use  the  maximum 
value  of  the  total  pressure  field  associated  with  the  first-order  rigid-body 
scattering  problem.  The  value  of  this  quantity  is  approximately  1.7*10'’  Pa. 
Using  these  values,  the  above  expression  for  the  terms  under  the  D' Alembert lan 
operator  gives  a  magnitude  of  l.l  Pa*.  This  is  less  than  5%  of  the  values 
predicted  for  the  difference-frequency  pressure  in  all  directions  at  distances 


*It  is  interesting  to  compare  this  value  with  that  obtained  for  pseudosound  in 
Section  IV.  F.  (0.1  Pa).  We  note  that  the  D'Alembertian  terms  are  of  the  same 
general  magnitude  as  pseudosound. 
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greater  than  approximately  10  cm  from  the  cylinder's  center  (compare  with 
Figs.  12-14).  Hence,  these  terms  are  sufficiently  small  so  as  to  not  severely 
affect  the  theoretical  treatment  at  all  distances  greater  than  a  few 
centimeters  from  the  surface  of  the  cylinder. 

Of  course,  at  shorter  distances,  this  value  becomes  a  more  significant 
fraction  of  the  second-order  pressure.  Hence,  the  solutions  represented  by 
Eq.  (15)  are  not  likely  to  be  accurate  at  distances  where  the  predicted 
pressure  amplitude  is  small  (say  of  the  order  of  pseudosound). 

One  more  consideration  is  of  interest  in  demonstrating  the  fact  that  the 

influence  of  the  D' Alembertian  terms  on  the  solution  is  quite  small.  If,  in 

solving  Eq.  (47),  the  D' Alembert ian  terms  are  not  grouped  with  the  second- 

order  pressure,  they  must  be  included  as  a  part  of  the  virtual  source  term  in 

the  Green's  function  solution  represented  by  Eq.  (84).  In  this  case,  the 

9  9 

virtual  source  term  will  include  quantities  of  the  form  □  Z  ,  which  vanishes 
identically  for  Z  =  plane  wave.  Hence,  any  nonvanishing  contributions  must 
arise  from  the  nonplanar  geometry  of  the  actual  primaries,  and  such  effects 
may  reasonably  be  expected  to  be  small. 

Another  consideration  of  some  concern  is  the  effect  that  the 
D'Alembert ian  terms  might  have  on  the  neglected  surface  integral  of  Eq . 

(83).  This  is  very  difficult  to  estimate  accurately,  owing  to  the  complicated 
form  of  the  Green's  function  and  the  first-order  fields  involved.  An 
extremely  simple  estimate  that  may  give  some  indication  of  the  order  of 
magnitude  of  the  error  may  be  obtained  by  inserting  for  the  D'Alembert ian 
terms  in  this  surface  integral  the  constant  value  of  1.1  Pa  obtained  above, 
multiplied  by  the  difference-frequency  wavenumber,  k^  (since  it  is  the 
gradient  of  the  source  terms  which  appears  in  this  surface  integral).  We 
replace  the  Green's  function  by  its  asymptotic  form,  but  simply  analyze  Its 
value  on  the  cylindrical  surface  [see  Eq.  (7.3.17)  of  Reference  29].  Lastly, 
the  surface  integral  itself  reduces  to  simply  2tt  times  the  cylindrical 
radius.  Hence,  the  D' Alembertian  terras  result  in 

(2na)(kd)(l.l  Pa)[  ^  /T7iTk"[a]  »  2.2  Pa 

of  influence  on  the  surface  integral.  This  too  may  be  regarded  as 
sufficiently  small  to  neglect,  provided  the  second-order  pressure  solution  is 
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tW. 


not  expected  to  be  accurate  at  distances  extremely  close  to  the  cylindrical 
surface  (where,  once  again,  "close"  refers  to  distances  at  which  the  second- 
order  pressure  approaches  the  value  of  pseudosound). 


At  this  point  we  can  also  put  the  discussion  of  the  neglect  of  the 
Green’s  function  integral  between  r  and  “  discussed  in  Chapter  II  on  a 
somewhat  firmer  mathematical  ground  for  this  particular  case  by  the  following 
semi-quantitative  argument:  The  integral  of  concern  represents  an  integration 
over  a  product  of  the  virtual  source  term  and  the  Green's  function.  We  may 
symbolically  represent  the  major  contributions  to  this  integral  by  the 
following  two  associated  integrals 


I 


1 


k.2r  2* 

(— .)  /  d0. 


k  2r  2n  00 

l2  =  (~*~2)  /  d0’  /  dr'r'P  (r’)P  (r')G(r,r'), 

P  c  0  r 

o  o 


where  the  coefficient  arises  from  the  formulation  of  the  simple  source,  Pp 
represents  the  incident  plane  wave,  PQ  represents  the  radiated  cylindrical 
wave,  Pg  represents  the  rigid-body  scattered  waves,  and  G(r,r')  represents  the 
Green's  function.  We  consider  1^  first. 


>  i  1c  r ’cosO*  + 

We  represent  Pp(r')  as  P^e  p  and  Pc(r')  (neglecting  an 

_ _ r ' 

unimportant  phase  factor)  as  A  /2/nk  r*  e  c  (we  have  used  the  asymptotic 


form  of  the  Hankel  function  since  k^c  ui  4.18  cm  *r,  for  which  the  asymptotic 
form  is  reasonably  accurate).  For  the  Green's  function  we  use 


G(r, r’)  =  i  HQ(l)(kdw), 

2  2  2 

where  w  ■  (x-x')  +  (y-y' )  .  This  is  the  infinite-space  Green's  function. 

We  will  discuss  the  errors  in  using  this  instead  of  the  proper  Green's 
function  following  the  present  estimates.  This  error  will  prove  to  be 
insubstant lal . 
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Hence,  we  will 


It  is  clear  that  G(r,r')  has  a  singularity  at  r  =  r*. 

r+e  » 

break,  the  Integral  1^  Into  two  parts:  I.  =  /  +  /  ,  where  e  is 

r  r+e 

infinitesimal . 

We  now  note  that  in  I ^ ,  both  an  integral  O'  and  r'  are  required.  We 
consider  the  angular  Integral  first.  It  should  be  noted  that  only  the 
functions  Pp  and  G  have  a  O'  dependence.  Hence,  we  consider  the  associated 
integral 

T  i  f2"  ,  ik  r'cosQ’  (1),.  x 

I0.  “  4  J  dQ  e  P  Ho  <kdw)‘ 

There  is  a  singularity  in  this  integral  at  O’  =0  (since  w  vanishes  at 

this  value  of  0').  Hence,  we  choose  to  resolve  Iq»  as  follows: 

2n-6  6 

Ifl,  =  /  +  /  ,  where  6  Is  an  infinitesimal  angle.  In  the  integral  between 

0  6  -6 

-6  and  6 ,  it  is  reasonable  to  replace  the  factor  etkpr  cosQ  by  eikpr  ,  since 
0'  has  infinitesimal  values  over  the  entire  range  of  integration.  We  also 
note  that  H  ^  ^  ^ (k^w)  may  be  resolved  into  functions  of  r  and  r'  alone  (see 
Ref.  43,  p.  979,  Eq .  8.531.2).  Hence,  this  Integral  becomes 

6  6 

/  ■  i  « V  /  d0'iVkdr)H0<1)(kdr'> 

-6  -0 


+  2  l  J  (k  r)H  ^(k  r’)cos(mO’)] 
m  d  ra 


The  integrals  are  elementary  and  give 


/  -  i  eV’[Vkdr>Ho<1)(kdr,)(26) 


+  4  l  j  (k  r)H  (1)(k  r*)sin(k  5)1 
m-l  m  d  m  d  d 


Both  terms  in  the  brackets  vanish  in  the  limit  5+0.  Therefore, 

6 

lira  /  ■  0.  This  demonstrates  the  singularity  at  0*  ■  0  in  Iq»  is 

6+0  -6 
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Inconsequent lal .  We  consider  next  the  Integral  /  .  Since  the  singularity 

6 

has  been  removed,  no  error  should  result  (for  purpose  of  the  current  estimate) 

if  the  function  H  ^^(k^w)  is  treated  casually.  We  note  that  w  may  be 

rewritten  as  w  =  (r^  +  r’^  -  2rr'cos0’)^.  This  quantity  varies  from  0  to  2r' 
as  O'  varies  from  0  to  2n.  We  choose  to  replace  w  by  r'  (which  is,  in  a 
certain  sense,  an  average  value).  Thus,  Iqi  becomes  simply 


T  l  u  ( l)/i  ,  x  r  *  jn)  ilt  t1 

I0<  =7Ho  <V  >  [  d°  6  P 


cos©' 


We  note  that  this  integral  can  be  performed  in  the  limit  6*0.  This 
result  is  (see  Ref.  43,  p.  482,  Eq.  3.915.2) 


/  dO'e^*<'pr ,C08G'  *  m[J  (k  r’)  +  J  (-k  r')]. 


o  p 


We  choose  to  replace  the  Bessel  functions  by  their  asymptotic  limits, 
gives 


!qi  -  j  Ho<l>(kdr*)  /iff/kdr'  [cos(kdr'  -  ^-)  +  icos^r’  *-  £-)  ] . 


This  can  now  be  used  to  give  the  following  expression  for  the  r '  integral  in 
the  region  between  r*  »  r  and  r'  *  r+e 


rfe  k  TAP.  r+e  . 

/  -J - J  dr’r'H(1)(kr)Hn(1)(k.r-) 


j - 2“  U  > 

r  poco 


x  k7r  -  j  )  +  icos(kdr’  +  £)] 

d 


Each  cosine  terra  can  be  resolved  into  exponentials.  If  we  also  replace 

the  Hankel  functions  with  asymptotic  forms  discussed  above,  the  Integral 
r*e  r-e  eikr* 

/  will  become  a  sum  of  integrals  of  the  general  form  /  dr'  -  , 
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where  k  represents  various  combinations  of  the  relevant  wavenumbers.  We  can 
relate  each  of  these  integrals  to  a  Fresnel  integral  (see  Ref.  32,  p.  300), 


„  ,  .  1  c  cost  .. 

C.  (x)  =  — ■—  J  — 3-  dt  , 

/  2  jt  0  / 1 


r-e  0  r=e 

since  /  =  /  +  /  .  Using  the  asymptotic  expression  for  C^(x)  (Ref.  32, 

r  r  0 

^co  s  t  1 

pp  300  and  302)  we  can  easily  show  that  /  — —  dt  III  —  +  — •  Hence,  the 


0  / 1 


1  11 
integrals  /  are  proportional  to  the  quantity  — —  -  —  ,  which  clearly 

r+c  r 
r 

vanishes  in  the  limit  e-K). 

In  order  to  conclude  the  analysis  of  1^,  we  next  consider  the  r'  integral 

00 

/  .  It  is  clear  that  the  same  analysis  used  on  the  O'  integral  during  the 

r+e  r+c 

discussion  of  the  integral  /  is  still  appropriate.  In  this  case,  however, 

r 

the  rosult  of  integration  on  r'  will  not  vanish,  and  hence  we  must  treat  the 

sum  cos(k^r'  -  ^-)  +  lcos(k^r'  +  ^-)  somewhat  more  carefully.  We  note  that  the 

absolute  maximum  value  that  the  magnitude  of  this  sum  can  have  is  /2  .  In 

00 

order  to  estimate  the  r'  integral  /  we  replace  this  sum  by  this  numerical 

r+e 

quantity  (this  will  furnish  a  conservative  estimate  of  the  upper  bound  for 
this  Integral).  Again  replacing  the  Hankel  functions  by  their  asymptotic 
expressions  results  in: 


oo  k  TAP  00  ik  r 

rj  =  — - -  r  /27V iT  /277k7  /2T7kT  /2  j  dr*  — — 2- 

2  4  c  d  d  /— r 

r+e  p  c  r=e  /r’ 

o  o 


The  Integral  on  r '  can  now  be  analyzed  by  again  using  the  asymptotic  form 
of  a  Fresnel  Integral.  If  we  now  let  e+0,  the  final  result  for  1^  Is 
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P  c  /SIT  k 
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Hence,  JljJ  falls  off  even  more  quickly  than  would  be  expected  on  the  basis  of 
cylindrical  spreading.  This  is  consistent  with  the  argument  of  the 
negligibility  of  the  Green's  function  integral  between  r  and  00  presented  in 
Chapter  II,  which  was  based  on  phase  cancellation  arguments.  Using  the 
present  experimental  values,  we  get  jl^j  3.37/r,  where  the  answer  will  be  in 
Pa  If  r  is  in  cm.  Hence,  at  a  propagation  distance  of  just  10  cm,  the 
contribution  of  1^  is  0.337  Pa.  This  is  very  small  when  compared  with  the 
numerical  values  presented  in  Chapter  III,  which  were  approximately  38  Pa  at  0 
*  0°  and  25  Pa  at  0  =  90°  at  this  radius. 

We  next  consider  12*  The  easiest  way  to  handle  Pr(r’)  in  this  integrand 
is  to  determine  empirical  coefficients  a  and  B  for  the  series  representation 
of  the  field  up  to  the  dipole  term.  This  representation  is 

Ik  r 

P  (r ' )  *  <q  +  P 

p  /iTf 

P 


Using  the  numerical  values  presented  in  Chapter  III  for  the  linearly  scattered 

5  s 

field  at  r  =  15  cm,  these  constants  become:  a  =  3.4x10  Pa,  3  =  1-56x10  Pa 

(where  values  at  0=0°  and  0  =  90°  have  been  used  to  analyze  the 

constants).  The  singularity  at  0'  =  0  in  I2  can  be  shown  to  be 

inconsequential  via  the  same  technique  as  used  for  Hence,  we  can  proceed 

to  analyze  the  angular  integral  without  regard  for  this  singularity.  As 

before,  we  resolve  H  ^^(k.w)  into  a  sura  over  Bessel  functions  of  the 

o  d 

Individual  radii.  Thus,  the  angular  Integral  tnvolved  in  I2  becomes: 

/  d0'  -  i  /  dO'  elkpr,{J  (k  r)H  (l)(k  r') 

*  o  yirvr  o  d  o  d 


+  2  }  J  (k,r)H  V i;(k  ,r' JcosmO’ ] . 

**  m  d  ra  d 
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The  angular  Integrals  are  again  simple  and  result  In 


/  d0' 


ik  r  * 

—  — P  laJo(kdr)Ho(1)(kdr’)  +  BJ^kjOH,  (1)(k,r '  )  ] 


2  A T71 
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If  we  replace  the  Bessel  functions  In  this  expression  by  their  asymptotic 
forms  (as  well  as  the  remaining  HankeL  functions  In  I2)*  and  choose  to  neglect 
phase  factors  (which  can  only  reduce  the  value  of  the  Integral  of  Interest), 
we  obtain  the  further  estimate  for  I2 


I 


2 


Co+fl)elkdr2k.rtA  * 
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We  again  compute  the  r*  Integral  with  the  aid  of  the  asymptotic  form  of  a 
Fresnel  Integral,  giving 


(a+0)2k  1*A 
a 

/k"p  c  2k  2 
coo  p 


Again  we  note  the  falloff  Is  greater  than  that  based  on  cylindrical 
spreading,  due  once  again  to  phase  cancellations.  Using  the  present 
experimental  values,  we  obtain  ■  1.68/r^2,  where  again  the  answer 

is  In  Pa  If  r  Is  In  cm.  This  gives  just  5.3*10"2  pa  at  a  propagation  distance 
of  just  10  cm.  Hence,  we  once  again  note  the  negligibility  of  these  integrals 
relative  to  the  numerical  values  presented  in  Chapter  III,  thus  reinforcing 
their  neglect  during  the  computation. 

We  consider  lastly  the  consequences  of  using  the  infinite  space  Green's 
function  Instead  of  the  correct  rigid-body  Green's  function  for  the  purposes 
of  this  estimate.  We  note  that  the  effect  (In  the  case  of  a  rigid  plane)  of 
the  boundary  terra  Is  essentially  a  doubling  of  the  pressure  associated  with 
the  free-space  Green's  function.  (This  can  be  understood  from  the  fact  that  a 
simple  Image  source  behind  the  plane  Is  all  that  is  required  to  satisfy  the 
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boundary  condition.)  This  case  may  reasonably  be  regarded  as  an  upper  limit 
to  the  potential  scattered  pressure  level  in  any  other  scattering  geometry. 
Since  both  contributions  discussed  above  are  well  under  1  Pa  at  r  =  10  cm,  it 
is  resonable  to  assume  that  even  Including  the  surface  terms  in  the  Green's 
function  would  no  more  than  double  the  result.  Since  this  still  gives  less 
than  1  Pa,  no  great  harm  is  done  via  use  of  the  free-space  Green's  function  in 
making  the  above  estimates. 

Even  with  the  choices  made  following  the  careful  procedure  detailed 
above,  several  difficulties  arise  in  a  nearfield  nonlinear  scattering 
experiment  that  have  not  been  encountered  in  previously  published 
experiments.  These  difficulties  will  be  described  next. 

C.  Difficulties  in  the  Present  Experiment  not  Encountered  in  Previous 
Research 

In  investigating  the  nonlinear  scattering  of  acoustic  waves  by  vibrating 
obstacles,  several  fundamental  experimental  difficulties  arise.  Previous 
experiments  involving  two  high-amplitude  primaries  interacting  (i.e.,  mixing) 
nonlinearity  in  a  fluid  medium  could  be  designed  to  avoid  these 
difficulties.  The  present  experiment  could  not  be  designed  to  avoid  them. 

The  most  important  of  these  difficulties  are: 

l .  Inadvertent  direct  radiation  of  the  sources  at  the  difference  frequency. 

Since  the  sources  are  finite  in  extent,  measurements  must  be  made  in  the 
extreme  nearfield*  of  the  sources  in  order  to  approximate  an  infinite  plane 
wave  and  infinitely  long  cylinder.  Hence,  direct  radiation  of  the  sources  at 
the  difference  frequency  will  tend  to  be  a  greater  source  of  error  when  the 
measuring  hydrophone  is  near  the  sources. 


*In  acoustics  the  term  farfield  (nearfleld)  is  a  relative  term  that  describes 
observation  points  at  distances  large  (small)  compared  with  the  dimensions  of 
the  source  and  the  wavelengths  involved.  In  the  case  of  a  piston  source,  the 
farfield  (nearfleld)  is  defined  to  be  distances  greater  (lesser)  than  the 
distance  to  the  last  maximum  in  the  on-axis  diffraction  pattern.  This  last 
maximum  is  located  approximately  at  a  distance  equal  to  the  square  of  the 
piston  radius  divided  by  the  wavelength.  The  farfield  (nearfield)  is 
frequently  referred  to  as  the  region  in  which  Fraunhofer  (Fresnel)  diffraction 
effects  occur. 
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2 .  Electrical  filtering  problems  due  to  experimental  constraints. 

The  difference  frequency  was  only  about  one  half  the  lowest  primary 
frequency.  In  addition,  the  pulse  lengths  had  to  be  less  than  about  10  cycles 
at  the  difference  frequency  to  avoid  Interfering  reflections  from  neighboring 
surfaces.  Hence,  the  usual  passive  methods  employed  for  electrical  filtering 
In  previous  farfield,  nonlinear  measurements  were  inappropriate. 

3 .  Difference-frequency  voltage  generated  nonllnearly  in  the  hydrophone's 
sensitive  element. 

This  effect,  due  to  nonlinear  mixing  of  the  primaries  in  the  hydrophone, 
provided  larger  difference-frequency  voltages  than  those  produced  by  the 
difference-frequency  pressure  generated  by  nonlinearities  of  the  fluid 
medium.  The  effect  was  observed  for  a  wide  range  of  available  hydrophones. 

Previous  experimental  Investigations  [27,31,50-53)  have  not  had  to 
contend  with  the  limitations  enumerated  above  because: 

(a)  Circular  piston  sources  were  used  in  these  investigations  enabling 
virtually  any  choice  of  frequency.  Then  the  difference  frequency  was 
chosen  very  much  lower  than  the  average  primary  frequency  (high  downshift 
ratio) . 

(b)  Previously  measurements  were  carried  out  in  the  acoustic  farfield  of 
the  sources  (or  at  least  not  in  the  extreme  nearfield). 

Recently,  measurements  were  made  in  the  nearfield  of  circular  and 
rectangular  piston*  sources  (55,56).  In  one  case  [55],  the  average  frequency 
of  the  primaries  was  so  high  (1.435  MHz)  that  the  primaries  were  absorption 
limited.  This  means  that  the  primaries  were  strongly  attenuated  (mostly  by 
viscous  absorption)  by  the  water  through  which  they  propagated  prior  to 
reaching  the  hydrophone  (the  acoustic  absorption  coefficient  varies 
approximately  as  the  square  of  the  frequency  (57)).  In  the  second  case  (56), 


*In  acoustics  the  term  ’’piston"  refers  to  a  planar  surface,  all  points  of 
which  are  moving  at  the  same  velocity.  The  cross-sect  ion  of  a  piston  is 
generally  circular  or  rectangular. 
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although  the  hydrophone  was  not  in  the  farfield,  it  was  not  in  the  extreme 
nearfield  either.  (In  the  current  sense,  the  "extreme  nearfield"  is  defined 
in  the  case  of  a  piston  source  as  being  within  the  Rayleigh  length  of  the 
piston  face.  The  Rayleigh  length  is  the  square  of  the  radius  of  the  piston 
face  divided  by  the  wavelength.)  The  closest  measurement  made  in  this  case 
was  at  a  hydrophone  position  2.72  times  the  Rayleigh  length  from  the  piston 
face.  Hence,  neither  of  these  studies  had  to  face  the  difficulties  inherent 
in  a  nonlinear  scattering  experiment. 

Fortunately,  the  first  difficulty  listed  above  (direct  radiation  at  the 
difference  frequency)  was  fairly  straightforward  to  eliminate.  By  adding  an 
appropriate  pulse-shaping  network  to  the  electromagnetic  driving  pulses  (prior 
to  amplification),  95%  of  the  directly  radiated  difference-frequency  component 
was  eliminated.  Prior  to  the  addition  of  this  network,  the  directly  radiated 
difference-frequency  component  was  comparable  in  magnitude  to  that  generated 
nonlinearly  by  the  fluid  medium.  Thus,  after  the  addition  of  the  network,  the 
direct  radiation  became  a  small  component  of  the  total  difference-frequency 
pressure  field. 

The  second  difficulty  noted  above  (the  problem  of  electrical  filtering) 
was  also  resolved  by  the  selection  of  extremely  linear  state-of-the-art  active 
electrical  filters.  There  are  two  reasons  why  previous  nonlinear  measurements 
were  not  faced  with  this  problem.  First,  earlier  measurements  were  made  at 
least  moderately  far  from  the  sources  and  did  not  involve  the  possibility  of 
unwanted  single  or  multiple  reflections.  Hence,  long  pulse  lengths  could  be 
used  without  interfering  reflections  being  received  and  a  passive  filter  could 
be  used.  Since  passive  filters  are  generally  far  more  linear  than  active 
filters,  nonlinear  generation  in  these  filters  did  not  present  a  problem.  In 
addition,  highly  effective  low-pass  filters  can  be  designed  passively — even 
more  effective  than  active  filters  due  to  lower  Internal  noise— when  one  can 
tolerate  the  attendant  long  turn-on  transients.  Secondly,  since  very  large 
downshift  ratios  were  used  (typically  50  thru  100),  the  filters  did  not  have 
to  be  able  to  separate  a  very  small  amplitude  difference-frequency  component 
from  the  electrically  large  amplitude  primaries  that  were  very  close  in 
frequency  to  it. 

In  order  to  assure  that  no  appreciable  difference-frequency  component  was 
being  generated  nonlinearly  in  the  active  filters  chosen  in  this  work,  a 


mixing  amplifier  was  used  to  electrically  add  two  electrical  signals  of  the 
same  frequencies  as  those  of  the  two  primaries  used  in  the  experiment.  The 
voltage  amplitudes  of  these  electrical  signals  were  chosen  to  be  comparable  to 
the  voltage  amplitudes  arising  from  the  hydrophone's  linear  response  to  the 
primary  pressures. 

No  dlf ferenc.e-f requency  component  could  be  measured  under  these  condi¬ 
tions  except  for  that  associated  with  the  noise  floor  of  the  receiver  de¬ 
vices.  This  noise  floor  corresponds  to  27  Pa,  which  Is  half  the  theoretically 
predicted  difference  frequency  in  the  forward  direction  at  just  6  cm  from  the 
center  of  the  cylinder.  The  nonlinear  signal  continues  to  grow  approximately 
linearly  in  this  direction  frc.  .his  point.  Although  it  is  clear  that  this 
noise  level  is  large  enough  to  prevent  precise  definitive  measurements  of  the 
difference-frequency  component  close  to  the  cylinder,  it  is  small  enough  to 
demonstrate  that  filter  nonlinearity  was  not  the  source  of  the  difference 
frequency  measured  (which  was  typically  about  10  times  the  theoretical  value). 

Unfortunately,  the  third  difficulty  (nonlinearity  of  the  hydrophone) 
could  not  be  eliminated.  At  the  time  the  experiment  described  here  was 
started,  no  one  even  suspected  that  hydrophones  were  non linear  to  a  measurable 
degree.  Well  after  the  start  of  the  present  experiment,  however,  a  study  of 
hydrophone  nonlinearity  was  performed  jointly  by  scientists  at  the  Naval 
Undersea  Systems  Center  and  the  Underwater  Sound  Reference  Detachment  of  the 
Naval  Research  Labortory.  The  results  of  this  study  became  available  in 
preliminary  form  [48,49].  Initially  it  was  hoped  that  the  results  (measured 
in  the  nearfield  at  a  high  downshift  ratio  for  a  piston  source  driven  at  two 
frequencies)  could  be  extrapolated  to  the  case  of  a  nonlinear  nearfield 

scattering  measurement.  This  did  not  prove  to  be  the  case.  In  fact,  when  the 

hydrophones'  nonlinear  responses  were  measured  by  a  more  appropriate  technique 

(different  than  that  used  in  Refs.  48  and  49),  they  proved  to  be  too  nonlinear 

to  make  a  correct  measurement  of  the  theoretically  predicted  difference 
frequency  generated  in  a  nonlinear  scattering  experiment.  Hence,  the 
nonlinearity  of  the  hydrophones  proved  to  be  the  limiting  factor  for  this 
experiment  prohibiting  valid  measurements  from  being  carried  out.  A 
description  of  nonlinearity  measurements  performed  on  the  hydrophones 
considered  for  use  in  this  experiment  is  provided  In  Section  IV.  F.  The  new 
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method  developed  la  this  work,  for  measuring  hydrophone  nonlinearity  and 
results  obtained  using  the  method  will  also  be  presented  in  that  section. 

D.  Selection  and  Design  of  Sound  Sources 

The  initial  factor  in  choosing  the  sound  sources  or  transducer  elements 
was  the  availability  of  a  suitable  cylindrical  transducer.  Having  such  a 
device  produced  commercially  would  have  been  time  consuming  as  well  as 
costly.  Fortunately,  standard  Navy  cylindrical  transducers  (Type  TR-127/WQM) 
were  readily  available  that  proved  to  be  adequate  for  the  experiment. 

The  TR-127/WQM  transducer  is  constructed  of  eight  cylindrical  rings  (O.D. 
=  2  cm;  I.D.  =  0.938  cm;  length  =  1.905  cm)  made  of  Type  1  Ceramic  [PZT-4,  MIL 
STANDARD  1376  (ships)]  mounted  coaxially  to  create  a  line  15.24-cra  long.  As 
used  by  the  Navy,  the  active  elements  are  contained  in  an  oil-filled,  butyl- 
rubber  boot  or  covering.  This  boot  was  removed  for  the  present  experiment  in 
order  to  avoid  any  interfering  effects  that  may  have  been  caused  by  these 
materials.  Unfortunately,  the  exposed  ceramic  became  deteriorated  by  the 
chlorine  used  to  prevent  algae  growth  in  the  test  pool.  To  avoid  this,  a 
second  cylinder  was  coated  with  Krylon,  a  commercially  available  clear  plastic 
spray  coating.  No  further  deterioration  was  noted. 

This  transducer  was  determined  to  have  a  breathing  mode  resonance 
frequency  at  102  kHz.  When  driven  at  a  signal  level  of  200-V  amplitude,  the 
cylinder  produces  an  acoustic  pressure  in  the  water  that  corresponds  to  a 
value  of  3.5  atmospheres  for  the  constant  A  in  the  asymptotic  expression  for 
the  cylinder  field: 


P  =  AH  ( l)(k  r)  ~  -  /X/r’. 
c  O  C  IT 

This  corresponds  to  a  pressure  at  the  outer  surface  of  the  cylinder 
(r  =  a  =  1  cm)  of  1.35  atm. 

Once  the  frequency  of  the  cylinder  was  chosen,  this  constrained  the 
design  of  the  piston  plane-wave  source  for  essentially  three  resons: 
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1.  In  order  to  most  effectively  filter  the  primary-wave  components  from 
the  hydrophone  electrical  response,  it  was  necessary  that  the  difference 
frequency  v  =  Vp  ”  vc  be  8*8ni f tcantly  less  than  both  the  plane-wave 
frequency  (vp)  and  the  cylindrical-wave  frequency  (v£).  (This  enables 
the  use  of  low-pass  filtering.) 

2.  The  piston  source  had  to  be  capable  of  producing  an  approximately 
planar  wave  over  a  reasonable  propagation  distance.  The  greatest 
distance  from  a  uniformly  vibrating  piston  source  where  the  wave  may  be 
considered  to  be  approximately  planar  is  at  the  position  of  the  last 
maximum  value  of  the  on-axis  piston  pressure.  This  position  occurs  at  X 
=  a2/X  where  a  =  piston  radius  [58].  Since  X  varies  inversely  as  X,  a 
small  wavelength  and  hence  a  high  frequency  is  desirable  for  this 
purpose. 

3.  As  can  be  seen  from  an  analysis  of  the  expression  for  the  difference- 
frequency  pressure  [Eq.  (115)],  the  amplitude  of  the  difference-frequency 
wave  varies  approximately  as  the  difference  frequency.  Hence,  as  the 
difference  frequency  becomes  smaller,  the  amplitude  of  the  difference- 
frequency  pressure  wave  becomes  more  difficult  to  measure. 


A  value  of  vp  ~  1.5  vc  tends  to  satisfy  simultaneously  all  three  of  the 
above  requirements.  This  results,  for  the  chosen  cylinder  frequency,  in  a 
value  of  Vp  of  approximately  150  kHz. 

It  turned  out  that  two  cylindrical  disks  of  PZT-4  were  commercially 
available  with  a  thickness  of  1.27  cm  and  a  radius  a  of  6.35  cm.  If  one  of 
the  disks  Is  mounted  so  it  Is  air  backed,  the  sound  pressure  produced  on  the 
side  opposite  from  the  air-backed  side  Is  nearly  doubled  (since  acoustic 
radiation  into  water  is  far  more  efficient  than  Into  air).  When  the  disk  is 
operated  at  its  fundamental  thickness-mode  resonance,  the  wavelength  in  the 
ceramic  Am  is  equal  to  twice  the  thickness  t  of  the  ceramic.  To  calculate  the 
corresponding  resonance  frequency  of  such  a  disk,  it  is  necessary  to  know  the 
longitudinal  sound  speed  in  the  material  in  the  thickness  direction;  i.e., 
the  Z  axis.  This  is  given  for  piezoelectric  materials  by  [59] 
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D=  tic* 

I  fc33‘ 


/P 


where  C 3^  ts  the  elastic  stiffness  at  constant  electric  displacement, 
measured  along  the  z  axis,  and  p  Is  the  mass  density  of  the  ceramic.  To 
evaluate  this  expression  as  well  as  several  others  given  below,  the  following 
values  of  certain  constants  [59]  for  PZT-4  will  be  needed: 

/e  =  635;  relative  dielectric  constant,  clamped; 

C3§  =  15.9  x  1010/m2; 

p  =  7.3  x  103  kg/m3; 

e33  =  15.1  C/ra2;  piezoelectric  stress  constant. 

Using  these  values  In  the  above  equation  for  V3°  gives  V3°  = 

4.6*103  m/sec.  This  results  in  a  resonance  frequency  of  about  181  kHz  for  the 
disk  (which  Is  close  to  the  desired  value).  Since  clamping  the  ceramic  in 
place  was  expected  to  lower  this  frequency,  the  available  ceramic  was  deemed 
adequate.  (Tne  measured  resonance  frequency  turned  out  to  be  162  kHz.) 

At  a  frequency  of  162  kHz,  the  maximum  distance  from  the  piston  face 
within  which  the  waves  could  reasonably  be  expected  to  be  planar  is  equal  to 
a2/X  ~  37  cm.  This  distance  was  considered  large  enough  to  allow  an 
experimental  test  of  the  theory.  However,  it  still  remained  to  be  seen 
whether  the  piston  could  provide  an  acoustic  wave  with  sufficient  pressure 
amplitude  (at  least  103  N/ra2)  at  a  reasonable  operating  power  and  voltage. 

The  peak  acoustic  pressure  P  produced  in  the  near field  of  an  air-backed 
piston  at  its  fundamental  thickness-mode  resonance  is  given  approximately  by 
(60! 

2e33V 


where  V  is  the  amplitude  of  the  sinusoidal  voltage  applied  across  the  two 
faces  of  the  piston.  For  a  peak  pressure  of  3xl03  N/m2,  this  expression  gives 
a  required  applied  voltage  amplitude  of  126  V. 
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To  calculate  the  acoustic  power  radiated  by  the  piston,  the  following 
equivalent  circuit  [53]  (at  resonance)  is  useful: 


O 

o 

where  the  acoustic  radiation  impedance  ZR  =  pgcS  (c  =  sound  speed  in  water, 

S  =  piston  face  area),  the  transformation  factor  a  =  e^-jS/t,  and  Cq  = 
e0e33^t’  che  usual  equation  for  the  capacitance  of  a  disk.  Using  these 
expressions  gives  Rg  =*  20.8  ohm.  The  Input  electrical  impedance  of  the  above 
equivalent  circuit  is  Z  *  (20.5  +  2.43j)  ohra;  l.e.,  the  impedance  is  primarily 
real  (resistive). 

The  average  acoustic  power  radiated  W  by  the  disk  at  an  applied  voltage 
amplltide  of  126  V  is  given  by  [60]: 

2 

2  2  V 

W  -  2a  V  /ZR  »  2^--  =  382  W. 

e 

The  corresponding  average  electrical  power  into  the  piston  is 
(384  -  5.70xl0-3j)W. 

A  power  amplifier  capable  of  producing  the  above  power  (in  a  pulsed  mode) 
and  voltage  was  available.  Therefore,  the  available  ceramic  disks  were  deemed 
appropriate  to  the  experiment. 

Unfortunately,  these  disks  suffered  the  same  damage  as  the  cylinder  when 
placed  in  chlorinated  water.  In  addition,  they  suffered  minor  pitting  due  to 
the  high  electric  field  existing  at  the  points  of  contact  of  the  electrical 
conductors.  One  of  them  was  repaired  by  filling  the  damaged  areas  with  a 
conducting  epoxy  and  by  also  epoxylng  on  two  thin  stainless-steel  circular 
plates  of  the  same  radius  as  the  ceramic  disks,  one  each  to  the  front  and  rear 
surfaces  (to  prevent  anomalously  high  electric  fields  at  the  points  of 
electrical  contact).  Although  this  procedure  worked  well  in  remedying  the 
above  problems,  the  characteristic  impedance  of  the  disk-ceramic  combination 
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was  significantly  different  than  the  ceramic  alone.  Fortunately,  the 
available  high-power  amplifier  was  capable  of  driving  the  disk  to  a  high 
enough  voltage  (360-V  amplitude)  to  produce  an  acoustic  pressure  amplitude  of 
3.16xl05  N/ra2. 

K.  Calibration  of  Selected  Hydrophones  to  Determine  First-Order  (Linear) 

Sens  it i vi t y 

From  a  rather  broad  range  of  available  hydrophones,  three  were  ultimately 
selected  for  first-order  calibration  based  on  their  inobtrusi veness  to  the 
sound  field,  as  well  as  the  expectation  that  their  nonlinearity  would  be  small 
(this  latter  aspect  will  be  discussed  in  detail  in  the  next  section).  These 
three  hydrophones  were  all  fabricated  at  the  USRD.  They  shall  be  here 
referred  to  as:  1)  small  spherical  hydrophone,  2)  F42D  hydrophone,  3)  lead 
metaniobate  hydrophone.  A  brief  description  of  each  follows. 

1 .  Small  Spherical  Hydrophone. 

This  hydrophone  was  made  of  a  0.38-cm-O.D. ,  0. 159-cm-I .D. ,  PZT-4 
spherical  shell  in  a  small  rubber  boot.  Although  the  linear  sensitivity  was 
expected  to  be  relatively  small  for  this  hydrophone,  due  to  its  relatively 
sm  1 1 1  size,  it  was  believed  this  would  be  compensated  for  by  a  very  low 
hydrophone  nonlinearity.  Also,  the  small  size  of  the  hydrophone  assured  that 
it  would  be  especially  inobtrusive  to  the  measured  sound  field. 

2 .  F42D  Hydrophone. 

The  F42D  is  a  standard  hydrophone  available  for  use  by  customers  of  the 
USRD.  Its  active  element  is  PZT-4  ceramic,  configured  in  a  spherical  shell 
design  with  an  O.D.  of  1.28  cm  and  an  I.D.  of  1.08  cm.  The  spherical  shell  is 
encapsulated  in  polyurethane.  Although  somewhat  larger  than  the  small 
spherical  hydrophone,  it  is  not  unacceptably  so.  Its  larger  size  gives  a 
corresponding  increase  in  first-order  sensitivity,  easing  analysis  of  the 
received  signals. 


3 .  Lead  Metaniobate  Hydrophone. 

This  hydrophone's  active  element,  as  Its  name  implies,  is  made  of  lead 
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mot -inlo bate  fashioned  Into  a  cylindrical  shape  and  encapsulated  in  a 
transparent  rubber  boot.  Being  of  single-piece  design,  it  was  believed  the 
nonlinearity  of  this  hydrophone  would  be  exceptionally  low  (since  glue  joints 
and  other  bonds  can  be  contributing  factors  to  this  nonlinearity). 
Additionally,  being  of  approximately  the  same  diameter  as  the  F42D  hydrophone, 
it  also  was  expected  to  have  a  higher  first-order  sensitivity  than  the  small 
spherical  hydrophone,  as  well  as  being  sufficiently  inobtrusive  acoustically. 

The  first-order  (linear)  sensitivity  of  a  hydrophone  in  V/pPa  is  defined 
as  the  ratio  of  the  voltage  developed  across  the  open-circuited  hydrophone 
terminals  when  a  plane  wave  is  incident  on  the  hydrophone  to  the  acoustic 
pressure  of  the  plane  wave.  The  sensitivity  is  a  complex  quantity;  it 
contains  both  an  amplitude  and  a  phase  angle.  The  phase  angle  is  not  required 
for  the  present  work.  The  sensitivity  is  also  a  function  of  frequency.  It 
must  be  separately  determined  for  each  frequency  component  that  is  to  be 
measured.  It  can  also  depend  on  the  direction  of  the  incoming  plane  wave, 
especially  at  higher  frequencies.  For  low  frequencies,  however,  it  is  often 
uniform  over  all  4x  steadians.  The  method  chosen  for  first-order  calibration 
was  one  recently  developed  at  the  USRD  161]-  Although  phase  calibration  was 
not  required  fnr  this  experiment,  this  procedure  also  gave  the  magnitude  of 
the  first-order  sensitivity,  was  simple  to  implement,  and  the  rigging  and 
necessary  equipment  were  readily  available. 

Since  the  frequencies  of  interest  were  already  established  by  the 
selection  of  the  cylinder  and  piston  at  102,  162,  and  60  kHz  (difference 
frequency),  these  hydrophones  had  to  be  calibrated  only  at  these 
frequencies.  Table  II  gives  the  results  of  this  calibration. 

F.  Determination  of  Hydrophone  Nonlinearity 

l .  General  Considerations. 

In  two  recent  reports  [48,49]  it  is  noted  that  the  response  of  ceramic 
hydrophone  elements  is  not  strictly  linear  but  may  be  more  accurately 
represented  by  the  parabolic  relationship 
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where  c  --  voltage  produced  by  the  hydrophone  In  pressure  field  P 

P  -=  pressure  wave  incident  on  the  hydrophone 

m  =  first-order  (linear)  response  of  the  hydrophone 
(in  V/uPa) 

n  =  second-order  (nonlinear)  response  of  the  hydrophone  [in 
V/(pPa)2]. 

They  define  the  first-order  sensitivity  (M)  in  dB  re  1  V/pPa  as 

M  =  20  log(m)  (119) 

and  the  second-order  sensitivity  in  dB  re  1  V/(pPa)2  as 

k  =  20  log(n).  (120) 


Table  II.  Hydrophone  Linear  Sensitivities 


HYDROPHONE 

FREQUENCY  SENSITIVITY 

SENSITIVITY 

(kHz) 

(V/pPa) 

(dB  re  1  V/pPa) 

162 

1.45xl0"12 

-236.8 

SMALL  SPHERE 

102 

1.5 5* 10” 12 

-236.2 

60 

1.33xlO-12 

-237.5 

162 

4.35xlO-11 

-207.2 

F42D 

102 

2.48xlO-11 

-212.1 

60 

2.91xlO-11 

-210.7 

LEAD 

162 

1 . 54x10“* 1 

-216.2 

MF.TANIOBATE 

102 

2.20xl0_u 

-213.2 

60 

5.80xlO-12 

-224.7 

The  unit  dB  refers  to  a  decibel ,  which  in  this  case  is  20  times  the 
logarithm  to  the  base  10  of  the  sensitivity  divided  by  the  reference 
sensitivity  of  1  V/pPa. 
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In  the  treatment  given  in  the  Moffett-Blue  and  Mof fett-Henriquez  reports 
for  incident  waves  containing  two  primary  components,  the  primary  pressure 
waves  are  assumed  to  be  of  equal  amplitude.  Since  this  was  clearly  not  to  be 
the  case  in  the  present  experiment,  the  theory  is  generalized  here  for  the 
case  of  unequal  primary  pressure  fields. 

Let  the  total  pressure  wave  at  the  hydrophone  due  to  primary  waves  of 
angular  frequencies  id  ^  ,0)2  and  amplitudes  ,P2  be  represented  as 

P  =  P^cos  (tD^t )  +  P2cos(<D2t).  (121) 

The  contribution  to  the  quadratic  term  in  Eq.  (118)  due  to  difference- 
frequency  components  may  be  calculated  by  squaring  Eq.  (121)  and  suitably 
identifying  terms.  By  use  of  the  trigonometric  identity  cos9  +  cos<J>  =  2[cos 
(9-<t>)/2]  [cos  (9+$)/2],  P^  may  be  readily  put  into  the  form 

P2  =  (Pt2/2  )[1  +  cos(2(Dlt)]  +  (P22/2)[l  +  cos(2(D2t)] 

+  riP2[cos((Dl+<D2)t  +  cos((*)l“ o»2)t] .  (122) 


According  to  Eq.  (118)  the  contribution  to  the  voltage  produced  by  the 
hydrophone  due  to  its  self  nonlinearity  is  nP^.  The  contribution  to  this 
voltage  at  the  difference  frequency  may  be  readily  seen  in  Eq.  (122)  to  be 
r|^>l*>2  since  this  is  the  coefficient  of  the  difference-frequency  term.  This 
Indicates  that  the  pressure  amplitudes  for  each  of  the  two  primary  waves  must 
be  measured  in  addition  to  the  voltage  produced  at  the  difference  frequency 
due  to  the  hydrophone  nonlinearity.  Only  then  can  the  nonlinear  response  n  of 
the  hydrophone  be  determined. 

In  the  experimental  portion  of  the  Moffett-Blue  and  Mof fett-Henriquez 
research,  the  two  primary  frequencies  were  chosen  close  together  to  produce  a 
large  downshift  ratio  (the  downshift  ratio  is  the  quotient  of  the  average 
value  of  the  primary  frequencies  to  the  difference  frequency).  This  avoids 
several  sources  of  error  in  the  nonlinear  measurement  (to  be  described  in 
detail  later  on).  Unfor-tunately ,  for  the  frequencies  of  Interest  in  the 
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current  nonlinear  scattering  measurement,  the  downshift  ratio  is  only  about 
2.2,  com-pared  with  more  usual  values  of  50  to  100.  In  order  to  circumvent 
difficulties  associated  with  this  unusually  small  downshift  ratio,  an 
approximate  approach  to  the  measurement  of  n  and  K  was  adopted. 

The  following  quantities  are  useful  In  describing  the  approach  taken: 

K162,nl62:  The  hydrophone  second-order  sensitivity  at 
162  kHz  and  the  second-order  response  at  a 
2-kHz  difference  frequency  In  the  presence 
of  primaries  at  161  and  163  kHz. 

K102,r*l02:  The  hydrophone  second-order  sensitivity  at 
102  kHz  and  the  second-order  response  at  a 
2-kHz  difference  frequency  In  the  presence 
of  primaries  at  101  and  103  kHz. 

K60,n6o  the  hydrophone  second-order  sensitivity  at  60  kHz  and  the 

second-order  response  at  a  60-kHz  difference  frequency  in  the 
presence  of  primaries  at  102  and  162  kHz. 


It  was  initially  assumed  that  would  be  approximate! y  the  average  of  Kjg? 
and  i*e.  , 


v  -  ^  1 02  _+  162 

*60  2 


(123) 


or  equivalently  that 


/nl02n162 


(124) 


Equations  (123)  and  (124)  are  reasonable  since  the  nonlinear  effect  Is 
quadratic.  Measurement  of  0jf)2  nn<*  ^1^2  woll'd  not  prove  as  difficult  is 
directly  measuring  n^Q,  because  each  of  these  has  a  reasonably  large  downshift 
ratio  (51  and  81,  respectively).  K[02  and  K162  iirc  directly  calculable  from 
r'l02  an<*  nl62  hy  the  relationship  K  =  20  Iog|nj,  and  Eqs.  (123)  and  (124)  can 
be  used  to  obtain  K^q. 

There  are  essentially  three  sources  of  error  in  the  nonlinear  calibration 
of  a  hydrophone,  all  of  which  are  minimized  by  the  selection  of  a  large 
downshift  ratio.  Each  of  these  will  now  be  considered  in  detail. 


a.  Generation  of  difference  frequency  by  nonlinear  interaction  of  the 
primaries  In  the  water  -  The  entire  basis  for  the  present  scattering 
experiment  Is  the  fact  that  when  sound  waves  of  two  different  frequencies 
are  simultaneously  present  In  water,  a  difference-frequency  component  is 
generated.  The  problem  in  a  nonlinear  hydrophone  calibration  is  to 
minimize  this  so  that  only  the  difference-frequency  components  generated 
by  a  hydrophone’s  self  nonlinearity  are  present  in  the  hydrophone 
electrical  output.  This  can  be  done  in  two  ways.  First,  It  is  known 
that  the  difference-frequency  component  generated  In  the  water  tends  to 
grow  with  distance  from  the  source.  Therefore,  to  minimize  the  effect, 
the  hydrophone  should  be  placed  as  close  as  practicable  to  the  source. 

The  problem  of  calculating  the  difference-frequency  pressure  generated  in 
the  extreme  nearfield  of  a  piston  source  is  not  trivial.  However,  the 
source  levels  can  be  estimated  from  data  presented  by  Moffett  and  Mellen 
[62].  Use  of  this  data  indicated  that  the  levels  generated  for  the 
selected  downshift  ratios  would  be  on  the  same  level  as  pseudosound  (see 
Section  F.  1.  c  below),  if  the  hydrophone  is  placed  5-10  cm  from  the 
source.  Hence,  this  source  of  error  can  be  made  negligibly  small. 


b .  Direct  radiation  at  the  difference  frequency  by  the  piston  source  - 
Since  any  amplifier  will  be  nonlinear  to  some  extent,  it  can  be  expected 
that  when  time-varying  signals  are  applied  to  the  amplifier  input  at  two 
different  frequencies  that  a  voltage  will  appear  at  the  output  at  the 
difference  frequency.  This  will  in  turn  be  applied  across  the  piston 
source  and  will  be  directly  radiated  into  the  water.  This  must  be 
calculated  beforehand,  therefore,  to  determine  if  it  can  have  a  serious 
effect  on  the  received  voltage  at  the  difference  frequency. 

In  order  to  do  so,  we  first  note  that  the  equation  for  the  normal 
surface  velocity  of  an  air-backed  piston  face  in  contact  with  a  fluid 
medium  (chosen  to  be  water)  is  given  by  [63] 


2aV 


Ze  "  JC3“S/(u)t) 

where  V  “  voltage  applied  across  the  piston 
S  ■  surface  area  of  the  piston  face  = 

<u  -  angular  frequency  of  the  applied  voltage 


(125) 


t  =  piston  thickness 

Ze  =  acoustic  radiation  impedance 

and  where  a  and  have  the  same  meaning  as  in  Section  D  of  this 

chapter.  The  expression  for  the  acoustic  radiation  impedance  of  a  piston 
face  is  given  by  [64]: 

Ze  =  Sp  c(0  -  ix  ) 

o  o  o 

12  2 

V2  k  a  ka  +  0 
1  ka  +  “ 

8ka/(3it)  ka  ♦  0 

X  *  (126) 

2/(nka)  ka  +  00 

(■here  pQ  is  the  water  density. 

It  is  clear  from  Eq  .  (126)  that  the  impedance  becomes  very  small 
relative  to  SpQC  as  ka*0. 

Since  in  Eq .  (125)  this  impedance  is  combined  with  another  term 
(- jC3§S/<Dt ) ,  which  is  much  larger,  the  impedance  term  is  negligible.  Equating 
it  to  zero  in  Eq .  (125)  we  obtain 

u  =  2out jV/(C3°S).  (127) 

Equation  (127)  can  be  combined  with  the  expression  for  the  pressure 
radiated  a  distance  r  by  a  periodic  simple  source  [65] 

P  .  S  eik(r'ct)  (128) 

4itr  (11 

where  the  volume  velocity  of  the  source  is  equal  to  uS  in  the  present 
case.  This  produces  the  radiated  pressure  per  unit  applied  voltage 
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(129) 


where  a  has  been  replaced  by  the  expression  given  for  It  In  Section  IV. 

D.  Using  a  frequency  of  2  kHz  (as  would  be  present  In  measuring  or 

nlfe2) »  and  a  distance  of  1  cm,  Eq.  (129)  yields  a  value  of  3.02  Pa/V. 

The  measured  amplifier  output  voltage  at  the  difference  frequency  was  at 
most  6  V.  This  gives  a  maximum  radiated  pressure  of  about  18.1  Pa.  To 
determine  what  effect  this  would  have  on  the  measurement,  the  first-order 
sensitivities  (Table  II)  must  be  used  to  calculate  the  voltages  produced 
at  the  difference  frequency  by  this  pressure.  They  are 

Small  Sphere:  2.40  x  10”^  V 
F42D:  5.30  *  10~4  V 
Lead  Metaniobate:  1.05  x  10-4  V. 

These  voltages  are  comparable  to  the  measured  voltages  at  the  difference 
frequency.  Although  this  indicates  that  a  directly  radiated  difference- 
frequency  component  Is  to  a  certain  extent  significant,  it  nonetheless 
does  not  invalidate  the  nonlinearity  calibrations  but  rather  sets  a  bound 
below  which  the  nonlinearity  of  the  hydrophones  in  question  cannot  he 
accurately  determined.  This  is  sufficient  for  the  purpose  of  this 
experiment . 

c.  Pseudosound  -  In  any  acoustical  measurement,  a  hydrophone's  motion  is 
uncertain  by  an  amount  of  the  same  order  as  the  difference  between 
Lagrangian  and  Euler lan  coordinates  (since  the  hydrophone  cannot  be 
completely  free  to  move  with  the  fluid  nor  can  it  be  completely  rigidly 
held  in  a  fixed  position).  Hence,  there  Is  an  uncertainty  in  the  meaning 
of  the  value  of  the  measured  acoustic  pressure  that  corresponds  to  the 
difference  between  the  Lagrangian  and  Kulertan  frames  of  reference.  This 
Is  presently  a  fundamental  experimental  limitation  In  that  it  cannot  be 
fully  eliminated  using  present-day  technology  (although  It  can  be  reduced 
by  using  a  massive  receiver  well  below  Its  lowest  mechanical 
resonance).  The  difference  in  acoustic  pressure  between  a  Lagrangian  and 
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Ku  I  <•  r  lan  frame  of  reference  is  known  as  pseudosound .  An  estimate  can  ha 
made  of  the  level  of  pseudosound  by  expanding  any  Lagrangian  quantity 
qL(a,t)  in  terms  of  the  appropriate  Eulerian  quantity  qE(x,t)  as  (66) 


qL(a,t)  =  qE(x,t) 


x=a-K(a ,t ) 


q  (x.t) 


x=a 


3q  (x,t) 
3x 


x=a 


C(x  ,t )  +  ... 


(130) 


For  definitions  of  a,  x,  £  see  Section  II.  B.  Using  Eq.  (130)  to  expand 
the  pressure  gives 


p‘-(a,t) 


PE(x,t) 


+  (x,t) 

x=a  3x 


x=o 


C(x,t). 


(131) 


To  obtain  an  upper  bound  for  pseudosound,  the  pressure  present  may  he 
represented  as  plane  waves  as  follows 


PE(x,t)  =  P^x.t)  +  Pj  (x,t) 


(132) 


where  Pi0>P20  are  rea^  consfants*  Equation  (132)  may  be  combined  wit 
Eq.  (51)  to  obtain  the  following  expression  for  S(x,t)  in  terms  of 
pressure 


f.(x,t)  =  f.,(x,t)  +  f,2(x,t) 
Ik 


lk„ 


Re  [ - Z"El  (  x  » 1 + - 2-P2  (x,t)). 


(133) 
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Equation  (132)  may  also  be  used  to  express  [ 3PE(x ,t ) ] /3x  as 


g* 

-P-  =  Re(ik1PE(x,t)  +  ik2P2R(x,t)] 


(134) 


Hence,  pseudosound  (the  difference  between  the  Lagrangian  and  Eulerian 
expression  for  the  pressure)  may  be  written 


PL  ~  PE  -  RelikjP^Cx.t)  +  ik2P2E(x,t)l 


k  P  E(x,t)  k  P  E(x,t) 
Reli(-i-i-j - +  -----y - )]. 


po“l 


PoW2 


(135) 


At  this  point,  it  is  helpful  to  make  a  few  definitions: 


h  =  ikipi 
al  5  1/(P0W12) 

&2  s  l/(Pou2 

Equation  (135)  may  now  be  expressed  as 
PL  -  PE  =  Rc(Z1+Z2)Re(a1Zl  +  a2Z2) 

=  y  (Zj  +  Zx*  +  Z2  +  Z2*)(alZl  +  a^*  +  a2Z2  +  a^*). 


(136) 


Contributions  to  the  difference  frequency  can  arise  only  from  the  cross¬ 
terms;  therefore,  pseudosound  at  the  difference  frequency  is  given  by 


(PL  -  PE)dlff  *  \  <a2ZiZ2*  +  a2Zl*Z2  +  alZ2Zl*  +  alW 


or 


<pL  -  pE)diff  ■  TT  (  A2  +  -2  >  P10P20' 

o  (2>l  a»2 


(137) 


(138) 


Using  the  experimental  values  involved  In  the  nonlinear  calibration  fn 
Eq.  (141)  yields  less  than  0.1  Pa  in  both  measurements.  Hence, 
pseudosound  is  a  rather  negligible  source  of  error  in  this  measurement. 
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2 .  Measurement  Setup  and  Results. 

Figure  29  is  a  schematic  representation  of  the  measurement  configuration 
used  in  the  nonlinear  calibration.  Table  III  gives  the  results.  The 
quantities  n^Q  and  K^q  are  the  values  obtained  using  F,qs.  (123)  and  (124). 

The  results  shown  In  Table  III  can  be  misleading.  Apparently  the  small 
spherical  hydrophone  Is  less  nonlinear  than  the  F42D  hydrophone  (due  to  Its 
smaller  value  of  H^q).  However,  cannot  be  considered  independently  of  the 

first-order  sensitivity  m.  An  important  quantity  in  determining  the 
desirability  of  a  hydrophone  for  use  in  a  nonlinear  experiment  is  the  a pp i ren t 
pressure  present  at  the  hydrophone  at  the  difference  frequency  that 
corresponds  to  the  electrical  signal  produced  by  the  nonlinearity  of  the 
hydrophone.  Let  the  ratio  of  this  quantity  to  the  product  of  the  primary 
pressures  be  called  the  "nonlinearity"  and  be  given  the  symbol  D.  The 
apparent  pressure  (in  Pa)  may  be  calculated  from  the  product  of  the  primaries 
in  terms  of  the  first-  and  second-order  sensitivities  as 

106nP.P, 

P  « - .  (139) 

app  m 

Therefore,  the  quantity  D  (measured  in  units  of  apparent  Pa  per  squared  Pa  ot 
primary)  may  be  calculated  from 


D  =  I06n/m.  (14n; 

[The  factors  of  10^  in  Eqs.  (139)  and  (140)  above  are  required  to  render 
consistent  the  units  chosen.)  Using  the  data  in  Tables  II  and  III  in  Eq . 

(140)  gives  for  the  small  spherical  hydrophone  and  the  F42D  hydrophone  the 
following  results: 

DSmall  =  4.74  x  10-8  _ 8a(aPP> 

Sphere  Pa2(primary)  (141) 

°F42D  *  5.43  *  10"9  - P«(?P-P>- - 

Pa^( primary) .  (142) 
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Fig.  29  -  Measurement  configuration  for  nonlinear  hydrophone 
calibration  via  the  Moffett~Blue  method 
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Table  III.  Nonlinear  Calibration  of  Hydrophones  by  Moffett-Blue  and 


Although  n6c,  for  the  F42D  hydrophone  Is  about  2.3  times  greater  than  that  for 
the  small  spherical  hydrophone,  it  Is  not  the  more  nonlinear  of  the  two. 
Results  in  Eqs.  (141)  and  (142)  show  that  the  small  spherical  hydrophone  is 
more  than  8.7  times  more  nonlinear  than  the  F42D  hydrophone. 

Although  the  value  of  D  for  the  F42D  hydrophone  is  apparently 
sufficiently  low  to  perform  the  nonlinear  scattering  measurements,  preliminary 
results  gave  anomalously  large  values  for  the  difference  frequency.  This 
indicated  that  the  values  of  D  obtained  by  the  above  method  are  not  applicable 
to  this  experiment  for  two  possible  reasons: 

a.  The  value  of  D  obtained  by  the  Moffett-Blue  and  Mof fet-Henrique2 
approaches  at  high  downshift  ratio  may  not  apply  to  low  downshift  ratio 
(indicating  the  hydrophone  nonlinearity  is  extremely  sensitive  to 
downshift  ratio). 

b.  The  hydrophone  nonlinearity  is  a  function  of  the  angle  at  which  the 
primaries  intersect  at  the  hydrophone  (in  the  Moffett-Blue  and  Moffett- 
Henriquez  approaches,  this  is  always  exactly  0°;  but  in  a  near  field 
scattering  experiment,  a  wide  continuum  of  angles  is  simultaneously 
present). 

In  order  to  determine  whether  either  or  both  of  the  above  possibilities  were 
present  during  the  experiment ,  a  new  approach  to  determine  hydrophone  non¬ 
linearity  (which  more  closely  approximated  the  experiment  of  interest)  was 
attempted.  The  geometry  of  this  measurement  is  shown  in  Fig.  30,  and  a  sche¬ 
matic  representation  of  the  measurement  configuration  is  presented  in  Fig.  31. 

In  this  measurement,  the  hydrophone  of  interest  is  placed  at  a  distance  X 
from  the  center  of  the  active  surfaces  of  the  piston  and  the  cylinder. 

Acoustic  pulses  are  sent  from  each  source  that  are  of  sufficient  lengths  to 
overlap  at  the  position  of  the  hydrophone  at  the  moment  of  measurement  (i.e., 
of  length  just  greater  than  X).  This  ensures  that  no  difference  frequency 
will  be  generated  by  the  fluid  medium,  since  the  waves  will  not  have  had  time 
to  overlap  In  the  fluid.  (Hence,  all  measured  difference  frequency  will  be 
generated  in  the  hydrophone  itself.)  The  piston  and  cylinder  were  driven  at 
the  same  frequencies  as  were  to  be  used  in  the  nonlinear  scattering 
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Fig.  30  -  Geometry  of  new  method  of 
nonlinear  hydrophone  calibration 
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ex  per linent  ,  L62  and  102  kHz,  respectively.  The  product  of  the  pressure  levels 


of  the  primaries  was  also  similar  to  that  in  the  nonlinar  scattering 
experiment,  specifically  10^®  Pa^.  The  values  of  D  obtained  using  this  method 
on  the  same  hydrophones  tested  using  the  Moffett-Blue  and  Mof fett-Henriquez 
approaches  are  presented  in  Table  IV. 


Table  IV.  Values  of  D  Obtained  from  Second  Hydrophone 
Nonlinearity  Calibration  Method 


Hydrophone 

^  Pa (apparent) 

(Pa2)  (primary) 

Small  Sphere 

2.47  x  10-7 

F42D 

6.51  x  10“8 

Lead  Metaniobate 

5.34  x  10-8 

(Values  accurate  to  within  7%  as  determined  by  reproducibility 
of  experimental  data) 


It  can  be  seen  from  Table  IV  that  the  results  were  considerably  greater 
in  all  cases  than  the  values  obtained  by  the  Moffett-Blue  and  Moffett- 
Henriquez  methods.  In  an  attempt  to  find  a  suitable  hydrophone,  several 
others  were  tested.  The  results  are  given  in  Table  V. 
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Table  V.  Values  of  D  for  Some  Additional  Hydrophones 


Pa  (a 

pparent ) 

Hydrophone 

D 

(Pa2) 

(primary) 

F42C30 

4.81 

oo 

1 

2 

X 

F42D36 

5.92 

X  10~8 

F42D14  ' 

7.40 

X  10~8 

(Values  accurate  to  within  7%  as  determined  by  reproducibility 
of  experimental  data) 

Since  the  product  of  the  pressure  values  of  the  primaries  in  the 

10  7 

nonlinear  scattering  experiment  everywhere  gave  values  of  order  10  Pa  and 
since  typical  pressures  expected  theoretically  at  the  dif terence-f requency 
component  were  of  the  order  of  100  Pa,  it  was  clear  that  no  hydrophone 
available  was  sufficiently  linear  to  perform  the  experiment. 

G.  Suggestions  for  Future  Measurements 

As  we  have  seen,  it  is  the  hydrophone  self-nonlinearity  that  precludes  a 
successful  measurement  of  the  difference-frequency  pressure  generated  in  the 
water  when  a  plane-wave  scatters  from  a  vibrating  cylinder.  Until  such  time 
as  a  sufficiently  linear  hydrophone  is  available,  it  would  appear  that 
confirmation  of  the  theory  is  not  possible.  However,  if  a  method  can  be  found 
for  Increasing  the  difference-frequency  pressure  generated  in  the  water  while 
maintaining  the  product  amplitude  of  the  primary  waves  approximately  constant, 
it  may  he  possible  to  make  the  difference-frequency  voltage  generated  by  the 
hydrophone  self-nonlinearity  a  small  part  of  the  total  difference-frequency 
signal.  It  will  be  recalled  that  during  the  discussion  regarding  the  choice 
of  the  experiment,  it  was  remarked  that  having  the  difference  frequency  lower 
than  either  primary  frequency  was  highly  desirable  because  it  greatly 
simplified  the  necessary  electrical  filtering  (as  well  as  avoiding  harmonics 
of  the  primaries).  Unfortunately,  due  to  the  severe  constraints  on  the 
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possible  frequency  of  the  cylindrical  source  (recall  that  100  kHz  was  judged 
to  be  about  the  maximum  possible),  the  difference-frequency  choice  was  also 
severely  constrained.  The  60-kHz  difference  frequency  chosen  was  about  the 
greatest  possible  that  also  allowed  electrical  filtering  to  distinguish  it 
from  the  primaries,  and  yet  it  resulted  in  a  pressure-amplitude  too  small  to 
measure  once  the  hydrophone  nonlinearity  was  discovered.  However ,  it  may  bo 
possible  to  choose  a  difference  frequency  above  the  cylinder  frequency  and  yet 
sufficiently  removed  from  the  higher  cylinder  harmonics,  such  that  their 
disturbing  Influence  on  the  measurement  becomes  acceptably  small.  For 
example,  if  we  retain  a  100-kHz  cylinder  frequency  but  allow  a  1-MHz  plane 
frequency,  the  difference  frequency  of  900  kHz  represents  the  ninth  harmonic 
of  the  cylinder.  For  a  cylinder  pressure-amplitude  of  10^  Pa,  the  nintj? 
harmonic  is  approximately  6  Pa  (in  the  plane-wave  case).  However,  since  the 
difference  frequency  in  this  case  scales  approximately  as  the  square  root  of 
the  primary  frequencies  and  directly  as  the  difference  frequency  (see  the 
expression  for  the  farfield  limit  given  in  Chapter  II),  the  difference- 
frequency  pressure  should  be  approximately  37.5  times  greater  than  in  the  case 
experimentally  studied  here.  This  gives  (for  directions  away  from  the 
forward-scattering  direction)  a  difference-frequency  pressure  of  about 
900  Pa.  Of  course,  even  this  suggested  measurement  is  not  without  great 
difficulties.  The  worst  of  these  are:  l)  The  900-kHz  difference  frequency  is 
only  about  1  dB  below  the  1-MHz  planar  primary  frequency,  resulting  in  even 
greater  electrical  filtering  difficulties.  2)  Once  again  the  plane-wave 
frequency  is  above  the  resonance  frequency  of  the  cylinder.  Hence,  the 
cylinder  will  not  act  as  a  rigid  body,  and  the  primary  field  will  again  have 
to  be  determined  empirically.  This  is  far  more  difficult  to  do  in  this  case, 
however,  since  the  wavelength  of  the  scattered  primary  field  (~  0.15  cm)  is 
becoming  small  compared  with  any  currently  available  hydrophone.  Hence, 
accurate  phase  measurements  needed  to  determine  the  scattering  coefficients 
would  be  very  difficult  to  perform. 

Another  possibility  is  to  operate  the  cylinder  at  a  higher  frequency 
(using  a  higher  harmonic  or  a  different  vibrational  mode).  Although  the 
radiated  pressure  amplitude  of  the  cylindrical  wave  will  be  substantially 
lower,  if  a  significant  Increase  in  difference-frequency  pressure  is  possible 
due  to  frequency  scaling,  then  a  more  favorable  overall  experimental  condition 
might  result.  For  example,  if  the  cylinder  can  be  operated  at  a  frequency  of 


I  MHz,  and  the  plane-wave  frequency  Is  chosen  to  be  1.5  MHz,  the  25-Pa 
difference-frequency  pressure  computed  for  the  present  research  would  scale  to 
over  2000  Pa.  Hence,  even  If  the  pressure  amplitude  of  the  radiated 
cylindrical  wave  is  two  orders  of  magnitude  less  than  that  at  its  102-kHz 
resonance,  the  difference-frequency  pressure  will  remain  comparable  to  that 
predicted  with  the  parameters  chosen  for  the  experimental  measurement 
attempted  here.  This  would  result  In  approximately  two  orders  of  magnitude  ol 
improvement  in  the  ratio  of  the  nonlinear  fluid  signal  to  the  hydrophone ' s 
self-nonlinearity  signal.  On  the  other  hand,  20  Pa  Is  still  not  a  very  large 
signal  and  may  remain  below  the  noise  floor  of  the  receiving  equipment 
(although  measurement  at  angles  near  0=0°  might  nonetheless  be  possible 
since  the  difference-frequency  pressure  at  these  favorable  angles  may  be  of 
sufficient  amplitude  to  measure).  Of  course,  if  no  cylinder  resonance 
frequency  can  be  found  that  generates  a  cylindrical  primary  pressure  amplitude 
greater  than  an  amplitude  three  orders  of  magnitude  less  than  that  produced  by 
the  102-kHz  frequency  used,  it  is  doubtful  even  at  0  =  0°  that  a  sufficient 
propagation  distance  exists  to  simultaneously  yield  difference-frequency 
pressure  amplitude  large  enough  to  measure;  and  also  satisfy  the  constraints  of 
Lh<  theory  In  regard  to  the  region  of  validity. 

Thus,  we  see  that  there  Is  no  obvious  or  straightforward  way  to  improve 
the  experimental  measurement.  However,  it  may  nonetheless  be  possible  to  in 
some  manner  take  advantage  of  the  direct  scaling  of  the  difference-frequency 
pressure  with  the  difference  frequency  to  avoid  the  significant  Interference 
of  the  difference-frequency  signal  associated  with  hydrophone  self 
nonlinearity. 

V.  CONCLUSIONS 

The  objectives  of  tills  research  were  threefold:  l)  To  re-derivo  tin* 
simple-source  formulation  of  the  second-order  nonlinear  wave  equation  for 
arbitrary  primary  fields.  2)  To  Investigate  theoretical ly  the  solutions  oi 
this  equation  for  three  cases  Involving  the  scattering  of  acoustic  waves  by 
vibrating  obstacles.  3)  To  subject  the  case  of  plane-wave  scattering  by  a 
vibrating  cylindrical  obstacle  to  an  experimental  analysis.  The  first  two  of 
these  objectives  were  successfully  achieved,  but  accomplishment  of  the  third 
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objective  was  Inhibited  by  hydrophone  nonllnearltles  that  were  larger  than 
were  previously  thought  to  exist. 

In  deriving  the  simple-source  formulation  of  the  second-order  nonlinear 
wave  equation,  care  was  taken  to  avoid  using  perturbation  analysis  until  the 
final  step.  While  the  acuracy  of  a  solution  obtained  via  a  perturbative 
approach  becomes  questionable  when  second-order  fields  become  comparable  to 
first-order  fields,  this  Is  not  so  for  analyses  based  on  the  Inherent  physical 
dependence  of  the  acoustic  variables  on  Mach  number.  Since  the  derivation  of 
the  simple-source  formulation  of  the  second-order  nonlinear  wave  equation 
presented  here  was  based  on  the  Inherent  dependence  of  the  variables  on  the 
Mach  number  (In  all  but  the  final  step),  the  validity  of  all  second-order 
equations  up  to  this  point  are  not  restricted  when  the  second-order  quantities 
become  comparable  to  the  first-order  quantities.  Also,  since  third-order 
acoustic  quantities  are  never  significant  in  water  and  since  no  restriction 
was  placed  on  the  geometry  of  the  fields,  this  equation  may  be  viewed  as  being 
valid  In  a  very  general  sense. 

The  problem  of  nonlinear  scattering  of  acoustic  waves  by  vibrating 
obstacles  was  solved  via  a  perturbation  solution  of  the  simple  source 
formulation  of  the  second-order  nonlinear  wave  equation  for:  1)  plane-wave 
scattering  by  a  vibrating  plane,  2)  cylindrical-wave  scattering  by  a  vibrating 
cylinder,  and  3)  plane-wave  scattering  by  a  vibrating  cylinder.  Since  the 
solutions  were  obtained  via  a  perturbation  method,  they  are  restricted  to 
situations  in  which  the  second-order  fields  remain  small  relative  to  the 
first-order  fields.  This  means  the  solutions  are  restricted  to  the  nearfield 
of  the  scattering  obstacles.  It  was  further  demonstrated  that  the  solutions 
to  this  problem  obtained  via  the  Censor  approach  for  the  sum-  and  difference- 
frequency  pressures  are  of  the  order  of  pseudosound.  Since  the  solutions  of 
the  second-order  nonlinear  wave  equation  for  the  sum-  and  difference-frequency 
pressures  tend  to  grow  with  increasing  distance  from  the  scatterer's  surface, 
they  overhelm  the  effect  predicted  by  Censor  within  a  fraction  of  a  wavelength 
of  propagation  distance  (as  was  conjectured  by  Rogers  (24]).  Graphical 
results  were  presented  in  the  case  of  plane-wave  scattering  from  a  vibrating 
cylinder  for  both  Censor's  theory  and  the  nonlinear  theory. 

Although  a  successful  comparison  of  theory  and  experiment  was  not 
achieved,  several  significant  observations  were  made.  No  previous 
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experimenters  had  performed  measurements  of  nonlinearly  generated  difference- 
frequency  acoustic  signals  in  the  extreme  nearfield  of  the  sources.  Several 
potential  sources  of  error  in  such  measurements  were  identified.  They 
included:  1)  direct  radiation  of  difference-frequency  pressure  by  the 

sources,  2)  electrical  filtering  problems  due  to  experimental  constraints,  and 
3)  difference-frequency  voltage  generated  nonlinearly  in  the  hydrophone. 
Although  solutions  to  the  first  two  difficulties  were  found,  the  third 
difficulty  proved  to  be  unresolvable .  However,  several  significant  measure¬ 
ments  of  hydrophone  nonlinearity  were  made  via  a  new  technique.  This  new 
technique  represents  a  significant  advance  over  the  technique  of  Moffett-Blue 
and  Mof fett-Henriquez.  As  of  the  writing  of  this  thesis,  a  positive  effect  of 
the  new  nonlinearity  measurement  presented  herein  has  been  a  significant 
expansion  of  the  effort  to  develop  a  linear  hydrophone  at  the  USRD. 

It  is  hoped  that  this  work  has  firmly  established  the  rather  general 
validity  of  the  simple-source  formulation  of  the  second-order  nonlinear  wave 
equation.  Similarly,  it  is  hoped  that  the  new  technique  of  hydrophone 
nonlinearity  calibration  developed  during  the  course  of  this  work  will  form 
the  basis  for  a  standard  method  of  second-order  hydrophone  calibration. 
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APPENDIX 


Integrals  of  the  form  (AI)  are  considered  by  several  other  authors  [67- 
74]  when  the  functions  R^^(x)  are  Bessel  functions.  In  relations  A2 ,  aM,  bu , 
Cy,  and  d ^  are  known  functions  corresponding  to  The  symbol  D  represents 

d/dx.  The  function  f(x)  and  the  product  nR^  are  both  assumed  bounded  and 
continuous  (or  with  at  most  a  finite  number  of  discontinuities)  over  an 
interval  [x^^],  insuring  that  the  integral  I  exists  in  the  same  interval. 

Recurrence  relations  [Eq.  (A2)]  may  be  combined  to  show  that  the 
functions  satisfy  the  differential  equation 
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Equation  (A3)  la  a  special  case  of  the  Sturm-Liouville  differential  equation 


D[p(x)DKx)]  +  (S(x)  +  Yr(x/]  iji(x)  =  0  , 


where  r  (x)  =  0 

p(x)  =  exp  j/dx  ~  -  Vl  -  Dd 

c  d 

S(x)  =  p(x)  c  c  ,  -  Dc  +  -jii  Dd  -  ~ — •  (d 
W  p— 1  U  d  u  b  u 

M  P-1 

'I'(X)  -  Ry1)(x)  ‘ 

If  either  Eq.  (A2a)  or  (A2b)  Is  a  two-term  recurrence  relation  (i.e.,  if 
bu  or  dy  is  equal  to  zero  for  all  u) ,  then  the  above  expressions  are  undefined 
and  does  not  satisfy  the  Sturra-Liouvll le  differential  equation.  In  this 

case  Ry*^  satisfies  instead  a  first-order  differential  equation  and  is  in  the 
form  of  an  exponential.  This  may  be  readily  seen  by  letting  dy  =  0  in 
Eq.  (A2b),  in  which  case 


R<l)(x)  -  exp[/dxcw!  .  (A5) 

On  the  other  hand,  if  by  =*  0,  we  obtain  from  Eq.  (A2a):  R^*^  “  * 

which  when  combined  with  Eq.  (A2b)  yields 
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(A6) 


Rp1J(x)  =  expj/dxjc^  + 

An  extensive  search  of  the  literature  indicated  that  functions  satisfying 
the  recurrence  relations  [Eq.  (A2)]  have  not  previously  been  named.  For  the 
purposes  of  this  Appendix  we  shall  refer  to  them  as  birecurrent  functions. 

Most  of  the  special  functions  of  physics  fail  into  this  category  (including 
all  Bessel  functions,  Legendre  functions,  Hermite  polynomials,  etc.).  We 
exclude  the  special  cases  given  in  Eqs.  (A5)  and  (A6)  from  this  category, 
preferring  to  call  them  exponential  terms  instead. 

At  this  point  it  will  be  noted  that  the  assumed  form  [Eq.  (Al)]  may  be 
applied  to  the  integrals  occurring  in  Eq.  (115)  with  the  following 
Identifications  made 


f(x)+r' 


♦  4m>(k„r') 


» 


where  zj^m^(knr ’ )  is  the  2th  order  Bessel  function  of  the  mth  kind.  For  m  = 

1,  Z £  may  be  taken  to  be  J £,  the  Bessel  function  of  the  first  kind.  For  in  = 

2,  Zjj  may  be  taken  to  be  either  or  the  Hankel  functions  of  the 

first  and  second  kinds,  respectively.  The  quantity  kn  is  taken  to  be  any  of 
the  appropriate  wavenumbers  occurring  in  Eq.  (115).  Hence,  the  technique 
herein  described  could  potentially  be  applied  directly  to  the  integrals 
occurring  in  this  equation  (as  will  be  noted  later,  this  turns  out  to  be 
somewhat  Impractical). 

The  Integration  technique  presented  in  this  Appendix  involves  a  general¬ 
ization  of  the  method  (described  by  Watson  [68])  used  by  Sonlne  [b7]  to  evalu¬ 
ate  certain  indefinite  Integrals  of  Bessel  functions.  The  integral  to  be 
evaluated  in  Sonine’s  method  is  replaced  by  a  differential  equation.  A  par¬ 
ticular  solution  of  the  differential  equation  is  then  sufficient  to  express 
the  result  of  integration.  In  the  present  work  we  generalize  the  method  to 
Include  all  functions  obeying  the  relations  [Eq.  (A2)].  In  addition  we 
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describe  an  approach  for  obtaining  and  solving  the  appropriate  differential 
equations . 


2 .  The  Details. 

We  assume  the  integral  (Eq.  (Al)]  may  be  expressed  in  the  form. 


1  1 
l  •••  l 


P1=0  P2=0 


p  **0 

ra 


P1,P2 


(x)  n  R'  , 

•**p-  i-l  Vpl 


(A7 ) 


1  m 


where  the  2m  coefficients  A  (x)  are  functions  to  be  determined. 

pl,p2 . Pm 

For  convenience,  we  will  represent  the  multiple  summation  and  the  coefficients 


in  Eq.  (A7)  by  the  shorthand  notations  £  and  A  ,  respectively.  In  order  to 

{p>  P 

determine  the  functions  Ap,  we  differentiate  Eq.  (A7),  substitute  for  DI  the 
integrand  from  Eq.  (Al),  and  obtain 

f(x)  n  r(1)  =  l  [a  d  n  r(1|  +  (da  )  n  R  1  • 

i-l  “i  (pi  L  p  i-i  Vi  »  i-i  ui+pJ  <A8> 


Due  to  the  recurrence  relations  (Eq.  (A2)],  it  is  always  possible  to 
express  the  first  sum  on  the  right-hand  side  of  Eq.  (A8)  in  the  form 
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(A9) 


where  the  2^™  coefficients  BM  =  B_  _  „  „  (x)  are  known 

Pq  Pi »P2 » • • • »Pmqi »q2 * • • • »qm 

functions  resulting  from  repeated  applications  of  the  relations  (Eq.  (A2)J  and 

m  ( i) 

the  regrouping  of  terms  in  the  form  n  R  .  . 

i-l  Ui  pi 
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Using  Kq.  (A9)  we  can  rewrite  Eq .  (A8)  to  obtain 


f(x) 
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n 

1=1 
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(A10) 


We  can  now  obtain  a  coupled  set  of  differential  equations  for  the 
functions  Ap  by  Imposing  the  sufficient  condition  that  the  coefficients  of 
like  special  functions  on  each  side  of  Eq.  (A10)  be  equal.  Doing  this,  we 
obtain  the  following  coupled  set  of  linear  inhomogeneous  differential 
equations  of  first  order 


m 

f(x)  n 
j=i 


DA  + 
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l 

(q) 


B  A  , 

pq  q 


(All) 


where  6  is  the  Kronecker  delta. 

In  solving  the  set  [Eq.  (All)]  of  2®  equations  in  the  2m  unknown 
functions  Ap,  one  normally  proceeds  by  differentiation  and  algebraic 
manipulation  to  uncouple  a  particular  function  from  the  remainder.  This 
results  in  a  differential  equation  of  order  2m+^.  A  particular  solution  of 
this  uncoupled  equation  involves  a  particular  choice  of  2m+^  constants.  Since 
this  is  exactly  the  number  of  arbitrary  constants  that  the  original  set 
[Eq.  (All)]  involves,  one  must  be  careful  not  to  introduce  any  further 
arbitrary  constants.  In  this  case  we  obtain  the  remaining  functions  by 
expressing  them  in  terms  of  derivatives  of  the  initial  function  that  has  been 
calculated,  rather  than  in  terms  of  integrals  of  it.  Regardless  of  the  method 
used  in  obtaining  a  particular  solution  of  Eq.  (All),  one  must  avoid 
Introducing  more  than  2m+^  arbitrary  constants.  Otherwise,  the  solution  so 
obtained  will  neither  satisfy  Eq.  (All)  nor  provide,  via  Eq.  (A7),  a  proper 
representation  of  the  integral  Eq.  (Al). 

When  the  integrand  of  Eq.  (Al)  contains  more  than  one  blrecurrent 
function,  it  may  be  desirable  to  move  one  (or  possibly  more)  of  the 
blrecurrent  functions  out  of  the  product  term  and  treat  it  as  part  of  f(x). 
Each  blrecurrent  function  appearing  in  the  product  term  doubles  the  number  of 
unknown  coefficient  functions  Ap  and,  hence,  doubles  the  number  of  coupled 
differential  equations  to  be  solved.  Thus,  we  halve  the  number  of 

1  A  7 
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differential  equations  each  time  we  move  a  birecurrent  function  out  of  the 
product  term  and  group  It  with  f(x).  However,  each  function  so  grouped  will 
appear  In  the  Inhomogeneous  term  of  the  final  set  of  differential  equations. 
Conversely,  none  of  the  blrecurrent  functions  grouped  In  the  product  term  will 
appear  explicitly  in  the  final  set  of  differential  equations. 

Any  particular  solution  of  Eq.  (All)  will  give  a  set  of  functions  Ap  that 
can  be  used  in  Eq.  (A7)  to  express  the  result  of  the  integration.  We  can  see 
this  if  we  differentiate  the  expression  that  results  by  substitution  of  this 
particular  set  into  Eq.  (A7).  The  resulting  Eq.  (A8)  is  obviously  satisfied 
since  the  Ap  are  a  particular  solution.  The  fact  that  only  a  particular, 
rather  than  a  general,  solution  is  required  is  a  powerful  aspect  of  the 
technique. 

The  coupled  set  [Eq.  (All)]  is  a  standard  form  of  linear  inhomogeneous 
differential  equations  of  first  order  that  may  be  solved  by  well-known 
methods.  A  particular  solution  of  Eq.  (AH)  is  easier  to  obtain  than  one  may 
suspect  since  each  equation  contains  exactly  one  term  involving  the  derivative 
of  a  particular  function  Ap  and  the  derivative  of  each  of  the  functions  Ap 
appears  in  only  one  equation. 

The  technique  described  above  is  equally  applicable,  with  slight 
modifications,  when  one  or  more  of  the  birecurrent  functions  in  the  product 
term  of  Eq.  (Al)  is  replaced  by  an  exponential  term  of  the  form  Eq.  (A5)  or 
(A6).  Because  the  exponential  terms  satisfy  two-term  recurrence  relations,  we 
do  not  have  to  include  p  **  1  terms  for  them  in  Eq.  (A7).  This  reduces  the 
number  of  unknown  coefficients  Ap  and  the  resulting  coupled  differential 
equations  by  a  factor  of  2.  Moving  the  exponential  term  from  the  product  term 
into  f(x)  does  not  change  the  number  of  differential  equations  to  be  solved. 

3.  An  Example. 

To  illustrate  the  technique,  we  obtain  the  result  to  the  following  well- 
known  integral:  I  *  /dxxsinpx.  In  this  case,  f(x)  -  x  and  Rp(x)  -  sinpx. 

The  simplest  way  to  apply  the  technique  to  this  problem  is  to  consider 
sinpx  as  the  imaginary  part  of  the  exponential  Rp(x)  “  exp(ipx)  and  assume 
that  I  -  A(x)exp( lux) .  Only  one  term,  And  hence  only  one  unknown  coefficient 
A(x),  is  required  in  l  in  this  case  because  differentiation  of  the  exponential 
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does  not  produce  new  functions  (i.e.,  recurrence  relation  [Eq.  (A2b)]  reduces 
to  a  2-term  equation  relating  DR^  to  Ry). 

We  shall  use  Instead  a  somewhat  more  complicated  approach  that  requires 
two  unknown  coefficients  in  order  to  illustrate  several  important  aspects  of 
the  technique.  This  approach  is  based  on  the  fact  that  the  set  of  two 
functions  sinux  and  cosux  is  closed  under  differentiation  so  that  it  is 
convenient  to  choose  =■  sinux  and  R,j+i  *  cosux.  We  now  proceed  to  obtain 
the  integral  following  a  step-by-step  procedure: 

a.  We  assume  I  may  be  expressed  in  the  form 

I  *  Ag(x)sinux  +  A^(x)cosux  .  (A12) 

b.  Differentiation  produces  DI  *  mAq(x)cos|ix  +  [DAg(x) ]sinux 
-  uA^(x)sinux  +  [DAj^x) ]cosux. 

c.  Equating  Dl  to  the  integrand  xsinpx  and  separately  equating  coeffi¬ 
cients  of  sinux  and  cosux,  we  obtain  the  following  differential  equations 

DAq  -  uA^  =*  x,  (A13a) 

UAQ  +  DA^  -  0.  (A13b) 

d.  We  now  uncouple  Ag  and  A^  by  substituting  into  Eq.  (Al3a)  the 
expression  for  Ag  obtained  from  Eq.  (A13b).  This  gives 

D2A1  +  u2At  =  -ux.  (A14) 

e.  A  particular  solution  of  Eq.  (A14)  is 

A^(x)  ■  -x/u.  (A1S) 

Substitution  of  Eq.  (A15)  into  Eq.  (A13b)  yields 

AQ(x)  -  1/u2.  (Alb ) 
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f.  Substitution  of  Eq.  (A15)  and  Eq.  (A16)  into  Eq.  (A12)  gives  the 
result  /dxxsinpx  ■  (sinpx  -  pxcospx)/p2. 


tn  order  to  investigate  the  consequences  of  using  a  different  particular 
solution  than  the  one  chosen,  we  first  obtain  the  general  solution  of  Eq. 
(A14).  It  is 


A^(x)  -  CjSinux  +  C^cospx  -  x/p,  (A17) 

where  and  C2  are  arbitrary  constants.  Substitution  of  Eq.  (A17)  into  Eq. 
(A13b)  yields 


Aq(x)  *  -Cjcospx  +  C^sinpx  +  1/p2.  (A18) 

Substitution  of  these  general  solutions  for  Ag  and  Aj  into  Eq.  (A12)  produces 
the  following  expression  for  the  desired  integral: 

/dxxsiniix  =*  (sinpx  -  pxcospx)/p2  +  C 2  .  (A19) 

Equation  (A19)  involves  only  a  single  arbitrary  constant  to  which  the 
indefinite  integral  under  consideration  is  entitled.  Since  the  completely 
general  solution  to  Eq.  (A14)  was  used  in  obtaining  Eq.  (A19),  it  is  clear 
that  any  particular  solution  to  Eq.  (A14)  would  have  sufficed,  with  only  the 
constant  C2  *n  Eq.  (A19)  being  affected  by  a  different  choice. 

4 .  A  Second  Example. 

Although  the  first  example  provides  a  succinct  illustration  of  the 
present  integration  technique,  an  integration-by-parts  approach  would 
certainly  have  been  more  straightforward. 

Another  example  is  presented  that  is  less  susceptible  to  standard 
techniques  and  is  more  relevant  to  the  nonlinear  scattering  problem  at  hand. 

We  note  again  that  the  integrals  contained  in  Eq.  (115)  may  all  be  put 
Into  the  form  Eq.  (Al).  Unfortunately,  if  this  were  done  directty,  it  would 
result  in  a  set  of  eight  coupled  first-order  differential  equations. 


Uncoupling  of  this  equation  has  not,  as  yet,  proved  to  be  a  tractable 
problem.  However,  if  the  Zeroth-order  Hankel  function  (which  occurs  in  all  of 
the  integrals)  and  one  of  the  other  of  the  two  remaining  Bessel  functions  are 
replaced  by  their  asymptotic  forms  (for  large  argument  as  compared  to  order), 
the  following  general  type  of  Integral  is  obtained 

I  -  /drrMZv(r)exp(ir),  (A20) 

where  Zv(r)  is  an  arbitrary  Bessel  function  of  one  of  the  first  three  kinds  of 
real  argument  r,  and  the  range  of  integration  is  restricted  to  r>0.  For 
generality,  both  the  order  v  and  the  exponent  y  are  chosen  to  be  complex. 

To  begin  the  technique,  we  assume  that  I  may  be  written  as 

I  -  AQ(r)Zv(r)  +  A1(r)Zy+1(r)  .  (A21) 

Differentiating  Eq.  (A21),  expressing  the  resulting  Bessel  function  deriva¬ 
tives  in  terms  of  Zy  and  Zv+1  by  use  of  the  appropriate  recurrence  relations, 
equating  the  result  to  the  integrand  of  Eq.  (A20),  and  imposing  the  sufficient 
eondltlor'.  that  coefficients  of  like-order  Bessel  functions  be  equal,  we  obtain 
the  following  coupled  set  of  differential  equations 

AQ  -  DAt  +  [(v  +  l)/r]A1  -  0, 

DAq  +  (v/r)AQ  +  Al  -  rwexp(ir).  (A22) 

These  equations  may  be  uncoupled  to  yield 

D2Al  -  (l/r)DAL  +  (1  +  (l  -  v2)/r2JA1  =  rMexp(ir).  (A23) 

We  now  define  a  function  p(r)  such  that 

Ar(r)  -  rp( r ) .  (A2 4) 

Substitution  Into  Eq.  (A23)  yields 
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(A25) 


r2D2p  +  rDp  +  [r2  -  v2]p  -  rM+lexp(ir). 

One  particular  solution  to  Eq.  (A25)  is  an  associated  Bessel  function  as 
defined  by  Luke  169],  namely,  p(r)  =>  iu+lH„  ..(-ir),  so  that 

Aj(r )  -  i^rH^-ir).  (A26) 

Use  of  the  known  properties  of  the  associated  Bessel  functions  results  in 

A0(r)  =  iM+1r  [(p  -  v)(n  -  v  -  1>Hu_1>NH.1(“ir) 

+  (-ir)wexp(ir)]/(2ii  +  1),  (A27) 

when  the  expression  for  Aj  given  by  Eq.  (A26)  is  substituted  into  the  first  of 
the  Eqs.  (A22).  We  thus  have 

/drrPZv(r)exp(i.r)  »  iP-1r{[(p  -  v)(u  -  v  -  l)Htj_1  ^(-ir) 

+  (~ir )yexp( ir ) ] /( 2p  +  l)}Zv(r) 

+  i^^rH^  ^v(-ir)Zv+1(r) .  (A28) 

This  result  is  identical  to  that  obtained  by  Luke  [73]  using  a  specialized 
Integration  technique  developed  by  McLachlan  and  Meyers  [74]  for  certain  inte¬ 
grals  involving  Bessel  and  Struve  functions.  Luke  [73]  provides  formulas  by 
which  the  associated  Bessel  functions  appearing  in  Eq.  (A28)  may  be  evaluated. 

The  above  example  resulted  in  a  differential  equation  that  was 
recognizable.  However,  the  technique  is  still  applicable  even  if  no 
previously  known  solution  to  the  differential  equation  exists.  We  first  try 
to  obtain  a  particular  solution  to  our  inhomogeneous  differential  equation  by 
using  standard  methods  [75]  such  as  the  method  of  Lagrange  or  the  method  by 
Cauchy.  If  none  of  these  methods  proves  satisfactory,  we  can  always  obtain  a 
solution  in  the  form  of  an  infinite  series.  As  an  example  of  this,  we  again 
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return  to  Eq.  (A23)  and  assume  that  may  be  expressed  as 


A 


1 


exp(ir)  l  Bmrm+W, 

B«-°» 


(A29) 


where  the  coefficients  Bm  are  constants  to  be  determined.  To  simplify 
subsequent  calculations,  we  chose  the  form  of  the  expansion  to  be  compatible 
with  the  inhomogeneous  term.  If  Eq.  (A29)  is  substituted  into  Eq.  (A23)  and 
coefficients  of  like  powers  of  r  are  equated,  the  following  recursion  relation 
is  obtained  for  the  Bm: 

(m  +  P  +  v)(m  +  p  -  v)B  +  l[2(m  +  p)  -  1]B  -  6  (A30) 

mxi  m  m,l 


We  can  obtain  a  particular  solution  to  Eq.  (A30)  and,  hence,  to  Eq.  (A23) 
by  setting  Bm  “  0  for  m>l  and  solving  for  the  nonzero  coefficients  Bm,  m  *  1 , 
0,  -1,  ...,.  The  resulting  descending  power  series  representation  for  A^  can 
be  expressed  in  terms  of  a  hypergeometric  function  as 


Al  -  ip'rlr{(-l)M+1(ir)wexP(ir)3F1 

[l,-g+v,-u+v,-u-v;  (l/2)-u;  (l/2)ir]/(2p+l)} . 


(A31) 


In  view  of  Eq.  (A26),  it  is  not  surprising  that  the  quantity  in  brackets 
is  identical  to  the  series  representation  of  v(-ir)  given  by  Luke  [69], 

The  solution  [Eq.  (A31)]  is  not  defined  if  p  -  -1/2.  It  is  also  not  defined 
if  p  is  an  odd  multiple  of  1/2  unless  both  p  ±  v  are  positive  integers  or 
zero.  The  solution  is  a  terminating  series  if  either  p+vorp-visa 
positive  Integer.  The  infinite  series  obtained  otherwise  is  an  asymptotic 
representation  of  Aj  that  is  valid  for  r-**>. 

We  can  obtain  a  second  particular  solution  for  A^  by  setting  Bm  =  0  for 
nKl  and  solving  Eq.  (A30)  for  Bm,  m  *  2,3,...,.  The  hypergeometric 
representation  of  the  resulting  ascending  power  series  in  r  is 
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i“+1r 


«-ir)M+l2F2 
(1,  m+j;  m-v  +  2. 


y  +  u  +  2;  -2ir)/ 


[(y  -  v  +  i)(y  +  V  +  1)1). 


(A32) 


The  quantity  appearing  in  braces  in  Eq.  (A32)  is  the  series  representation  of 
the  associated  Bessel  function  h^  v(-ir)  as  defined  by  Luke  [69J.  This 
function  is  a  second  particular  solution  to  the  differential  equation 
satisfied  by  Hy  v  and  hence  can  be  used  in  place  of  H  in  the  solution 
[Eq.  (28)]  to  the  original  integral. 

The  solution  [Eq.  (A32)]  is  a  terminating  series  if  p  is  a  positive  odd 
multiple  of  -1/2  (other  than  -1/2)  and  if  both  y  ±  v  are  not  positive 
integers.  It  is  not  defined  if  either  y+vory-visa  negative  integer 
and  y  is  not  a  positive  odd  multiple  of  -1/2  (other  than  1/2). 

5 .  Additional  Illustrative  Examples. 

The  two  previous  examples  illustrate  the  power  and  versatility  of  the 
current  integration  technique,  but  two  objections  may  be  raised:  The  first 
example  can  be  handled  trivially,  and  the  second  example  is  one  involving 
products  of  Bessel  functions  and,  hence,  is  amenable  to  the  original  approach 
proposed  by  Sonine.  The  examples  that  follow  will  serve  to  illustrate  the 
applicability  of  the  current  technique  to  integrals  that  are  not  of  the  Bessel 
function  type.  Although  some  of  these  may  be  solved  by  standard  techniques, 
they  nonetheless  illustrate  the  broad  range  of  integrands  that  can 
successfully  be  handled  via  this  technique  (and  to  the  author's  knowledge, 
several  of  these  Integrals  have  not  been  previously  tabulated). 

a.  Some  integrals  Involving  Legendre  functions.  We  now  consider  some 

examples  of  integrals  of  the  general  form 

I  -  JdxPv(x)f(x),  (A33) 
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where  Pv(x)  is  the  Legendre  function  of  order  v  and  f(x)  has  the  same  meaning 
as  in  Eq.  (Al).  We  assume  the  integral  I  may  be  represented  in  the  usual  way 


A(x)P  (x)  +  B(x)P  ,,(x).  Following  the  procedure  outlined  in  Section  2,  we 

V  vtl 

obtalOwo  coupled  equations  for  A  and  B  which  can  be  uncoupled  to  yield 


-xB(x)  + 


(1-x2) 

(v+l) 


B*(x), 


(A34) 


where 


(1-x  )B"(x)  -  2xB’(x)  +  v(v+l)B(x)  *  (v+l)f(x). 


(A35) 


As  a  first  example  of  an  integral  of  the  form  Eq.  (A33),  we  consider  the 
case  where  f(x)  ■  1.  A  particular  solution  to  Eq.  (A35)  for  this  case  is  B(x) 
=  (1/v).  Equation  (A34)  now  determines  A(x)  to  be  A(x)  *  -(x/v).  Substitution 
of  A  and  B  into  the  representation  for  I  now  produces 


/dxPv(x)  -  -  £pv  +ipv+l*  v  *  0 


(A36) 


We  next  consider  a  more  challenging  integrand  [i.e.,  one  which  cannot  be 
handled  by  direct  manipulation  of  the  recurrence  relations  for  Pv(x)].  Let 
f(x)  be  ln(l±x).  In  obtaining  a  particular  solution  to  Eq.  (A35)  we  use  the 
inhomogeneous  term  as  a  guide  and  assume  B(x)  «*  K^ln(l±x)  +  ,  where  and 

K2  are  undetermined  constants.  Direct  substitution  into  Eq.  (A35)  gives 
*  1/v  and  K2  *  l/[v^(v+l)].  Equation  (A34)  may  next  be  used  to  show 
that  A(x)  *  {-(x/v)ln(l±x)  ±  l/[v(v+l)]  -  x/v^}.  Substitution  of  A  and  B 
into  the  representation  for  I  results  in 


/ dx  lnClixjPfx) 


ln(l±x  )  ± 


V( V+l ) 


in  -  }\  v 


iln(l±x)  + 

L  v*(v+l)J 


v  *  0,-1 


(A37) 


155 


b.  Some  integrals  involving  Hermite  functions .  We  next  consider 


Integrals  of  the  general  form 

I  »  JdxHy(x)f(x)  ,  (A38) 

where  Hv(x)  is  the  Hermite  function  of  order  v,  and  once  again  f(x)  has  the 
same  meaning  as  was  used  in  connection  with  Eq.  (Al).  We  assume  the  integral 
of  Eq.  (A38)  may  be  represented  as  A(x)Hy(x)  +  B(x)Hv_^(x).  (Note  that  Hv  and 
are  used  to  represent  the  integral  as  opposed  to  Rv  and  This 

difference  is  inconsequential.  Any  two  orders  separated  by  one  integral  value 
will  be  adequate  to  Implement  the  procedure.)  Omitting  the  details,  we 
uncouple  the  resulting  coupled  set  to  obtain 

A(x)  -  -  i  B(x)  -  (A39) 

where 


2vf(x)  =*  -  B"(x)  -  2xB'(x)  -  2(v+1)B(x). 


(A40) 


As  a  first  example  of  this  general  form,  we  let  f(x)  *  e~x  .  [This  is 
the  usual  weighting  function  used  in  the  orthogonality  integral  for  Hy(x).] 

If  we  assume  a  particular  solution  of  Eq.  (A40)  of  the  form  B(x)  =  Ke-X  (with 
K  an  undetermined  constant),  direct  substitution  gives  K  ■  -1.  Equation  (A39) 
gives  A(x)  **  0  so  that 


_  2  _  2 
/ dxe  X  Hy(x)  -  -e  X  Hvl(x), 


(A41) 


As  a  second  example  we  let  f(x)  ■  x_lJ,  where  for  the  moment  y  is  an 

arbitrary  exponent.  If  we  assume  a  particular  solution  of  the  form  B(x)  * 

-K, 

K^x  ,  where  and  K2  are  constants,  direct  substitution  into  (A42)  produces 


2vx  M  -  -K^d^+Ux 


~(K,+2)  -K 

+  2*^  -  (v+1) Jx 


(A42) 
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There  are  two  sets  of  values  for  the  constants  ,  K2 ,  and  p  for  which 

— 2v 

Eq.  (A42)  Is  satisfied:  *  (~v+l ) ( \>+2)  *  *  V+1  *  y  *  v+^  ’  an<* 

K1  ”  ~  0&2)  »  *2  ”  ’  11  *  ”*•  *  Usin8  the  first  set  of  values,  we  obtain 

the  result 

/dxHv(x)x 


- ( v+3) 


r  2x~v  x_(u+2)i 

(v+l)(v+2)  (  v+2 ) 


Hv,(x> 


2u 

(v+l)(v+2) 


x“(w+1)Hv_1(x) 


(A43) 


v  *  -1,  -2. 


We  obtain  from  the  second  set  of  values 

/dxxHv(x)  -  »v(x)  -  H^Cx)  (A44) 

v  *  -2. 

For  a  final  example  involving  Hermite  functions,  we  let  f(x)  =  xelYx, 
where  y  Is  a  constant  Initially  assumed  to  be  arbitrary. 

A  particular  solution  to  Eq.  (A40)  can  be  obtained  with  B(x)  *  KeiYx, 
where  K  is  an  unknown  constant.  Direct  substitution  shows  that  a  solution 
exists  for  y  *  /2(vfl)  and  K  *  iv//2(v+l )  .  Using  the  resulting  solution  for 
B  in  Eq.  (A39)  to  obtain  A,  we  then  have 


r.  i/2(v+l)x  u  ,  . 
Jdxxe  Hy(x) 


i/2(v+l)x 


(7=!  +  j)  Hv(x)  + 
\/2Tv+1)  / 


iv 


/2(v+l) 


H  ,  (x) 
v-lv 


V  *  -l.  (A43) 
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c.  Some  examples  Involving  Laguerre  functions*  We  now  consider 
Integrals  of  the  general  form 

I  »  /dxf (x)Ly(x)  ,  (A46) 

where  Lv(x)  is  the  Laguerre  function  of  order  v,  and  f(x)  has  the  same  meaning 
as  in  Eq.  (Al).  As  usual,  we  represent  I  In  the  form  A(x)Ly(x)  +  B(x)Lv_^(x) 
and  obtain  the  following  uncoupled  equations 

A(x)  -  -  l)  B(x)  +£b'(x),  (A47) 


where 


vf(x)  -  xB”(x)  +  (x+l)B'(x)  +  (\H-l)B(x).  (A48) 

As  an  example  of  this  case  we  let  f(x)  «  xe~(v+l)x.  To  obtain  a 
particular  solution  of  Eq.  (A48),  we  assume  B(x)  *  Ke”^v+1)x.  Direct 
substitution  gives  K  *  l/(v+l).  Using  this  value  yields  the  integral 

r  -(\rt-Ux  -(vfl)x 

/ dxxe  '  }  Ly(x)  -  —  l-(1+x>Lv(x)  +  Lv-l(x)1 

v  *  -1.  (A49) 

As  a  final  example,  we  let  f(x)  ■  x(l+x)~^  .  A  particular  solution  to 

Eq.  (A48)  can  then  be  obtained  assuming  B(x)  *  K(l+x)~^v+^ ,  where  again  K  is 
an  unknown  constant.  Direct  substitution  into  Gq.  (A48)  yields 
K  *  v/  [  ( v+l )  (  v+2 )  ]  so  that 


/dxx(l+x)  ^3^Lv(x) 


(1+x) 


-(v+1) 


(v+2) 


Mx>  + 


V  *  -1,-2. 


(A50) 
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d.  Some  final  results.  The  following  is  a  tabulation  of  some  additional 
results  obtained  using  the  Integration  technique  described  In  this  paper.  For 
the  sake  of  brevity  the  derivations  have  been  omitted. 


2 

/ dxxe'x  Hv(x)Hm(x) 


)  (M+v+1) 

1  2[(u-v)2-1] 


Hy(x)Hv(x)  +  ( 


v=S«>Hv<*>Vi<*) 


+  ( 


(A51) 


where  H  is  a  Hermite  function  and  y-  M  i  1. 

/dx[P  (X)]2  «  f  IV*)2  +  P-jjC*)2!  -  Vx)p-*s(x) 


and 

[  --  (-+1(-^1)(--+-1--]  /d*Pv(x)xU  -  /dxPv(x)x 


(A52) 


